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Selected solutions of Einstein's field equations: 
their role in general relativity and astrophysics 



Jiff Bicak 

Institute of Theoretical Physics, 
Charles University, Prague 

1 Introduction and a few excursions 

The primary purpose of all physical theory is rooted in reality, and most rela- 
tivists pretend to be physicists. We may often be members of departments of 
mathematics and our work oriented towards the mathematical aspects of Ein- 
stein's theory, but even those of us who hold a permanent position on "scri" , 
are primarily looking there for gravitational waves. Of course, the builder of 
this theory and its field equations was the physicist. Jiirgen Ehlers has always 
been very much interested in the conceptual and axiomatic foundations of 
physical theories and their rigorous, mathematically elegant formulation; but 
he has also developed and emphasized the importance of such areas of rela- 
tivity as kinetic theory, the mechanics of continuous media, thermodynamics 
and, more recently, gravitational lensing. Feynman expressed his view on the 
relation of physics to mathematics as follows [Q : 

"The physicist is always interested in the special case; he is never inter- 
ested in the general case. He is talking about something; he is not talking 
abstractly about anything. He wants to discuss the gravity law in three di- 
mensions; he never wants the arbitrary force case in n dimensions. So a certain 
amount of reducing is necessary, because the mathematicians have prepared 
these things for a wide range of problems. This is very useful, and later on it 
always turns out that the poor physicist has to come back and say, 'Excuse 
me, when you wanted to tell me about four dimensions...' " Of course, this 
is Feynman, and from 1965... 

However, physicists are still rightly impressed by special explicit formulae. 
Explicit solutions enable us to discriminate more easily between a "physical" 
and "pathological" feature. Where are there singularities? What is their char- 
acter? How do test particles and fields behave in given background space- 
times? What are their global structures? Is a solution stable and, in some 
sense, generic? Clearly, such questions have been asked not only within gen- 
eral relativity. 

By studying a special explicit solution one acquires an intuition which, 
in turn, stimulates further questions relevant to more general situations. 
Consider, for example, charged black holes as described by the Reissner- 
Nordstrom solution. We have learned that in their interior a Cauchy horizon 
exists and that the singularities are timelike. We shall discuss this in greater 
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detail in Section 3.1. The singularities can be seen by, and thus exert an in- 
fluence on, an observer travelling in their neighborhood. However, will this 
violation of the (strong) cosmic censorship persist when the black hole is 
perturbed by weak ("linear") or even strong ("nonlinear") perturbations? 
We shall see that, remarkably, this question can also be studied by explicit 
exact special model solutions. Still more surprisingly, perhaps, a similar ques- 
tion can be addressed and analyzed by means of explicit solutions describing 
completely diverse situations - the collisions of plane waves. As we shall note 
in Section 8.3, such collisions may develop Cauchy horizons and subsequent 
timelike singularities. The theory of black holes and the theory of colliding 
waves have intriguing structural similarities wliic;li, first of all, stem from the 
circumstance that in both cases there exist two symmetries, i.e. two Killing 
fields. What, however, about more general situations? This is a natural ques- 
tion inspired by the explicit solutions. Then "the poor physicists have to come 
back" to a mathematician, or today alternatively, to a numerical relativist, 
and hope that somehow they will firmly learn whether the cosmic censor- 
ship is the "truth", or that it has boon a very inspirational, but in general 
false conjecture. However, even after the formulation of a conjecture about 
a general situation inspired by particular exact solutions, newly discovered 
exact solutions can play an important role in verifying, clarifying, modifying, 
or ruling out the conjecture. And also "old" solutions may turn out to act 
as asymptotic states of general classes of models, and so become still more 
significant. 

Exact explicit solutions have played a crucial role in the development of 
many areas of physics and astrophysics. Later on in this Introduction we 
will take note of some general features which are specific to the solutions of 
Einstein's equations. Before that, however, for illustration and comparison 
we shall indicate briefly with a few examples what influence exact explicit 
solutions have had in other physical theories. Our next introductory excur- 
sion, in Section 1.2, describes in some detail the (especially early) history of 
Einstein's route to the gravitational fleld equations for which his short stay in 
Prague was of great significance. The role of Ernst Mach (who spent 28 years 
in Prague before Einstein) in the construction of the first modern cosmolog- 
ical model, the Einstein static universe, is also touched upon. Section 1.3 is 
devoted to a few remarks on some old and new impacts of the other simplest 
"cosmological" solutions of Einstein's equations - the Minkowski, the de Sit- 
ter, and the anti de Sitter spacetimes. Some specific features of solutions in 
Einstein's theory, such as the observability and interpretation of metrics, the 
role of general covariance, the problem of the equivalence of two metrics, and 
of geometrical characterization of solutions are mentioned in Section 1.4. Fi- 
nally, in the last (sub) sections of the "Introduction" we give some reasons 
why we consider our choice of solutions to be "a natural selection" , and we 
briefly outline the main body of the article. 
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1.1 A word on the role of explicit solutions in other parts of 
physics and astrophysics 

Even in a linear theory like Maxwell's electrodynamics one needs a good 
sample, a useful kit, of exact fields like the homogeneous field, the Coulomb 
field, the dipole, the quadrupole and other simple solutions, in order to gain 
a physical intuition and understanding of the theory. Similarly, of course, 
with the linearized theory of gravity. Going over to the Schrodinger equa- 
tion of standard quantum mechanics, again a linear theory, consider what we 
have learned from simple, explicitly soluble problems like the linear and the 
three-dimensional harmonic oscillator, or particles in potential wells of var- 
ious shapes. We have acquired, for example, a transparent insight into such 
basic quantum phenomena as the existence of minimum energy states whose 
energy is not zero, and their associated wave functions which have a certain 
spatial extent, in contrast to classical mechanics. The three-dimensional prob- 
lems have taught us, among other things, about the degeneracy of the energy 
levels. The case of the harmonic oscillator is, of course, very exceptional since 
Hamiltonians of the same type appear in all problems involving quantized os- 
cillations. One encounters them in quantum electrodynamics, quantum field 
theory, and likewise in the theory of molecular and crystalline vibrations. It is 
thus perhaps not so surprising that the Hamiltonian and the wave functions 
of the harmonic oscillator arise even in the minisuperspace models associated 
with the Hartle-Hawking no-boundary proposal for the wave function of the 
universe (2|, and in the minisuperspace model of homogeneous spherically 
symmetric dust filled universes j^. 

In nonlinear problems explicit solutions play still a greater role since to 
gain an intuition of nonlinear phenomena is hard. Landau and Lifshitz in 
their Fluid Mechanics (Volume 6 of their course) devote a whole section to 
the exact solutions of the nonlinear Navier-Stokes equations for a viscous 
fluid (including Landau's own solution for a jet emerging from the end of a 
narrow tube into an infinite space filled with fluid) . 

Although Poisson's equation for the gravitational potential in the classical 
theory of gravity is linear, the combined system of equations describing both 
the field and its fluid sources (not rigid bodies, these are simple!) character- 
ized by Euler's equations and an equation of state are nonlinear. In classical 
astrophysical fluid dynamics perhaps the most distinct and fortunate example 
of the role of explicit solutions is given by the exact descriptions of ellipsoidal, 
uniform density masses of self-gravitating fluids. These "ellipsoidal figures of 
equilibrium" Q include the familiar Maclaurin spheroids and triaxial Jacobi 
ellipsoids, which are characterized by rigid rotation, and a wider class dis- 
covered by Dedekind and Riemann, in which a motion of uniform vorticity 
exists, even in a frame in which the ellipsoidal surface is at rest. The solutions 
representing the rotating ellipsoids did not only play an inspirational role in 
developing basic concepts of the theory of rigidly rotating stars, but quite un- 
expectedly in the study of inviscid, differentially rotating polytropes. These 
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closely resemble Maclaurin spheroids, although they do not maintain rigid 
rotation. As noted in the well-known monograph on rotating stars "the 
classical work on uniformly rotating, homogeneous spheroids has a range of 
validity much greater than was usually anticipated" . It also influenced galac- 
tic dynamics Q: the existence of Jacobi ellipsoids suggested that a rapidly 
rotating galaxy may not remain axisymmetric, and the Riemann ellipsoids 
demonstrated that there is a distinction between the rate at which the matter 
in a triaxial rotating body streams and the rate at which the figure of the 
body rotates. Since rotating incompressible ellipsoids adequately illustrate 
the general feature of rotating axisymmetric bodies, they are also used in 
the studies of double stars whose components are close to each other. The 
disturbances caused by a neighbouring component are treated as first order 
perturbations. Relativistic effects on the rotating incompressible ellipsoids 
have been investigated in the post-Newtonian approximation by various au- 
thors, recently with a motivation to understand the coalescence of binary 
neutron stars near their innermost stable circular orbit (see |^] for the latest 
work and a number of references). 

As for the last subject, which has a more direct connection with exact 
explicit solutions of Einstein's equations, we want to say a few words about 
integrable systems and their soliton solutions. Soliton theory has been one 
of the most interesting developments in the past decades both in physics 
and mathematics, and gravity has played a role both in its birth and recent 
developments. It has been known from the end of the last century that the 
celebrated Korteweg-de Vries nonlinear evolution equation, which governs 
one dimensional surface gravity waves propagating in a shallow channel of 
water, admits solitary wave solutions. However, it was not until Zabusky and 
Kruskal (the Kruskal of section 2.4 below) did extensive numerical studies 
of this equation in 1965 that the remarkable properties of the solitary waves 
were discovered: the nonlinear solitary waves, named solitons by Zabusky 
and Kruskal, can interact and then continue, preserving their shapes and 
velocities. This discovery has stimulated extensive studies of other nonlin- 
ear equations, the inverse scattering methods of their solution, the proof 
of the existence of an infinite number of conservation laws associated with 
such equations, and the construction of explicit solutions (see ^] for a re- 
cent comprehensive treatise). Various other nonlinear equations, similar to 
the sine-Gordon equation or the nonlinear Schrodinger equation, arising for 
example in plasma physics, solid state physics, and nonlinear optics, have 
also been successfully tackled by these methods. At the end of the 1970s sev- 
eral authors discovered that Einstein's vacuum equations for axisymmetric 
stationary systems can be solved by means of the inverse scattering meth- 
ods, and it soon became clear that one can employ them also in situations 
when both Killing vectors are spacelike (producing, for example, soliton-type 
cosmological gravitational waves). Dieter Maison, one of the pioneers in ap- 
plying these techniques in general relativity, describes the subject thoroughly 
in this volume. We shall briefly meet the soliton methods when we discuss 
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the uniformly rotating disk solution of Neugebauer and Meinel (Section 6.3), 
colliding plane waves (Section 8.3), and inhomogeneous cosmological models 
(Section 12.2). Our aim, however, is to imderstand the meaning of solutions, 
rather than generation techniques of finding them. From this viewpoint it 
is perhaps first worth noting the interplay between numerical and analytic 
studies of the soliton solutions - hopefully, a good example of an interaction 
for numerical and mathematical relativists. However, the explicit solutions 
of integrable models have played important roles in various other contexts. 
The most interesting muZti-dimensional integrable equations are the four- 
dimensional self-dual Yang-Mills equations arising in field theory. Their so- 
lutions, discovered by R. Ward using twistor theory, on one hand stimulated 
Donaldson's most remarkable work on inequivalent differential structures on 
four-manifolds. On the other hand, Ward indicated that many of the known 
integrable systems can be obtained by dimensional reduction from the self- 
dual Yang-Mills equations. Very recently this view has been substantiated 
in the monograph by Mason and Woodhouse The words by which these 
authors finely express the significance of exact solutions in integrable systems 
can be equally well used for solutions of Einstein's equations: "they combine 
tractability with nonlinearity, so they make it possible to explore nonlinear 
phenomena while working with explicit solutions" .|] 

1.2 Einstein's field equations 

Since Jiirgen Ehlers has always been, among other things, interested in the 
history of science, he will hopefully tolerate a few remarks on the early his- 
tory of Einstein's equations to which not much attention has been paid in 
the literature. It was during his stay in Prague in 1911 and 1912 that Ein- 
stein's intensive interest in quantum theory diminished, and his systematic 
effort in constructing a relativistic theory of gravitation began. In his first 
"Prague theory of gravity" he assumed that gravity can be described by a 
single function - the local velocity of light. This assumption led to insur- 
mountable difficulties. However, Einstein learned much in Prague on his way 
to general relativity |]ll| : he understood the local significance of the principle 
of equivalence; he realized that the equations describing the gravitational field 
must be nonlinear and have a form invariant with respect to a larger group 

^ In 1998, in the discussion after his Prague lecture on the present role of physics 
in mathematics. Prof. Michael Atiyah expressed a similar view that even with 
more powerful supercomputers and with a growing body of general mathematical 
results on the existence and uniqueness of solutions of differential equations, the 
exact, explicit solutions of nonlinear equations will not cease to play a significant 
role. (As it is well known. Sir Michael Atiyah has made fundamental contributions 
to various branches of mathematics and mathematical physics, among others, to 
the theory of solitons, instantons, and to the twistor theory of Sir Roger Penrose, 
with whom he has been interacting "under the same roof" in Oxford for 17 years 

00 
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of transformations than the Lorentz group; and he found that ^^spacetime 
coordinates lose their simple physical meaning" , i.e. they do not determine 
directly the distances between spacetime points]^ In his "Autobiographical 
Notes" Einstein says: "Why were seven years ... required for the construction 
of general relativity? The main reason lies in the fact that it is not easy to 
free oneself from the idea that coordinates must have an immediate metrical 
meaning"... Either from Georg Pick while still in Prague, or from Marcel 
Grossmann during the autumn of 1912 after his return to Zurich (cf. 
Einstein learned that an appropriate mathematical formalism for his new the- 
ory of gravity was available in the work of Riemann, Ricci, and Levi-Civita. 
Several months after his departure from Prague and his collaboration with 
Grossmann, Einstein had general relativity almost in hand. Their work [p^ 
was already based on the generally invariant line element 

ds^ — g^i/dx'^dx'^ (I) 

in which the spacetime metric tensor (7^^(x''), ^^v, p ~ 0, 1, 2, 3, plays a dual 
role: on the one hand it determines the spacetime geometry, on the other it 
represents the (ten components of the) gravitational potential and is thus a 
dynamical variable. The disparity between geometry and physics, criticized 
notably by Ernst Mach£]had thus been removed. When searching for the field 
equations for the metric tensor, Einstein and Grossmann had already real- 
ized that a natural candidate for generally covariant field equations would 
be the equations relating - in present-day terminology - the Ricci tensor 
and the energy-momentum tensor of matter. However, they erroneously con- 
cluded that such equations would not yield the Poisson equation of Newton's 
theory of gravitation as a first approximation for weak gravitational fields 
(see both §5 in the "Physical part" in |l3| written by Einstein and §4, be- 
low equation (^6|), in the "Mathematical part" by M. Grossmann). Einstein 
then rejected the general covariance. In a subsequent paper with Grossmann 
[ p^ , they supported this mis-step by a well-known "hole" meta-argument and 
obtained (in today's terminology) four gauge conditions such that the field 
equations were covariant only with respect to transformations of coordinates 
permitted by the gauge conditions. We refer to, for example, [|l^ for more 
detailed information on the further developments leading to the final version 
of the field equations. Let us only summarize that in late 1915 Einstein first 
readopted the generally covariant field equations from 1913, in which the 
Ricci tensor R^^, was, up to the gravitational coupling constant, equal to the 
energy-momentum tensor T^^ (paper submitted to the Prussian Academy on 

^ At that time Einstein's view on the future theory of gravity are best summarized 
in his reply to M. Abraham ||l^, written just before departure from Prague. 

^ Mach spent 28 years as Professor of Experimental Physics in Prague, until 1895, 
when he took the History and Theory of Inductive Natural Sciences chair in 
Vienna. 
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November 4) . From his vacuum field equations 

= 0, (II) 

where R^i, depends nonhnearly on gaf} and its first derivatives, and linearly 
on its second derivatives, he was able to explain the anomalous part of the 
perihelion precession of Mercury - in the note presented to the Academy 
on November 18. And finally, in the paper jl^ submitted on November 25 
(published on December 2, 1915), the final version of the gravitational field 
equations, or Einstein's field equations appeared^ 

1 SttG 
Rp.f - 2^iivR = -^^T^v, (III) 

where the scalar curvature R = g^^Rfiv Newton's gravitational constant G 
and the velocity of light c are the (only) fundamental constants appearing 
in the theory. If not stated otherwise, in this article we use the geometrized 



units in which G = c — 1 , and the same conventions as in 1 18 and [g^ . 

Now it is well known that Einstein further generalized his field equations 
by adding a cosmological term +Ag^,y on the left side of the field equations 
(III). The cosmological constant A appeared first in Einstein's work "Cos- 
mological considerations in the General Theory of Relativity" |Q submitted 
on February 8, 1917 and published on February 15, 1917, which contained 
the closed static model of the Universe (the Einstein static universe) - an 
exact solution of equations (HI) with yl > and an energy-momentum ten- 
sor of incoherent matter ("dust"). This solution marked the birth of modern 
cosmology. 

We do not wish to embark upon the question of the role that Mach's prin- 
ciple played in Einstein's thinking when constructing general relativity, or 
upon the intriguing issues relating to aspects of Mach's principle in present- 
day relativity and cosmologj?^ - a problem which in any event would far 
exceed the scope of this article. Although it would not be inappropriate to 

* David Hilbert submitted his paper on these field equations five days before 
Einstein, though it was published only on March 31, 1916. Recent analysis 
of archival materials has revealed that Hilbert made significant changes in the 
proofs. The originally submitted version of his paper contained the theory which 
is not generally covariant, and the paper did not include equations (III). 

^ It was primarily Einstein's recognition of the role of Mach's ideas in his route 
towards general relativity, and in his christening them by the name "Mach's 
principle" (though Schlick used this term in a vague sense three years before 
Einstein), that makes Mach's Principle influential even today. After the 1988 
Prague conference on Ernst Mach and his influence on the development of physics 
PH , the 1993 conference devoted exclusively to Mach's principle was held in 
Tiibingen, from which a remarkably thorough volume was prepared [^, covering 
all aspects of Mach's principle and recording carefully all discussion. The clarity 
of ideas and insights of Jiirgen Ehlers contributed much to both conferences and 
their proceedings. For a brief more recent survey of various aspects of Mach's 
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include it here since exact solutions (such as Godel's universe or Ozsvath's 
and Schiicking's closed model) have played a prominent role in this context. 
However, it should be at least stated that Einstein originally invented the idea 
of a closed space in order to eliminate boundary conditions at spatial infinity. 
The boundary conditions "flat at infinity" bring with them an incrtial frame 
unrelated to the mass-energy content of the space, and Einstein, in accordance 
with Mach's views, believed that merely mass-energy can influence inertia. 
Field equations (111) are not inconsistent with this idea, but they admit as 
the simplest solution an empty flat Minkowski space (T^i/ = 0, g^i, = rj^^i, — 
diag (—1, -1-1, -f 1, -1-1)), so some restrictive boundary conditions are essential 
if the idea is to be maintained. Hence, Einstein introduced the cosmological 
constant A, hoping that with this space will always be closed, and the bound- 
ary conditions eliminated. But it was also in 1917 when de Sitter discovered 
the solution |25| of the vacuum field equations (II) with added cosmological 
term {A > 0) which demonstrated that a nonvanishing A does not necessarily 
imply a nonvanishing mass-energy content of the universe. 

1.3 "Just so" notes on the simplest solutions: the Minkowski, de 
Sitter, and anti de Sitter spacetimes 

Our brief intermezzo on the cosmological constant brought up three explicit 
simple exact solutions of Einstein's field equations - the Minkowski, Einstein, 
and de Sitter spacetimes. To these also belongs the anti de Sitter spacetime, 
corresponding to a negative A. The de Sitter spacetime has the topology x 

(with corresponding to the time) and is best represented geometrically 
as the 4-dimensional hyperboloid — + w'^ + x"^ + y'^ + z"^ = {3/ A) in 5- 
dimensional flat space with metric ds^ = —dv^ + dw^ + dx^ -f dy^ + dz^ . The 
anti de Sitter spacetime has the topology x i?-^, and can be visualized as 
the 4-dimensional hyperboloid -U'^ -V"^ + X'^ + Y'^ + Z"^ = (-3/yl), yl < 0, in 
flat 5-dimensional space with metric ds"^ = —dU^ — dV'^+dX^+dY^+dZ'^. As 
is usual (cf. e.g. |2^,^), we mean by "anti de Sitter spacetime" the universal 
covering space which contains no closed timelike lines; this is obtained by 
unwrapping the circle S^. 

These spacetimes will not be discussed in the following sections. Ocassion- 
ally, for instance, in Sections 5 and 10, we shall consider spacetimes which 
become asymptotically de Sitter. However, since these solutions have played 
a crucial role in many issues in general relativity and cosmology, and most 
recently, they have become important prerequisites on the stage of the theo- 
retical physics of the "new age" , including string theory and string cosmology, 
we shall make a few comments on these solutions here, and give some refer- 
ences to recent literature. 



principle in general relativity, see the introductory section in the work [g3|, in 
which Mach's principle is analyzed in the context of perturbed Robertson- Walker 
universes. Most recently, Mach's principle seems to enter even into M theory . 
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The basic geometrical properties of these spaces are analyzed in the Bat- 
telle Recontres lectures by Penrose [^^, and in the monograph by Hawking 
and Ellis [ p6| , where also references to older literature can be found. The im- 
portant role of the de Sitter solution in the theory of the expanding universe 
is finely described in the book by Peebles and in much greater detail 
in the proceedings of the Bologna 1988 meeting on the history of modern 
cosmology p9[ |. 

The Minkowski, de Sitter and anti de Sitter spacetimes are the simplest so- 
lutions in the sense that their metrics are of constant (zero, positive, and neg- 
ative) curvature. They admit the same number (ten) of independent Killing 
vectors, but the interpretations of corresponding symmetries differ for each 
spacetime. Together with the Einstein static universe, they all are confor- 
mally flat, and can be represented as portions of the Einstein static universe 
[ p6|]27|| . However, their conformal structure is globally different. In Minkowski 
spacetime one can go to infinity along timelike geodesies and arrive to the 
future (or past) timelike infinity i~^ (or i~); along null geodesies one reaches 
the future (past) null infinity {J~ ) ; and spacelike geodesies lead to spa- 
tial infinity i^ . Minkowski spacetime can be compactified and mapped onto 
a finite region by an appropriate conformal rescaling of the metric. One thus 
obtains the well-known Penrose diagram in which the three types of infini- 
ties are mapped onto the boundaries of the compactified spacetime - see for 
example the boundaries on the "right side" in the Penrose diagram of the 
Schwarzschild-Kruskal spacetime in Fig. 3, Section 2.4, or the Penrose com- 
pactified diagram of boost-rotation symmetric spacetimes in Fig. 13, Section 
11. (The details of the conformal rescaling of the metric and resulting dia- 
grams are given in |^ 27| and in standard textbooks, for example [ p^[jl9| , p0| . ) 
In the de Sitter spacetime there are only past and future conformal infinities 
JT", both being spacelike (cf. the Penrose diagram of the "cosmological" 
Robinson- Trautman solutions in Fig. 11, Section 10); the conformal infinity 
in anti de Sitter spacetime is timelike. 

These three spacetimes of constant curvature offer many basic insights 
which have played a most important role elsewhere in relativity. To give just 
a few examples (see e.g. [^,^): both the particle (cosmological) horizons 
and the event horizons for geodesic observers are well illustrated in the de 
Sitter spacetime; the Cauchy horizons in the anti de Sitter space; and the 
simplest acceleration horizons in Minkowski space (hypersurfaces t^ = in 
Fig. 12, Section 11). With the de Sitter spacetime one learns (by considering 
different cuts through the 4-dimensional hyperboloid) that the concept of 
an "open" or "closed" universe depends upon the choice of a spacelike slice 
through the spacetime. There is perhaps no simpler way to understand that 
Einstein's field equations are of local nature, and that the spacetime topology 
is thus not given a priori, than by considering the following construction 
in Minkowski spacetime. Take the region given in the usual coordinates by 

< I7I2/I 5: Ijkl < 1, remove the rest and identify pairs of boundary 
points of the form (t, l,j/, z) and (t, — 1,?;, z), and similarly for y and z. In 
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this way the spatial sections are identified to obtain a 3-torus - a flat but 
closed manifold.^ 

The spacetimes of constant curvature have been resurrected as basic are- 
nas of new physical theories since their first appearance. After the role of the 
de Sitter universe decreased with the refutation of the steady-state cosmol- 
ogy, it has inflated again enormously in connection with the theory of early 
quasi-exponential phase of expansion of the universe, due to the false- vacuum 
state of a hypothetical scalar (inflaton) fleld(s) (see e.g. |^). We shall men- 
tion the de Sitter space as the asymptotic state of cosmological models with 
a nonvanishing A (so verifying the "cosmic no-hair conjecture") in Section 
10 on Robinson- Trautman spacetimes. Motivated by its importance in infla- 
tionary cosmologies, several new useful papers reviewing the properties of de 



Sitter spacetime have appeared 1 33 3J] ; they also contain many references to 
older literature. For the most recent work on the quantum structure of de 
Sitter space, see 

In the last two years, anti de Sitter spacetime has come to the fore in light 
of Maldacena's conjecture |Q relating string theory in (asymptotically) anti 
de Sitter space to a non-gravitational conformal field theory on the boundary 
at spatial infinity, which is timelike as mentioned above (see, e.g. |3^, where 
among others, in the Appendix various coordinate systems describing anti de 
Sitter spaces in arbitrary dimensions are discussed). 

Amazingly, the Minkowski spacetime has recently entered the active new 
area of so called pre-big bang string cosmology |Q. String theory is here 
applied to the problem of the big bang. The idea is to start from a simple 
Minkowski space (as an "asymptotic past triviality" ) and to show that it is in 
an unstable false-vacuum state, which leads to a long pre-big bangian infia- 
tionary phase. This, at later times, should provide a hot big bang. Although 
such a scenario has been criticized on various grounds, it has attractive fea- 
tures, and most importantly, can be probed through its observable relics |]38| . 

Since it is hard to forecast how the roles of these three spacetimes of 
constant curvature will develop in new and exciting theories in the next mil- 
lennium, let us better conclude our "just so" notes by stating three "stable" 
results of complicated, rigorous mathematical analyses of (the classical) Ein- 
stein's equations. 

In their recent treatise ||39|] , Christodoulou and Klainerman prove that 
any smooth, asymptotically flat initial data set which is "near flat, Minkowski 
data" leads to a unique, smooth and geodesically complete solution of Ein- 



This very simple point was apparently unknown to Einstein in 1917, although 
soon after the publication of his cosmological paper, E. Freundlich and F. Klein 
pointed out to him that an elliptical topology (arising from the identification of 
antipodal points) could have been chosen instead of the spherical one considered 
by Einstein. Although topological questions have been followed with a great in- 
terest in recent decades, the chapter by Geroch and Horowitz in "An Einstein 
Centenary Survey" remains the classic; for more recent texts, see for example 
[B2| and references therein. 
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stein's vacuum equations with vanishing cosmological constant. This demon- 
strates the stability of the Minkowski space with respect to nonhnear (vac- 
uum) perturbations, and the existence of singularity- free, asymptoticaUy flat 
radiative vacuum spacetimes. Christodoulou and Klainerman, however, are 
able to show only a somewhat weaker decay of the field at null infinity than is 
expected from the usual assumption of a sufficient smoothness at null infinity 
in the framework of Penrose (see e.g. |^ for a brief account). 

Curiously enough, in the case of the vacuum Einstein equations with a 
nonvanishing cosmological constant, a more complete picture has been known 
for some time. By using his regular conformal field equations, Friedrich 
demonstrated that initial data sufficiently close to de Sitter data develop into 
solutions of Einstein's equations with a positive cosmological constant, which 
are "asymptotically simple" (with a smooth conformal infinity), as required 
in Penrose's framework. More recently, Friedrich |42 has shown the existence 
of asymptotically simple solutions to the Einstein vacuum equations with a 
negative cosmological constant. For the latest review of Friedrich's thorough 
work on asymptotics, see [Q. 

Summarizing, thanks to these profound mathematical achievements we 
know that the Minkowski, de Sitter, and anti de Sitter spacetimes are the 
solutions of Einstein's field equations which are stable with respect to general, 
nonlinear (though "weak" in a functional sense) vacuum perturbations. A 
result of this type is not known for any other solution of Einstein's equations. 



1.4 On the interpretation and characterization of metrics 

Suppose that a metric satisfying Einstein's field equations is known in some 
region of spacetime and in a given coordinate (reference) system x^. A funda- 
mental question, frequently "forgotten" to be addressed in modern theories 
which extend upon general relativity, is whether the metric tensor gapix^) 
is a measurable guantity. Classical general relativity offers (at least) three 
ways of giving a positive answer, depending on what objects are considered 
as "primitive tools" to perform the measurements. The first, elaborated and 
emphasized primarily by M0ller |Q| , employs standard rigid rods in the mea- 
surements. However, a "rigid rod" is not really a simple primitive concept. 
The second procedure, due to Synge accepts as the basic concepts a 
"particle" and a "standard clock" . If and + dx^ are two nearby events 
contained in the worldline of a clock, then the separation (the spacetime in- 
terval) between the events is equal to the interval measured by the clock. The 
main drawback of this approach appears to lie in the fact that it does not 
explain why the same functions gai3{x^) describe the behavior of the clock as 
well as paths of free particles, as explained in more detail by Ehlers, Pirani 
and Schild p6[ , in the motivation for their own axiomatic but constructive 
procedure for setting up the spacetime geometry. Their method, inspired by 
the work of Weyl and others, uses neither rods nor clocks, but instead, fight 
rays and freely falling test particles, which are considered as basic tools for 
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measuring the metric and determining the spacetime geometry. (For a simple 
description of how this can be performed, see exercise 13.7 in [T^ ; for some 
new devefopments which build upon, among others, the Ehlers-Pirani-Sachs 
approach, see ]4^.) After indicating that the metric tensor is a measurable 
quantity let us briefly turn to the role of spacetime coordinates. 

In special relativity there are infinitely many global inertial coordinate 
systems labelling events in the Minkowski manifold M*; they are related by el- 
ements of the Poincare group. The inertial coordinates labels X", X^, X^, X^ 
of a given event do not thus have intrinsic meaning. However, the spacetime 
interval between two events, determined by the Minkowski metric r]^,^, rep- 
resents an intrinsic property of spacetime. Since the Minkowski metric is 
so simple, the differences between inertial coordinates can have a metrical 
meaning (recall Einstein's reply to Abraham mentioned in Section 1.2). In 
principle, however, both in special and general relativity, it is the metric, 
the line element, which exhibits intrinsically the geometry, and gives all rel- 
evant information. As Misner puts it, if you write down for someone the 
Schwarzschild metric in the "canonical" form (equation in Section 2.2) and 
receive the reaction "that [it] tells me the g^^, gravitational potentials, now 
tell me in which (i, r, 6, ip) coordinate system they have these values?" , then 
there are two valid responses: (a) indicate that it is an indelicate and unneces- 
sary question, or (b) ignore it. Clearly, the Schwarzschild metric describes the 
geometrical properties of the coordinates used in (|^). For example, it implies 
that worldlines with fixed r, 0, ip are timelike at r > 2Af , orthogonal to the 
lines with t = constant. It determines local null cones (given by ds^ = 0), i.e. 
the causal structure of the spacetime. In addition, in Schwarzschild coordi- 
nates the metric (H) indicates how to measure the radial coordinate of a given 
event, because the proper area of the sphere going through the event is given 
just by the Euclidean expression Anr'^ {r is thus often called "the curvature 
coordinate"). On each sphere the angular coordinates 9,(p have the same 
meaning as on a sphere in Euclidean space. The Schwarzschild coordinate 
time t, geometrically preferred by the timelike (for r > 2M) Killing vector, 
which is just equal to d/dt, can be measured by radar signals sent out from 
spatial infinity (r ^ 2M) where t is the proper time (see e.g. fl^). The coor- 
dinates used in ^) are in fact "more unique" than the inertial coordinates in 
Minkowski spacetime, because the only possible continuous transformations 
preserving the form (||) are rigid rotations of a sphere, and t ^ t + constant. 
Such a simple interpretation of coordinates is exceptional. However, the sim- 
ple case of the Schwarzschild metric clearly demonstrates that all intrinsic 
information is contained in the line element. 

It is interesting, and for some purposes useful, to consider not just one 
Schwarzschild metric with a given mass M but the family of such metrics 
for all possible M. In order to cover also the future event horizon let us de- 
scribe the metrics by using Eddington-Finkelstein ingoing coordinates as in 
equation (|^), Section 2.3. This equation can be interpreted as a family of 
metrics with various values of M given on a fixed background manifold M. i , 
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with V ^ ]R,r ^ (0, oo), and 9 ^ [0,Tr], Lp E [0, 27r). Alternatively, however, we 
may use, for example, the Kruskal null coordinates U, V in which the metric 
is given by equation (||), Section 2.3, with U = V — U,V = V + U. We may 
then consider metrics on a background manifold Ai2 given hy U G 1R,V G 
(0,oo), 9 E [0,7r], and ip E [0, 27r), which corresponds to A^i. However, these 
two background manifolds are not the same: the transformation between the 
Eddington-Finkelstein coordinates and the Kruskal coordinates is not a map 
from Ml to 7W2 because it depends on the value of mass M. Therefore, 
the "background manifold" used frequently in general relativity, for example 
in problems of conservation of energy, or in quantum gravity, is not defined 
in a natural, unique manner. The above simple pedagogical observation has 
recently been made in connection with gauge fixing in quantum gravity by 
Hajicek in order to explain the old insight by Bergmann and Komar, that 
the gauge group of general relativity is much larger than the diffeomorphism 
group of one manifold. To identify points when working with backgrounds, 
one usually fixes coordinates in all solution manifolds by some gauge condi- 
tion, and identifies those points of all these manifolds which have the same 
value of the coordinates. 

Returning back to a single solution {Ai,gap), described by a manifold A4 
and a metric gap in some coordinates, a notorious (local) "equivalence prob- 
lem" often arises. A given (not necessarily global) solution has the variety of 
representations which equals the variety of choices of a 4-dimensional coor- 
dinate system. Transitions from one choice to another are isomorphic with 
the group of 4-dimensional diffeomorphisms which expresses the general co- 
variance of the theory.^ Given another set of functions g'^pix''^) which satisfy 
Einstein's equations, how do we learn that they are not just transformed com- 
ponents of the metric gaf}{x'^)l In 1869 E. B. Christoffel raised a more general 
question: under which conditions is it possible to transform a quadratic form 
gaf3{x'^)dx"dx^ in n-dimensions into another such form g'^p{x''*)dx'" dx'^ by 
means of smooth transformation x'^{x"^)7 As Ehlers emphasized in his paper 

As pointed out by Kretschmann soon after the birth of general relativity, one can 
always make a theory generally covariant by taking more variables and insert- 
ing them as new dynamical variables into the (enlarged) theory. Thus, standard 
Yang-Mills theory is covariant with respect to the transformations of Yang-Mills 
potentials, corresponding to a particular group, say SU{2). However, the theory 
is usually formulated on a fixed background spacetime with a given metric. The 
evolution of a dynamical Yang-Mills solution is thus "painted" on a given space- 
time. When the metric - the gravitational field - is incorporated as a dynamical 
variable in the Einstein- Yang-Mills theory, the whole spacetime metric and Yang- 
Mills field are "built-up" from given data (cf. the article by Friedrich and Rendall 
in this volume). The resulting theory is covariant with respect to a much larger 
group. The dual role of the metric, determined only up to 4-dimensional diffeo- 
morphisms, makes the character of the solutions of Einstein's equations unique 
among solutions of other field theories, which do not consider spacetime as being 
dynamical. 
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[ pO| on the meaning of ChristofFel's equivalence problem in modern field the- 
ories, Christoffel's results apply to metrics of arbitrary signature, and can be 
thus used directly in general relativity. Without going into details let us say 
that today the solution to the equivalence problem as presented by Cartan is 
most commonly used. For both metrics gap and one has to find a frame 
(four 1-forms) in which the frame metric is constant, and find the frame com- 
ponents of the Ricmann tensor and its covariant derivatives up to - possibly 
- the 10th order. The two metrics gap and g'^^ are then equivalent if and 
only if there exist coordinate and Lorentz transformations under which one 
whole set of frame components goes into the other. In a practical algorithm 
given by Karlhede (sil , recently summarized and used in |^2| , the number of 
derivations required is reduced. 

A natural first idea of how to solve the equivalence problem is to employ 
the scalar invariants from the Riemann tensor and its covariant derivatives. 
This, however, does not work. For example, in all Petrov type A'^ and /// 
nonexpanding and nontwisting solutions all these invariants vanish as shown 
recently (see Section 8.2), as they do in Minkowski spacetime. 

However, even without regarding invariants, at present much can be learnt 
about an exact solution (at least locally) in geometrical terms, without ref- 
erence to special coordinates. This is thanks to the progress started in the 
late 1950s, in which the group of Pascual Jordan in Hamburg has played the 
leading role, with Jiirgen Ehlers as one of its most active members. Ehlers' 
dissertation]^ Q from 1957 is devoted to the characterization of exact solu- 
tions. 

The problem of exact solutions also forms the content of his contribution 
to the Royaumont GR-conference , as well as his plenary talk in the Lon- 
don GR-conference Q . A detailed description of the results of the Hamburg 
group on invariant geometrical characterization of exact solutions by using 
and developing the Petrov classification of Weyl's tensors, groups of isome- 
trics, and conformal transformations are contained in the first paper [^7| in 
the (today "golden oldies" ) series of articles published in the "Abhandlungen 
der Akademie der Wissenschaften in Mainz" . An English version, in a some- 
what shorter form, was published by Ehlers and Kundt in the "classic" 
1962 book "Gravitation: An Introduction to Current Research" compiled by 
L. Witten. (We shall meet these references in the following sections.) In the 
second paper of the "Abhandlungen" |5^, among others, algebraically spe- 
cial vacuum solutions are studied, using the formalism of the 2-component 
spinors, and in particular, geometrical properties of the congruences of null 
rays are analyzed in terms of their expansion, twist, and shear. 



The English translation of the title of the dissertation reads: "The construction 
and characterization of the solutions of Einstein's gravitational field equations" . 
In |^3| the original German title is quoted, as in our citation |^|, but "of the 
solutions" is erroneously omitted. This error then reemerges in the references in 
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These tools became essential for the discovery by Roy Kerr in 1963 of the 
solution which, when compared with all other solutions of Einstein's equa- 
tions found from the beginning of the renaissance of general relativity in the 
late 1950s until today, has played the most important role. As Chandrasekhar 
eloquently expresses his wonder about the remarkable fact that all sta- 
tionary and isolated black holes are exactly described by the Kerr solution: 
"This is the only instance we have of an exact description of a macroscopic 
object. Macroscopic objects, as we see them all around us, are governed by 
a variety of forces, derived from a variety of approximations to a variety of 
physical theories. In contrast, the only elements in the construction of black 
holes are our basic concepts of space and time ..." The Kerr solution can also 
serve as one of finest examples in general relativity of "the incredible fact 
that a discovery motivated by a search after the beautiful in mathematics 
should find its exact replica in Nature..." |]60| . 

The technology developed in the classical works |5^,0 , and in a number 
of subsequent contributions, is mostly concerned with the local geometrical 
characterization of exact spacetime solutions. A well-known feature of the 
solutions of Einstein's equations, not shared by solutions in other physical 
theories, is that it is often very complicated to analyze their global proper- 
ties, such as their extensions, completeness, or topology. If analyzed globally, 
almost any solution can tell us something about the basic issues in general 
relativity, like the nature of singularities, or cosmic censorship. 

1.5 The choice of solutions 

Since most solutions, when properly analyzed, can be of potential interest, 
we are confronted with a richness of material which puts us in danger of men- 
tioning many of them, but remaining on a general level, and just enumerating 
rather than enlightening. In fact, because of lack of space (and of our under- 
standing) we shall have to adopt this attitude in many places. However, we 
have selected some solutions, hopefully the fittest ones, and when discussing 
their role, we have chosen particular topics to be analyzed in some detail, 
and left other issues to brief remarks and references. 

Firstly, however, let us ask what do we understand by the term "exact 
solution". In the much used "exact-solution-book" the authors "do not 
intend to provide a definition" , or, rather, they have decided that what they 
"chose to include was, by definition, an exact solution" . A mathematical 
relativist-purist would perhaps consider solutions, the existence of which has 
been demonstrated in the works of Friedrich or Christodoulou and Klainer- 
man, mentioned at the end of Section 1.3, as "good" as the Schwarzschild 
metric. Most recently, Penrose presented a strong conjecture which may 
lead to a general vacuum solution described in the complicated (complex) for- 
malism of his twistor theory. Although in this article we do not mean by exact 
solutions those just mentioned, we also do not consider as exact solutions only 
those explicit solutions which can be written in terms of elementary functions 
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on half of a page. We prefer, recalling Feynman, simple "special cases" , but we 
also discuss, for example, the late-time behaviour of the Robinson- Trautman 
solutions for which rigorously convergent series expansions can be obtained, 
which provide sufhciently rich "special information" . 

Concerning the selection of the solutions, the builder of general relativity 
and the gravitational field equations (III) himself indicates which solutions 
should be preferred p3[ |: "The theory avoids all internal discrepancies which 
we have charged against the basis of classical mechanics... But, it is similar 
to a building, one wing of which is made of fine marble (left part of the 
equation) , but the other wing of which is built of low grade wood (right side 
of equation). The phenomenological representation of matter is, in fact, only 
a crude substitute for a representation which would correspond to all known 
properties of matter. There is no difficulty in connecting Maxwell's theory... 
so long as one restricts himself to space, free of ponderable matter and free 
of electric density..." 

Of course, Einstein was not aware when he was writing this of Yang-Mills- 
Higgs fields, or of the dilaton field, etc. However, remaining on the level of 
field theories with a clear classical meaning, his view has its strength and 
motivates us to prefer (electro)vacuum solutions. A physical interpretation 
of the vacuum solutions of Einstein's equations have been reviewed in papers 
by Bonnor [Q, and Bonnor, Griffiths and MacCallum |^ five years ago. 
Our article, in particular in emphasizing and describing the role of solutions 
in giving rise to various concepts, conjectures, and methods of solving prob- 
lems in general relativity, and in the astrophysical impacts of the solutions, 
is oriented quite differently, and gives more detail. However, up to some ex- 
ceptions, like, for example, metrics for an infinite line-mass or plane, which 
are discussed in , and new solutions which have been discovered after the 
reviews |6^|65| appeared as, for example, the solution describing a rigidly 
rotating thin disk of dust, our choice of solutions is similar to that of ||6^ , |65| . 

In selecting particular topics for a more detailed discussion we will be 
led primarily by following overlapping aspects: (i) the "commonly acknowl- 
edged" significance of a solution - we will concentrate in particular on the 
Schwarzschild, the Kerr, the Taub-NUT, and plane wave solutions, and (ii) 
the solutions and their properties that I (and my colleagues) have been di- 
rectly interested in, such as the Reissner-Nordstrom metric, vacuum solu- 
tions outside rotating disks, or radiative solutions such as cylindrical waves, 
Robinson- Trautman solutions, and the boost-rotation symmetric solutions. 
Some of these have also been connected with the interests of Jiirgen Ehlers, 
and we shall indicate whenever we are aware of this fact. 

Vacuum cosmological solutions are discussed in less detail than they de- 
serve. A possible excuse - from the point of view of being a relativist, a rather 
unfair one - could be that a special recent issue of Reviews of Modern Physics 
(Volume 71, 1999), marking the Centennial of the American Physical Soci- 
ety, contains discussion of the Schwarzschild, the Reissner-Nordstrom and 
other black hole solutions, and even remarks on the work of Bondi et al. [pa 
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on radiative solutions, but among the cosmological solutions only the stan- 
dard models are mentioned. A real reason is the author's lack of space, time, 
and energy. In the concluding remarks we will try to list at least the most 
important solutions (not only the Priedmann models!) which have not been 
"selected" and give references to the literature in which more information 
can be found. 

1.6 The outline 

Since the titles of the following sections characterize the contents rather 
specifically, we restrict ourselves to only a few explanatory remarks. In our 
discussion of the Schwarschild metric, after mentioning its role in the solar 
system, we indicate how the Schwarzschild solution gave rise to such concepts 
as the event horizon, the trapped surface, and the apparent horizon. We pay 
more attention to the concept of a bifurcate Killing horizon, because this 
is usually not treated in textbooks, and in addition, Jiirgen Ehlers played 
a role in its first description in the literature. Another point which has not 
received much attention is Penrose's nice presentation of evidence against 
Lorentz-covariant field theoretical approaches to gravity, based on analysis 
of the c;ausal structure of the Schwarzschild spacetime. Among various astro- 
physical implications of the Schwarzschild solution we especially note recent 
suggestions which indicate that we may have evidence of the existence of 
event horizons, and of a black hole in the centre of our Galaxy. 

The main focus in our treatment of the Reissner-Nordstrom metric is 
directed to the instability of the Cauchy horizon and its relation to the cosmic 
censorship conjecture. We also briefly discuss extreme black holes and their 
role in string theory. 

About the same amount of space as that given to the Schwarzschild so- 
lution is devoted to the Kerr metric. After explaining a few new concepts 
the metric inspired, such as locally nonrotating frames and ergoregions, we 
mention a number of physical processes which can take place in the Kerr back- 
ground, including the Penrose energy extraction process, and the Blanclford- 
Znajek mechanism. In the section on the astrophysical evidence for a Kerr 
metric, the main attention is paid to the broad iron line, the character of 
which, as most reec^nt observations indicate, is best explained by assuming 
that it originates very close to a maximally rotating black hole. The dis- 
cussion of recent results on black hole uniqueness and on multi-black hole 
solutions concludes our exposition of spacetimes representing black holes. In 
the section on axisymmetric fields and relativistic disks a brief survey of var- 
ious static solutions is first given, then we concentrate on relativistic disks 
as sources of the Kerr metric and other stationary fields; in particular, we 
summarize briefly the recent work on uniformly rotating disks. 

An intriguing case of Taub-NUT space is introduced by a new constructive 
derivation of the solution. Various pathological features of this space are then 
briefly listed. 
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Going over to radiative spacetimes, we analyze in some detail plane waves 
- also in the light of the thorough study by Ehlers and Kundt j53|. Some 
new developments are then noted, in particular, impulsive waves generated 
by boosting various "particles" , their symmetries, and recent use of the 
Colombeau algebra of generalized functions in the analyses of impulsive 
waves. A fairly detailed discussion is devoted to various effects connected 
with colliding plane waves. 

In our treatment of cylindrical waves we concentrate in particular on two 
issues: on the proof that these waves provide explicitly given spacetimes, 
which admit a smooth global null infinity, even for strong initial data within 
a (2 + l)-dimensional framework; and on the role that cylindrical waves have 
played in the first construction of a midisuperspace model in quantum grav- 
ity. Various other developments concerning cylindrical waves are then sum- 
marized only telegraphically. 

A short section on Robinson- Trautman solutions points out how these 
solutions with a nonvanishing cosmological constant can be used to give an 
exact demonstration of the cosmic no-hair conjecture under the presence of 
gravitational radiation, and also of the existence of an event horizon which 
is smooth but not analytic. 

As the last class of radiative spacetimes we analyze the boost-rotation 
symmetric solutions representing uniformly accelerated objects. They play a 
unique role among radiative spacetimes since they are asymptotically flat, 
in the sense that they admit global smooth sections of null infinity. And as 
the only known radiative solutions describing finite sources they can provide 
expressions for the Bondi mass, the news function, or the radiation patterns in 
explicit forms. They have also been used as test-beds in numerical relativity, 
and as the model spacetimes describing the production of black hole pairs in 
strong fields. 

Vacuum cosmological solutions such as the vacuum Bianchi models and 
Gowdy solutions are mentioned, and their significance in the development of 
general relativity is indicated in the last section. Special attention is paid to 
their role in understanding the behaviour of a general model near an initial 
singularity. 

In the concluding remarks, several important, in particular noK- vacuum 
solutions, which have not been included in the main body of the paper, are 
at least listed, together with some relevant references. A few remarks on the 
possible future role of exact solutions ends the article. 

Although we give over 360 references in the bibliography, we do not at 
all pretend to give all relevant citations. When discussing more basic facts 
and concepts, we quote primarily textbooks and monographs. Only when 
mentioning more recent developments do we refer to journals. The complete 
titles of all listed references will hopefully offer the reader a more complete 
idea of the role the explicit solutions have played on the relativistic stage and 
in the astrophysical sky. 
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2 The Schwarzschild solution 

In his thorough "Survey of General Relativity Theory" Jiirgen Ehlers 
begins with an empirical motivation of the theory, goes in depth and detail 
through his favourite topics such as the axiomatic approach, kinetic theory, 
geometrical optics, approximation methods, and only in the last section turns 
to spherically symmetric spacetimes. As T. S. Eliot says, "to make an end is 
to make a beginning - the end is where we start from" , and so here we start 
with a few remarks on spherical symmetry. 



2.1 Spherically symmetric spacetimes 

In the early days of general relativity spherical symmetry was introduced in 
an intuitive manner. It is because of the existence of exact solutions which are 
singular at their centres (such as the Schwarzschild or the Reissner-Nordstrom 
solutions), and a realization that spherically symmetric, topologically non- 
trivial smooth spacetimes without any centre may exist psf , that today the 
group-theoretical definition of spherical symmetry is preferred (for a detailed 
analysis, see e.g. p9| , p6|,p7[ ) . 

Following Ehlers ipTfl, we define a spacetime {M,ga[3) to be spherically 
symmetric if the rotation group SO3 acts on {^A,ga|3) as an isometry group 
with simply connected, complete, spacelike, 2-dimensional orbits. One can 
then prove the theorem |^,^ that a spherically symmetric spacetime is the 
direct product Ai ~ S'^ x N, where S'^ is the 2-sphere manifold with the 
standard metric gs on the unit sphere; and iV is a 2-dimensional manifold 
with a Lorentzian (indefinite) metric g^, and with a scalar r such that the 
complete spacetime metric gap is "conformally decomposable" , i.e. r~^gap 
is the direct sum of the 2-dimensional parts gM and gs- Leaving further 



technicalities aside (see e.g. [g6|j67||6g[) we write down the final spherically 



symmetric line element in the form 

ds^ = -e^'I'dt^ + e^^dr^ + r^{d9^ + sin^ 9 dip^), (1) 



where (following |67|) we permit 0(r, t) and A(r, t) to have an imaginary part 
in/ 2 so that the signs of dt^ and dr^ in (|l|), and thus the role of r and t 
as space- and time- coordinates may interchange (a lesson learned from the 
vacuum Schwarzschild solutions - see below). The "curvature coordinate" 
r is defined invariantly by the area, Anr'^, of the 2-spheres r = constant, 
t = constant. There is no a priori relation between r and the proper distance 
from the centre (if there is one) to the spherical surface. 



2.2 The Schwarzschild metric and its role in the solar system 

Starting from the line element (|l|) and imposing Einstein's vacuum field equa- 
tions, but allowing spacetime to be in general dynamical, we are led uniquely 
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(cf. Birkhoff's theorem discussed e.g. in ]18| , |26[ ) to the Schwarzschild metric 

ds^ = -(^l- dt^ + (^1 - + {de^ + sin^ 9 dip^) , (2) 



where M — constant has to be interpreted as a mass, as test particle orbits 
show. The resuhing spacetime is static at r > 2M (no spherically symmetric 
gravitational waves exist), and asymptotically flat at r — + cx3. 

Undoubtedly, the Schwarzschild solution, describing the exterior gravi- 
tational field of an arbitrary - static, oscillating, collapsing or expanding - 
spherically symmetric body of (Schwarzschild) mass M, is among the most 
influential solutions of the gravitational field equations, if not of any type 
of field equations invented in the 20th century. It is the first exact solution 
of Einstein's equations obtained - by K. Schwarzschild in December 1915, 
still before Einstein's theory reached its definitive form and, independently, 
in May 1916, by J. Droste, a Dutch student of H. A. Lorentz (see for 
comprehensive survey). 

However, in its exact form (involving regions near r « 2M) the metric 
(H) has not yet been experimentally tested (a more optimistic recent sugges- 
tion will be mentioned in Section 2.6). When in 1915 Einstein explained the 
perihelion advance of Mercury, he found and used only an approximate (to 
second order in the gravitational potential) spherically symmetric solution. In 
order to find the value of the deflection of light passing close to the surface of 
the Sun, in his famous 1911 Prague paper, Einstein used just the equivalence 
principle within his "Prague gravity theory" , based on the variable velocity 
of light. Then, in 1915, he obtained this value to be twice as big in general 
relativity, when, in addition to the equivalence principle, the curvature of 
space (determined from (^ to first order in M/r) was taken into account. 

Despite the fact that for the purpose of solar-system observations the 
Schwarzschild metric in the form (||) is, quoting [Q, "too accurate", it has 
played an important role in experimental relativity. Eddington, Robertson 
and others introduced the method of expanding the Schwarzschild metric 
at the order beyond Newtonian theory, and then multiplying each post- 
Newtonian term by a dimensionless parameter which should be determined by 
experiment. These methods inspired the much more powerful PPN ("Para- 
metrized post-Newtonian" ) formalism which was developed at the end of the 
1960s and the beginning of the 1970s for testing general relativity and alter- 
native theories of gravity. It has been very effectively used to compare general 
relativity with observations (see e.g. [p"8|j7l| , [7^ and references therein). In or- 
der to gain at least some concrete idea, let us just write down the simplest 
generalization of (0), namely the metric 



ds' 



2M 

+ 2 (/3 - 7) — 



r 

+r^ {dd^ + sin"^ 9 d<p^) , (3) 



dt^-^(l + 2j^) dr^ 
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which is obtained by expanding the metric in A//r up to one order beyond 
the Newtonian approximation, and multiplying each post-Newtonian term 
by dimensionless parameters which distinguish the post-Newtonian hmits of 
different metric theories of gravity, and should be determined experimentally. 
(In general, one needs not just two but ten PPN parameters [|l8 
Einstein's theory: P = j — 1. Calculating from metric the advance of the 
pericentre of a test particle orbiting a central mass M on an ellipse with semi- 
major axis a and eccentricity e, one finds A(f> = ^{2 + 2j — (3)6TTM/[a{l — e'^)], 
whereas the total deflection angle of electromagnetic waves passing close to 
the surface of the body is Aip = 2(1 -I- 7)M/ro, where rg is the radius of 
closest approach of photons to the central body. 

Measurements of the deflection of radio waves and microwaves by the Sun 
(recently also of radio waves by Jupiter) at present restrict 7 to ^(1 -|- 7) = 
l.OOOlibO.OOl ||7l| , [72| . Planetary radar rangings, mainly to Mercury, give from 
the perihelion shift measurements the result {2j + 2 — = 1.00 ±0.002, so 
that /? = 1.000 ± 0.003, whereas the measurements of periastron advance for 
the binary pulsar systems such as PSR 1913-1-16 implied agreement with Ein- 
stein's theory to better than about 1% (see e.g. fzi] , !?^ for reviews). There are 
other solar-system experiments verifying the leading orders of the Schwarz- 
schild solution to a high accuracy, such as gravitational redshift, signal retar- 
dation, or lunar geodesic precession. A number of advanced space missions 
have been proposed which could lead to significant improvements in values of 
the PPN parameters, and even to the measurements of post-post-Newtonian 
effects [||. 

Hence, though in an approximate form, the Schwarzschild solution has 
had a great impact on experimental relativity. In addition, the observational 
effects of gravity on light propagation in the solar system, and also today 
routine observations of gravitational lenses in cosmological contexts [^3| , have 
significantly increased our confidence in taking seriously similar predictions 
of general relativity in more extreme conditions. 

2.3 Schwarzschild metric outside a collapsing star 

I recall how Roger Penrose, at the beginning of his lecture at the 1974 Erice 
Summer School on gravitational collapse, placed two figures side by side. The 
first illustrated schematically the bending of light rays by the Sun (surpris- 
ingly, Penrose did not write "Prague 1911" below the figure). I do not remem- 
ber exactly his second figure but it was similar to Fig. 1 below: the spacetime 
diagram showing spherical gravitational collapse through the Schwarzschild 
radius into a spherical black hole. 

It is in all modern books on general relativity that the Schwarzschild 
radius at Rg — 2M is the place where Schwarzschild coordinates t, r are 
unsuitable, and that metric (^) has a coordinate singularity but not a physi- 
cal one. One has to introduce other coordinates to extend the Schwarzschild 
metric through Rg. In order to describe all spacetime outside a collapsing 
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Fig. 1. The gravitational collapse of a spherical star (the interior of the star is 
shaded). The light cones of the three events, O, P, Q, at the centre of the star, and 
of the three events outside the star are illustrated. The event horizon, the trapped 
surfaces, and the singularity formed during the collapse are also shown. Although 
the singularity appears to lie in a "time direction" , from the character of the light 
cone outside the star but inside the event horizon it is seen that it has a spacelike 
character. 



spherical body it is advantageous to use ingoing Eddington-Finkelstein coor- 
dinates {v, r, 0, If) where v = t + r + 2M log(r/2M - 1). Metric (|) takes the 
form 

ds^ = - —j dv^ + 2dvdr + {d9^ + sin^ 6 dip"^) , (4) 

{v, 9, (p) = constant are ingoing radial null geodesies. Fig. 1, plotted in these 
coordinates, demonstrates well several basic concepts and facts which were 
introduced and learned after the end of 1950s when a more complete under- 
standing of the Schwarzschild solution was gradually achieved. The metric 
holds only outside the star, there will be another metric in its interior, for 
example the Oppenheimer-Snyder collapsing dust solution (i.e. a portion of a 
collapsing Friedmann universe), but the precise form of the interior solution 
is not important at the moment. Consider a series of flashes of light emitted 
from the centre of the star at events 0,P, Q (see Fig. 1) and assume that 
the stellar material is transparent. As the Sun has a focusing effect on the 
light rays, so does matter during collapse. As the matter density becomes 
higher and higher, the focusing effect increases. At event P a special wave- 
front will start to propagate, the rays of which will emerge from the surface 
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of the star with zero divergence, i.e. the null vector /c" = dx" /dw, w being 
an affine parameter, tangent to null geodesies, satisfies k"^ — 0. The wave- 
front then "stays" at the hypersurface r = 2M in metric (^) , and the area of 
its 2-dimensional cross-section remains constant. The null hypersurface rep- 
resenting the history of this critical wavefront is the (future) event horizon. 
Note that the light cones turn more and more inwards as the event horizon is 
approached. They become tangential to the horizon in such a way that radial 
outgoing photons stay at r = 2M whereas ingoing photons fall inwards, and 
will eventually reach the curvature singularity at r — 0. As Fig. 1 indicates, 
wavefronts emitted still later than the critical one, as for example that emit- 
ted from event Q, will be focused so strongly that their rays will start to 
converge, and will form (closed) trapped surfaces. The light cones at trapped 
surfaces are so turned inwards that both ingoing and outgoing radial rays 
converge, and their area decreases. 

Consider a family of spacelike hypersurfaces S{t) foliating spacetime (r 
is a time coordinate, e.g. v — r). The boundary of the region of S(t) which 
contains trapped surfaces lying in S{t) is called the apparent horizon in S{t). 

In general, the apparent horizon is different from the intersection of the 
event horizon with S{t), as a nice simple example (based again on an exact 
solution) due to Hawking Q shows. Assume that after the spherical collapse 
of a star a spherical thin shell of mass m surrounding the star collapses 



r=2 (M+m) 




Fig. 2. The "teleological" behaviour of the event horizon during the gravitational 
collapse of a star, followed by the collapse of a shell. The event horizon moves out- 
wards because it will be crossed b^ the shell. The apparent horizon moves outwards 
discontinuously (adapted from 
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and eventually crashes at the singularity at r = (Fig. 2). In the vacuum 
region inside the shell there is the Schwarzschild metric (Q) with mass M, 
and outside the shell with mass M + m. Hence the apparent horizon on S{ti) 
will be at r = 2M and will remain there until S{t2) when it discontinuously 
jumps to r = 2{M + m). One can determine the apparent horizon on a 
given hypersurface. In order to find the event horizon one has to know the 
whole spacetime solution. The future event horizon separates events which 
are visible from future infinity, from those which are not, and thus forms the 
boundary of a black hole. 

From the above example of a shell collapsing onto a Schwarzschild black 
hole we can also learn about the "teleological" nature of the horizon: the mo- 
tion of the horizon depends on what will happen to the horizon in the future 
(whether a collapsing shell will cross it or not). This teleological behaviour of 
the horizon has later been discovered in a variety of astrophysically realistic 
situations such as the behaviour of a horizon perturbed by a mass orbiting a 
black hole (see for enlightening discussions of such effects). 

By studying the Schwarzschild solution and spherical collapse it became 
evident that one has to turn to global methods to gain a full understanding of 
general relativity. The intuition acquired from analyzing the Schwarzschild 
metric helped crucially in defining and understanding such concepts as the 
trapped surface, the event horizon, or the apparent horizon in general situ- 
ations without symmetry. Nowadays these concepts are explained in several 
advanced textbooks and monographs (e.g. | 

Following from the example of spherical collapse one is led to ask whether 
generic gravitational collapses lead to spacetime singularities and whether 
these are always surrounded by an event horizon. The Pcnrose-Hawking sin- 
gularity theorems | p^ , p6[ show that singularities do arise under quite generic 
circumstances (the occurrence of a closed trapped surface is most significant 
for the appearance of a singularity). The second question is the essence of 
the cosmic censorship hypothesis. Various exact solutions have played a role 
in attempts to "prove" or "disprove" this "one of the most important issues" 
of classical relativity. We shall meet it in several other places later on, in 
particular in Section 3.1. There a more detailed formulation is given. 
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2.4 The Schwarzschild-Kruskal spacetime 

In the remarks above we considered the Schwarzschild solution outside a 
static (possibly oscillating, or expanding from r > 2M) star, and outside 
a star collapsing into a black hole. It is not excluded that just these situa- 
tions will turn out to be physically relevant. Nevertheless, in connection with 
the Schwarzschild metric it would be heretical not to mention the enormous 
impact which its maximal vacuum analytic extension into the Schwarzschild- 
Kruskal spacetime has had. This is today described in detail in many places 



(see e.g. |lq,M,E6 76 ). We need two sets of the Schwarzschild coordinates to 



cover the complete spacetime, and we obtain two asymptotically flat spaces. 
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i.e. the spacetime with two ("right" and "left") infinities. The metric in 
Kruskal coordinates V, related to the Schwarzschild r, t (in the regions 
with r>2M) by 

U = ±(r/2M - i)i/2gr-/4M ^Qgj^ (i/4M) , 

V = ±{r/2M - l)i/2e'-/4A^ sinh {t/4M) , (5) 

takes the form 

The introduction of the Kruskal coordinates which remove the singularity of 
the Schwarzschild metric (|^) at the horizon r = 2M and cover the complete 
spacetime manifold (every geodesic either hits the singularity or can be con- 
tinued to the infinite values of its affine parameter), was the most influential 
example which showed that one has to distinguish carefully between just a 
coordinate singularity and the real, physical singularity. It also helped us to 
realize that the definition of a singularity itself is a subtle issue in which the 
concept of geodesic completeness plays a significant role (see for a recent 
analysis of spacetime singularities). 

The character of the Schwarzschild-Kruskal spacetime is best seen in the 
Penrose diagram given in Fig. 3, in which the spacetime is compactified by 
a suitable conformal rescaling of the metric. Both right and left infinities are 
represented, and the causal structure is well illustrated because worldlines 
of radial light signals (radial null geodesies) are 45-degree lines in the dia- 
gram. In particular the black hole region // and a "newly emerged" (as a 
consequence of the analytical continuation) white hole region IV (with the 
white-hole singularity at r = 0) are exhibited. For more detailed analyses 
of the Penrose diagram of the Schwarzschild-Kruskal spacetime the reader 
is referred to e.g. |p^ , |l9|j2^ j7^ . Here we wish to turn in some detail to two 
very important concepts in black hole theory which were first understood 
by the analytic extension of the Schwarzschild solution, and which are not 
often treated in standard textbooks. These are the concepts of the bifurcate 
horizon and of the horizon surface gravity. Jiirgen Ehlcrs played a somewhat 
indirect, but important and noble part in their introduction into literature. 

These concepts were the main subject of the last work of Robert Boyer 
who became one of the victims of a mass murder on August 1, 1966, in 
Austin, Texas. Jiirgen Ehlers was authorized by Mrs. Boyer to look through 
the scientific papers of her husband, and together with John Stachel, pre- 
pared posthumously the paper ||7^ from R. Boyer's notes. Ehlers inserted his 
own discussions, generalized the main theorem on bifurcate horizons, but the 
paper jTSj was published with R. Boyer as the only author. 

In the Schwarzschild spacetime there exists the timelike Killing vector, 
d/dt, which when analytically extended into all Schwarzschild-Kruskal man- 
ifold, becomes null at the event horizon r = 2M, and is spacelike in the 



26 Jin Bicak 



.+ r = (future singularity) ^h- 



I 

r = const. > 2M />' 
r = 2M 'y/J- 












r = >» y/y/^ 


















— i — v^J^ 




















I// t = 









i r = (past singularity) ' 



Fig. 3. The Penrose diagram of tlie compactified Scliwarzscliild-Kruskal spacetime. 
Radial null geodesies are 45-degrees lines. Timelike geodesies reach the future (or 
past) timelike infinities i"*" (or null geodesies reaeh the future (or past) null 
infinities J'^ (or J^~) and spaeelike geodesies lead to spatial infinities i" . (Notice 
that at i° the lines t — constant are tangent to each other - this is often not taken 
into account in the literature - see e.g. p6|j3o| .) 



regions // and IV with r < 2M. In Kruskal coordinates it is given by 

k" = (fc^ = U/4M, = V/AM, k" = 0, fc*^ = 0) . (7) 

Hence it vanishes at all points with U = V = 0, 6* e [0,7r],(^ e [0, 27r). 
These points, forming a spacelike 2-sphcrc which we denote B (in Schwarz- 
schild coordinates given by r = 2M, t — constant), are fixed points of the 
1-dimensional group G of isometries generated by k°' (see Fig. 3). At the 
event horizon the corresponding 1-dimensional orbits are null geodesies, with 
k" being a tangent vector. However, since /c" vanishes at _B, these orbits are 
incomplete. 

This (and similar observations for other black hole solutions) motivated 
a general analysis of the bifurcate Killing horizons given in [Q. There it 
is proven for spacetimes admitting a general Killing vector field which 
generates a 1-dimensional group of isometries, that (i) a 1-dimensional orbit 
is a complete geodesic if the gradient of the square vanishes on the orbit, 
(ii) if a geodesic orbit is incomplete, then it is null and (^^),q 0. In addition, 
if ^" = dx°' / dv (v being the group parameter), the affine parameter along 
the geodesic is w — e"^, where k — constant satisfies 

(_^2),a ^ 2kC- (8) 

In the Schwarzschild case, with = fc" = {d/dt)°', and considering the 
part V — U oi the horizon, we get k — 1/4M. The relation w — e'^" is 
just the familiar equation V = e"/^^^, where V = V + U is the Kruskal null 
coordinate and v is the Eddington-Finkelstein ingoing null coordinate used 
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in (Q) . (Notice that V is indeed the afRne parameter along the nuU geodesies 
at the horizon V = U.) The quantity k, first introduced in jT^ , has become 
fundamental in modern black hole theory, and also in its generalizations in 
string theory. It is the well-known surface gravity of the black hole horizon. 

With K 7^ 0, the limit points corresponding to v —oo,w = are fixed 
points of G. (Unless the spacetime is incomplete, there exists a continuation 
of each null geodesic beyond these fixed points to u; < 0.) One can show 
that the fixed points form a spacelike 2-dimensional manifold B, given by 
U = V = in the Schwarzschild case; this "bifurcation surface" is a totally 
geodesic submanifold. By the original definition |^^, a Killing horizon is a G 
invariant null hypersurface N on which =0. (A recent definition |^ , ^ 
specifies a Killing horizon to be any union of such hypersurfaces.) If k 7^ 0, 
at each point of B there is one null direction orthogonal to B which is not 
tangent to N ^ N U B. The null geodesies intersecting B in these directions 
form another null hypersurface, N, which is also a Killing horizon. The union 
U is called a bifurcate Killing horizon (Fig. 4). 



Fig. 4. The bifurcate Killing horizon consisting of two null hypersurfaces N and 
N which intersect in the spacelike 2-dimensional "bifurcation surface" B. 



Bifurcate Killing horizons exist also in flat and other curved spacetimes. 
For example, in the boost-rotation symmetric spacetimes (Section 11), null 
hypersurfaces z = ±t form the bifurcate Killing horizon corresponding to the 
boost Killing vector; B, given by z = < = 0, is then not compact. (As in the 
Schwarzschild-Kruskal spacetime, a bifurcate Killing horizon locally divides 
the spacetime into four wedges.) However, the first motivation for analyzing 
Killing horizons came from the black hole solutions. 

Both Killing horizons and surface gravity play an important role in black 
hole thermodynamics and quantum field theory on curved backgrounds p3|, 
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in particular in their two principal results: the Hawking effect of particle cre- 
ation by black holes; and the Unruh effect showing that a thermal bath of 
particles will be seen also by a uniformly accelerated observer in flat space- 
time when the quantum field is in its vacuum state with respect to inertial 
observers. Recently, new results were obtained which support the view 
that a spacetime representing the final state of a black hole formed by collapse 
has indeed a bifurcate Killing horizon, or the Killing horizon is degenerate 
(« = 0). 

2.5 The Schwarzschild metric as a case against Lorentz-covariant 
approaches 

There are many other issues on which the Schwarzschild solution has made an 
impact. Some of astrophysical applications will be very briefly mentioned later 
on. As the last theoretical point in this section 1 would like to discuss in some 
detail the causal structure of the Schwarzschild spacetime including infinity. 
By analyzing this structure, Penrose js^ presented evidence against various 
Lorentz (Poincare)-covariant field theoretical approaches, which regard the 
physical metric tensor g to be not much different from any other tensor 
in Minkowski spacetime with flat metric 77 (see e.g. [^,|6j). I thought it 
appropriate to mention this point here, since Jiirgen Ehlers, among others, 
certainly does not share a field theoretical viewpoint. 

The normal procedure of calculating the metric g in these approaches is 
from a power series expansion of Lorentz-covariant terms (in quantum the- 
ory this corresponds to an infinite summation of Feynman diagrams). The 
derived field propagation has to follow the true null cones of the curved metric 
g instead of those of r/. However, as Penrose shows, in a satisfactory theory 
the null cones defined by g should not extend outside the null cones de- 
fined by 77, or "the causality defined by g should not violate the background 
77-causality" . Following [Q, let us write this condition as g < rj. Now at 
first sight we may believe that g < 77 is satisfied in the Schwarzschild field 
since its effect is to "slow down" the velocity of light (cf. "signal retarda- 
tion" mentioned in 2.2). However, in the field-theoretical approaches one of 
the main emphasis is in a consistent formulation of scattering theory. This 
requires a good behaviour at infinity. But with the Schwarzschild metric, null 
geodesies with respect to metric g "infinitely deviate" from those with respect 
to 77: for example, the radial outgoing g null geodesies 6,ip — constant, and 
u ^ t — r — 2Mlog(r/2M — f) = constant at r — > 00 go "indefinitely far" 
into the retarded time t — r oi rj, and hence, do not correspond to outgoing 
7y-null geodesies t — r = constant. One can try to use a different flat metric 
associated with the Schwarzschild metric g which does not lead to patho- 
logical behaviour at infinity, but then it turns out that < 77 is violated 
locally. In fact, Penrose Q| proves the theorem, showing that there is an 
essential incompatibility between the causal structures in the Schwarzschild 
and Minkowski spacetimes which appears either asymptotically or locally. 
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This incompatibility is easily understood with the exact Schwarzschild 
solution, but it is generic, since one is concerned only with the behaviour of 
the space at large distances from a positive-mass source, i.e. with the causal 
properties in the neighbourhood of spacelike infinity i^. 

In the present post-Minkowskian approximation methods for the gener- 
ation of gravitational waves by relativistic sources, a suitable (Bondi-type) 
coordinate system is constructed at all orders in the far wave zone, which 
in particular corrects for the logarithmic deviation of the true light cones with 
respect to the coordinate flat light cones (cf. contribution by L. Blanchct in 
this volume). 

2.6 The Schwarzschild metric and astrophysics 

In his introductory chapter "General Relativity as a Tool for Astrophysics" 
for the Seminar in Bad Honnef in 1996 |8^, Jiirgen Ehlers remarks that "The 
interest of black holes for astrophysics is obvious... The challenge here is to 
find observable features that are truly relativistic, related, for example, to 
horizons, ergoregions... Indications exist, but - as far as I am aware - no firm 
evidence." 

There are many excellent recent reviews on the astrophysical evidence for 
black holes (see e.g. |88| , ^^ ). It is true, that the evidence points towards 
the presence of dark massive objects - stellar-mass objects in binaries, and 
supermassive objects in the centres of galaxies - which are associated with 
deep gravitational potential wells where Newtonian gravity cannot be used, 
but it does not offer a clear diagnostic of general relativity. 

Many investigations of test particle orbits in the strong-gravity regions 
(r < lOM) have shown basic differences between the motion in the Schwarz- 
schild metric and the motion in the central field in Newton's theory (e.g. 
[ p^|]7^ , |9l[ ) . For example for 3M < r < 6M unstable circular particle orbits 
exist which are energetically unbound, and thus perturbed particles may es- 
cape to infinity; at r = 3M circular photon orbits occur and there are no 
circular orbits for r < 3M. Particles are trapped by a Schwarzschild black 
hole if they reach the region r < 3M. 

About ten years ago, the study of the behaviour of particles and gyro- 
scopes in the Schwarzschild field revived interest in the "classical" problem of 
the definition of gravitational, centrifugal, and other inertial "forces" acting 
on particles and gyros moving on the Schwarzschild or on a more general 
curved backgrounds, usually axisynimetric and stationary (see e.g. | |9^ , p3t , 
and many references therein). One would like to have a split of a covariantly 
defined quantity (like an acceleration) into non-covariant parts, the physical 
meaning of which would increase our intuition of relativistic effects in astro- 
physical problems. If, for example, we adopt the view that the "gravitational 
force" is velocity-independent, then we find that at the orbits outside the cir- 
cular photon orbit (r > 3M), the centrifugal force is as in classical physics. 
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repulsive, while it becomes attractive inside this orbit, being zero exactly at 
the orbit 

Rclativistic effects will, of course, play a role in many astrophysical situ- 
ations involving spherical accretion, the structure of accretion disks around 
compact stars and black holes, their optical appearance etc. They have be- 
come an important part of the arsenal of astrophysicists, and they have en- 
tered standard literature (see e.g. Though this whole field of science 
lies beyond the scope of this article, I would like to mention three recent issues 
which provide us with hope that we may perhaps soon meet the challenge 
noted in Jiirgen Ehlers' remarks made in Bad Honnef in 1996. 

The first concerns our Galactic centre. Thanks to new observations of stars 
in the near infrared band it was possible to detect the transverse motions of 
stars (for which the radial velocities are also observed) within 0.1 pc in our 
Galactic centre. The stellar velocities up to 2000 km/sec and their dependence 
on the radial distance from the centre are consistent with a black hole of mass 
2.5 X IO^Mq. In the opinion of some leading astrophysicists, our Galactic 
centre now provides "the most convincing case for a supermassive hole, with 
the single exception of NGC 4258" |8|. (In NGC 4258 a disk is observed 
whose inner edge is orbiting at 1080 km/sec, implying a black hole - "or 
something more exotic" - with a mass of 3.6 x IO^Mq.) Perhaps wc shall 
be able to observe rclativistic effects on the proper motions of stars in our 
Galactic centre in the not too distant future. 

The second issue concerns the fundamental question of whether observa- 
tions can bring convincing proof of the existence of black hole event hori- 
zons. Very recently some astrophysicists claimed that new observations, 
in particular of X-ray binaries, imply such evidence. The idea is that thin 
disk accretion cannot explain the spectra of some of X-ray binaries. One has 
to use a different accretion model, a so called advection-dominated accretion 
flow model (ADAF) in which most of the gravitational energy released in 
the infalling gas is carried (advected) with the flow as thermal energy, which 
falls on the central object. (In thin disks most of this energy is radiated out 
from the disk.) If the central compact object (for example a neutron star) 
has a hard surface, the thermal energy stored in the flow is re-radiated after 
the flow hits the surface. However, some of the X-ray binaries show such low 
luminosities that a very large fraction of the energy in the flow must be ad- 
vected through an event horizon into a black hole Although Rees Jssf , 
for example, considers this evidence "gratifyingly consistent with the high- 
mass objects in binaries being black holes", he believes that it "would still 
not convince an intelligent sceptic, who could postulate a different theory 
of strong-field gravity or else that the high-mass compact objects were (for 
instance) self-gravitating clusters of weakly interacting particles...". 

^ Curiously enough, Feynman in his 1962-63 lectures on gravitation |^] writes 
that "inside r — 2M [not 3Af!]... the 'centrifugal force' apparently acts as an 
attraction rather than a repulsion" . 
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For a sceptical optimistic relativist, the most challenging observational 
issue related to black holes probably is to find astrophysical evidence for a 
Kerr metric. We shall come to this point in Section 4.3. 

The last (but certainly not the least) issue lies more in the future, but 
eventually should turn out to be most promising. It is connected with both 
the Numerical Relativity Great Challenge Alliance and the "great challenge" 
of experimental relativity: to calculate reliable gravitational wave-forms and 
to detect them. When gravitational waves from stars captured by a super- 
massive black hole, or from a newly forming supermassive black hole, or, 
most importantly, from coalescing supermassive holes will be detected and 
compared with the predictions of the theory, we should learn significant facts 
about black holes • Are these so general remarks entirely inappropriate 

in the section on the Schwarzschild solution? 

One of the most important roles of the Schwarzschild solution in the 
development of mathematical relativity and especially of relativistic astro- 
physics stems from its simplicity, in particular from its spherical symmetry. 
This has enabled us to develop the mathematically beautiful theory of lin- 
ear perturbations of the Schwarzschild background and employ it in various 
astrophysically realistic situations (see e.g. [[75U7^ , |9l|| , and many references 
therein). Surprisingly enough, this theory does not only give reliable results 
in such problems as the calculation of waves emitted by pulsating neutron 
stars, or waves radiated out from stars falling into a supermassive black hole. 
Very recently we have learned that one can use perturbation theory of a single 
Schwarzschild black hole as a "close approximation" to black hole collisions. 
Towards the end of the collision of two black holes, they will not in fact be 
two black holes, but will merge into a highly distorted single black hole JOSf . 
When compared with the numerical results on a head-on collision it has been 
found that this approximation gives predictions for separations A as large as 
A/M - 7. 

3 The Reissner-Nordstrom solution 

This spherically symmetric solution of the Einstein-Maxwell equations was 
derived independently^ by H. Reissncr in 1916, H. Weyl in 1917, and G. 

^'^ In the literature one finds the solution to be repeatedly connected only with 
the names of Reissner and Nordstrom, except for the "exact-solutions-book" 
[ pH : there in four places the solution is called as everywhere else, but in one 
place (p. 257) it is referred to as the "Reissner- Weyl-solutions" . An enlightening 
discussion on p. 209 in ^] shows that the solution belongs to a more general 
"Weyl's electrovacuum class" of electrostatic solutions discovered by Weyl (in 
1917) which follow from an Ansatz that there is a functional relationship between 
the gravitational and electrostatic potentials. As will be noticed also in the case 
of cylindrical waves in Section 9, if "too many" solutions are given in one paper, 
the name of the author is not likely to survive in the name of an important 
subclass... 
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Nordstrom in 1918. It represents a spacetime with no matter sources except 
for a radial electric field, the energy of which has to be included on the right- 
hand side of the Einstein equations. Since Birkhoff 's theorem, mentioned in 
connection with the Schwarzschild solution in Section 2.2, can be generalized 
to the electrovacuum case, the Reissner-Nordstrom solution is the unique 
spherical electrovacuum solution. Similarly to the Schwarzschild solution, it 
thus describes the exterior gravitational and electromagnetic fields of an ar- 
bitrary - static, oscillating, collapsing or expanding - spherically symmetric, 
charged body of mass M and charge Q. The metric reads 

\ r ) \ r r'^ ) 

+ [dO^ + sin^ 61 , (9) 

the electromagnetic field in these spherical coordinates is described by the 
"classical" expressions for the time component of the electromagnetic poten- 
tial and the (only non-zero) component of the electromagnetic field tensor: 

At=-^, Ftr^^Frt^-^. (10) 

A number of authors have discussed spherically symmetric, static charged 
dust configurations producing a Reissner-Nordstrom metric outside, some of 
them with a hope to construct a "classical model" of a charged elementary 
particle (see |^ for references) . The main influence the metric has exerted on 
the developments of general relativity, and more recently in supersymmetric 
and superstring theories (see Section 3.2), is however in its analytically ex- 
tended electrovacuum form when it represents charged, spherical black holes. 



3.1 Reissner-Nordstrom black holes and the question of cosmic 
censorship 

The analytic extensions have qualitatively different character in three cases, 
depending on the relationship between the mass M and the charge Q. In 
the case > Af^ (corresponding, for example, to the field outside an elec- 
tron), the complete electrovacuum spacetime is covered by the coordinates 
(t,r,9,ip), < r < oo. There is a naked singularity (visible from infinity) at 
r = in which the curvature invariants diverge. If < M^, the metric (|) 
describes a (generic) Reissner-Nordstrom black hole; it becomes singular at 
two radii: 

r = r±^M±{M^-Q^)i. (11) 

Similarly to the Schwarzschild case, these are only coordinate singularities. 
Graves and Brill |9^ discovered, however, that the analytic extension and 
the causal structure of the Reissner-Nordstrom spacetime with > is 
fundamentally different from that of the Schwarzschild spacetime. There are 
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two null hypersurfaces, at r = ?'+ and r ~ r^, which are known as the outer 
(event) horizon and the inner horizon; the Killing vector d / dt is null at the 
horizons, timelike at r > r_|_ and r < r„, but spacelike at r_ < r < r+. The 
character of the extended manifold is best seen in the Penrose diagram in Fig. 
5, in which the spacetime is compactified by a suitable conformal rescaling of 
the metric (see, e.g. [ p"8|j2^p0t ) . As in the compactified Kruskal-Schwarzschild 
diagram in Fig. 3, the causal structure is well illustrated because worldlines 
of radial light signals are 45-degree lines. There are again two infinities il- 
lustrated - the right and left - in regions / and ///. However, the maximally 
extended Reissner- Nordstrom geometry consists of an infinite chain of asymp- 
totic regions connected by "wormholes" between the real singularities (with 
divergent curvature invariants) at r = 0. In Fig. 5, the right and left (past 
null) infinities in regions /' and ///' are still seen - the others are obtained 
by extending the diagram vertically in both directions. 




Fig. 5. The compactified Reissner-Nordstrom spacetime representing a non- 
extreme black hole consists of an indefinite chain of asymptotic regions ( "universes" ) 
connected by "wormholes" between timelike singularities. The worldline of a shell 
collapsing from "universe" I and re-emerging in "universe" I' is indicated. The 
inner horizon at r = r_ is the Cauchy horizon for spacelike hypersurface E. It is 
unstable and will thus very likely prevent such a process occuring. 

An important lesson one has learned is that the character of the singu- 
larity need not be spacelike as it is in the Schwarzschild case, or with the 
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big bang singularities in standard cosmological models. Indeed, the singu- 
larities in the Reissner-Nordstrom geometry are timelike: they do not block 
the way to the future. By solving the geodesic equation one can show that 
there are test particles which start in "our universe" (region 7), cross the 
outer horizon at r = r+ and the inner horizon at r = r_, avoid the singu- 
larity and through a "white hole" (the outer horizon between regions IV' 
and /') emerge into "another universe" /' with its own asymptotically flat 
region. Such a gravitational bounce can occur not only with test parti cles. The 
studies of the gravitational collapse of charged spherical shells ( |lO0|] and ref- 



erences therein) and of charged dust spheres ([ 101 and references therein) 
have shown that a bounce can take place also in fully dynamical cases.|^ The 
part of Fig. 5 which is "left" from the worldline of the surface of the sphere 
or the shell is "covered" by t he interior of the sphere or flat space inside 



the shell. As observed in |10C|, the outcome of the bounce of a shell can be 
different, depending on the value of the shell's total mass, charge and rest 
mass. The shell may crash into the "right" singularity or it may continue to 
expand and emerge in region /'. If the rest mass of the shell is negative the 
collapse may even lead to a naked singularity. 

Now even if the shell collapses into a black hole, and after a bounce, 
emerges in region /', a locally naked singularity is present: the timelike sin- 
gularity at r = 0, to the "right" from the wordline of the shell. An observer 
travelling into the future "between" the shell and the singularity can be sur- 
prised by a signal coming from the singularity (see Fig. 5). Penrose's strong 



cosmic censorship conjecture (see e.g. 102 ) suggests that this should not 
happen. In its physical formulation, as given by Wald (see also for the 
precise formulation), it says that all physically reasonable spacetimes are 
globally hyperbolic, i.e. apart from a possible initial (big bang-type) singu- 
larity, no singularity is visible to any observer. It was just the example of the 
Reissner-Nordstrom solution (and a similar property of the Kerr solution) 
which inspired Penrose to formulate the strong cosmic censorship conjec- 
ture, in addition to its weak version which only requires that from generic 
nonsingular initial data on a Cauchy hypersurface no spacetime singularity 



develops which is visible from infinity. As Penrose | 102 | puts it, "it seems to 
be comparatively unimportant whether the observer himself can escape to 
infinity" . 

It is evident from Fig. 5 that the Reissner-Nordstrom spacetime is not 
globally hyperbolic, i.e. it does not possess a Cauchy hypersurface E, the 



An intuitive explanation [101| of this bounce is that as the sphere (the shell) 
contracts, the volume of the exterior region increases, and hence also the total 
energy in the electric field, which eventually exceeds the energy in the sphere. 
However, the external plus internal energy does not change during collapse (there 
are no waves), and so in the neighborhood of a highly contracted charged object, 
the gravitational field must have a repulsive character corresponding to a negative 
mass-energy. 
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initial data on which (for a test field, say) would determine the development 
of data in the entire future. If data are given on the spacelike hypersurface 
S "connecting" left and right infinities of regions /// and /, the Cauchy 
development will predict what happens only in regions ///and / above 17, and 
in region //, i.e. not beyond the null hypersurfaces (inner horizons) r = r_ 
between region // and regions V in the figure. The inner horizons r = 
represent the Cauchy horizon for a typical initial hypersurface like S. As 
noticed above, what is happening at an event in regions V is in general 
influenced not only by data on S but also by what is happening at the 
(locally) naked singularities (which cannot be predicted since the physics at 
a singularity cannot be controlled). 

Penrose was also the first who predicted that the inner (Cauchy) horizon 



is unstable |27 . If this is true, a null singularity, or possibly even a spacelike 
singularity may arise during a general collapse, so preventing a violation of the 
strong cosmic censorship conjecture. The instability of the Cauchy horizon 
can in fact be expected by using first the following simple geometrical-optics 
argument. 

Introduce the ingoing Eddington-Finkelstcin null coordinate v by 

v^t + r*=t + j f{r)dr, 

/ = l-2M/r + QVr2, (12) 

which brings the metric (^ into the form as of equation in the Schwarz- 
schild case. Consider a freely falling observer who is approaching the inner 
horizon given by r = r-,v — oo. Denoting the observer's constant spe- 
cific energy parameter (see e.g. |l8|) hy E = —U^, where the Killing vector 
^ = d/dv, observer's four-velocity U = d/dr (t - observer's proper time), 
the geodesic equations imply 

r'+f = E\ V = f-'[E - (E' - f)h (13) 

where a dot denotes d/dr. Between the horizons, ^ is spacelike and E can 
be negative. Geodesic equations ( p^ ) imply dr/dv — ^/ for an observer with 
E < 0, approaching r — 7'_ from region //. Expanding / near r — r-, 
f — ,/'(^-)(^ ^ ^-) = — 2K(r — r_), where 

is the surface gravity of the inner horizon (as it follows from definition (||)), 
and integrating, we get / near r_. From the second equation in ( |l3|) we then 
obtain the "asymptotic formula" 

V ~ constant I E I e^'^ = \ , v ^ oo, At 0, (15) 

kAt 

where At is the amount of proper time the observer needs to reach the inner 
(Cauchy) horizon. 
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Imagine now two nearby events Aout and Bout in the outside world /, for 
example the emission of two photons from a given fixed r > r-|_, which are 
connected by the ingoing null geodesies v = constant and v + dv = constant 
with events Ain and Bin on the worldline of the observer approaching r_ . The 
interval of proper time between Ao^t and Bo^t is dTout ^ dv, whereas ( p^ 
implies that the interval of proper time between Ain and Bin , as measured 
by the observer approaching r_, is dr^n ~ e~™dv. Therefore, 

^e-«^ (16) 



so that d& V ^ +00 the events (clocks) in the outside world are measured to 
proceed increasingly fast by the inside observer approaching the inner hori- 
zon. In the limit, when the observer crosses the Cauchy horizon, he sees the 
whole future history (from some event as Aout) of the external universe to 
"proceed in one flash": drin — > with dTout 00. An intuitive explanation 
is given by the fact that the observers in region / need infinite proper time 
to reach v = 00, whereas inside observers only finite proper time. The in- 
falling radiation will thus be unboundedly blue-shifted at the inner horizon 
which in general will lead to a divergence of the energy density there. This 
infinite blueshift at the inner (Cauchy) horizon makes it generally unstable 
to perturbations ("the blue-sheet instability"). 

There exists an extensive literature analyzing the exciting questions of 



the black hole interiors (see for example [Q, the introductory review [103| 
in the proceedings of the recent workshop devoted entirely to these issues, 
and other contributions in the proceedings which give also many further 
references). Some crucial questions are still the subject of much debate. One 
of the following two approaches to the problem is usually chosen: (i) a linear 
perturbation analysis of the behaviour of fields at the Cauchy horizon, (ii) 
the simplified nonlinear, spherically symmetric models of black hole interiors. 

In the first approach one considers the evolution of linear perturbations, 
representing scalar, electromagnetic, or gravitational fields, on the Reissner- 
Nordstrom background. Since there is a nonvanishing background electric 
field, the electromagnetic and gravitational perturbations are cottpferf.pl It is 
a remarkable fact that "wave equations" for certain gauge-invariant combi- 
nations of perturbations can be derived from which all perturbations can 



eventually be constructed ^ , 105| , I06| . In the simplest case of the scalar field 



This leads to various interesting phenomena. For example, the scattering of in- 
cident electromagnetic and gravitational waves by the Reissner-Nordstrom black 
hole allows for the partial conversion of electromagnetic waves into gravitational 
waves and vice versa jo^. When studying stationary electromagnetic fields due 
to sources located outside the Reissner-Nordstrom black hole, one discovers that 
closed magnetic field lines not linking any current source may exist, since gravita- 
tional pertu rbations constitute, via the background Maxwell tensor, an efi'ective 
source [104]. 
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on the Reissner- Nordstrom background, after resolving the field into spher- 



ical harmonics and putting — , the wave equation has the form |107| 

'^,U - ^,r'r' + Fi{r*)^ = (^, (17) 
where the curvature-induced potential barrier is given by 



(18) 



where r is considered to be a function of r* (cf. (p^). In order to determine 
the evolution of the field below the outer horizon in a real gravitational col- 
lapse, one first concentrates on the evolution of the field outside a collapsing 
body (star). A nonspherically symmetric scalar test field (generated by a non- 
spherical distribution of "scalar charge" in the star) serves as a prototype for 
(small) asymmetries in the external gravitational and electromagnetic fields, 
which are generated by asymmetries in matter and charge distributions inside 
the star. Now when a slightly nonspherical star starts to collapse, the per- 
turbations become dynamical and propagate as waves. Their evolution can 
be determined by solving the wave equation ([l^). Because of the potential 
barrier (|l^) the waves get backscattered and produce slowly decaying radia- 
tive tails, as shown in the classical papers by Price ]10g ,109|, and generalized 



to the Reissner-Nordstrom case in |107[|11C |. The tails decay in the vicinity 



of the outer event horizon r+ (i.e. between regions / and II in Fig. 5) as 
^ ~ ^-2(i-i-i) for Z-pole perturbations .0 The decaying tails provide the initial 
data for the "internal problem" - the behaviour of the field near the Cauchy 
horizon. Calculations show (see |l03[ and references therein) that near the 
Cauchy horizon the behaviour of the field remains qualitatively the same: 
^{u^v^oo) ~ y^'^(^+'^) -\- {slowly varying function of u}, where u = 2r* — v 
is constant along outgoing radial null geodesies in region //. However, as a 
consequence of the "exponentially growing blueshift" , given by formula (|l5|) , 
the rate of change of the field diverges as the observer approaches the Cauchy 
horizon: d^/dr = {d^/dx°')U°' ~ {; ~ y-'^^-^e'^" . Therefore, the mea- 
sured energy density in the field would also diverge, causing an instability of 
the Cauchy horizon, which would be expected to create a curvature singular- 



ity. More detailed considerations |103| show that the singularity, at least for 



large \u\, is null and weak (the metric is well-defined, only the Riemann ten- 
sor is singular). Any definitive picture of the Cauchy horizon instability can 
come however only from a fully nonlinear analysis which takes into account 
the backreaction of spacetime geometry to the growing perturbations. 

The second approach to the study of the Cauchy horizon instabilities em- 
ploys a simplified, spherically symmetric model which treats the nonlineari- 



ties exactly [ 1 1 1 1 . The ingoing radiation is modelled by a stream of ingoing 

This result is true for a general charged Reissner-Nordstrom black hole with 
q2 < M^. In the extremal case, = M^, the field decays only as !i/ ~ 
[107|. 
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charged null dust [112 which is infinitely blueshifted at the inner horizon. 
There is, however, also an outgoing stream of charged null dust considered 
to propagate into region // towards the inner horizon. The outgoing flux may 
represent radiation coming from the stellar surface below the outer horizon, 
as well as a portion of the ingoing radiation which is backscattered in region 
J/, and irradiates thus the inner horizon. A detailed analysis based on exact 
spherically symmetric solutions revealed a remarkable effect: an effective in- 
ternal gravitational- mass parameter of the hole unboundedly increases at the 
inner (Cauchy) horizon (though the external mass of the hole remains finite). 
This "mass inflation phenome non" causes the divergence of some curvature 
scalars at the Cauchy horizon | 111 ] . In reality, the classical laws of general 
relativity will break down when the curvature reaches Planckian values. 

It is outside the scope of this review to discuss further the fascinating 
issues of black hole interiors. They involve deep questions of classical rela- 
tivity, of quantum field theory on curved background (as, for example, in 
discussions of electromagnetic pair production and vacuum polarization ef- 
fects inside black holes) , and they lead us eventually to quantum gravity. We 
refer again especially to |7^] and |102l for more information. Let us only add 
three further remarks. We mentioned above the work on the inner structure 
of Reissner- Nordstrom black holes because this is the most explored (though 
not closed) area. However, Kerr black holes (Section 4) possess also inner 
horizons and there are many papers concerned with the instabilities of the 
Kerr Cauchy horizons (see |7£,103| for references). Secondly, at the beginning 
of 1990s, it was shown that the inner horizons of the Reissner-Nordstrom-de 
Sitter and Kerr-de Sitter black holes are classically stable in the case when 
the surface gravity at the inner horizons is smaller than the surface gravity 
at the cosmological horizon (|103| and references therein, in particular, the 
review [113|). Penrose |114| even suggested that "it may well be that cosmic 
censorship requires a zero (or at least a nonpositi ve) cosmological constant" . 
Very recently, however, three experts in the field [ 115 have claimed that out- 
going modes near to the black hole (outer) event horizon lead to instability 
for all values of the parameters of Reissner-Nord strom -de Sitter black holes. 
Let me borrow again a statement from Penrose 114 : "My own feelings are 
left somewhat uncertain by all these considerations" . 

Finally, a new contribution [[ll6| to the old problem of testing the weak 
cosmic censorship by employing a Reissner-Nordstrom black hole indicates 
that one can overcharge a near extreme {Q^ — > Af^) black hole by throwing in 
a charged particle appropriately. However, the backreaction effects remain to 
be explored more thoroughly. The question of cosmic censorship thus remains 
as interesting as ever. 
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3.2 On extreme black holes, d-dimensional black holes, string 
theory and "all that" 

In the previous section we considered generic Reissner- Nordstrom black holes 
with > . They have outer and inner horizons given by (^ij), with 
nonvanishing surface gravities (cf. (Q) for the inner horizon). For Af^ = 
the two horizons coincide at r+ = r_ = M. Defining the ingoing null 
coordinate v as in (13), we obtain the ingoing extension of the Reissner- 
Nordstrom metric (p) in the form 

ds^ = -(\~^^ dr^ + 2dvdr + r^{d0'^ + sin^9dip^). (19) 

This is the metric of extreme Reissner-Nordstrom black holes. Frequently, 
these holes are called "degenerate". At the horizon r — AI, the Killing vector 
field k = d/dv obeys the equation {k^ka),/^ = 0, so that regarding the general 
relation m), the surface gravity k = 0, i.e. the Killing horizon is degenerate. 
Using {k^^i3 — and the Killing equation, we easily deduce that the hori- 
zon null generators with tangent fc" — dx°' /dv satisfy the geodesic equation 
with affine parameter v. The generators have infinite afRne length to the past 
given by w ^ — oo (in contrast to the generators of a bifurcate Killing horizon 
- cf. Section 2.4). This part of the extreme Reissner-Nordstrom spacetime, 
given by r = M^v — > — cxd, is called an "internal infinity". That there is no 
"wormhole" joining two asymptotically flat regions and containing a minimal 
surface 2-sphere like in the non-extreme case can also be seen from the met- 
ric in the original Schwarzschild-type coordinates. Considering an embedding 
diagram t = constant, 6* = 7r/2 in flat Euclidean space one finds that an in- 
finite "tube" , or an asymptotically cylindrical region on each t = constant 
hypersurface develops. The boundary of the cylindrical region is the internal 
infinity. It is a compact 2-dimensional spacelike surface. The hypersurfaces 
t = constant do not intersect the horizon but only approach such an inter- 
section at the internal infinity. (See for the conformal diagram and a de- 
tailed discussion, including analysis of the electrovacuum Robinson-Bertotti 
universe as the asymptotic limit of the extreme Reissner-Nordstrom geometry 
at the internal infinity.) 

There has been much interest in the extreme Reissner-Nordstrom black 
holes within standard Einstein-Maxwell theory. They admit surprisingly sim- 



ple solutions of the perturbation equations |117]. Some of them appear to be 



stable with respect to both classical and quantum processes, and there are 



attempts to interpret them as solitons |118 . Also, they admit supersymmetry 
TT9|. 



The quotation marks in the title of this section play a double role: the 
last two words are just "quoting" from the end of the title of a general re- 
view on string theory and supersymmetry prepared for the special March 
1999 issue of the Reviews of Modern Physics in honor of the centenary of 
the American Physical Society by Schwarz and Seiberg |120[, but they also 
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should "self-ironically" indicate my ignorance in these issues. In addition, 
unified theories of the type of string theory appear to be somewhat outside 
the direct interest of Jiirgen Ehlers, who has always emphasized the depth 
and economy of general relativity because it is a "background-independent" 
theory: string theories still suffer from the lack of a background-independent 
formulation. Nevertheless, they are beautiful, consistent, and very challeng- 
ing constructions, representing one of the most active areas of theoretical 
physics. Recently, string theory provided an explanation of the Bekenstein- 
Hawking prediction of the entropy of extreme and nearly extreme black holes. 
From the point of view of this review we should emphasize that many of the 
techniques that have been used to obtain exact solutions - mostly exact black 
hole solutions - in generalized theories like string theory were motivated by 
classical general relativity. There are also results in classical general relativity 
which are finding interesting generalizations to string theories, as we shall see 
with one example below. 

Before making a few amateurish comments on new results concerning 
extreme black holes in string theories, let us point out that in many pa- 
pers from the last 20 years, black hole solutions were studied in spacetimes 
with the number of dimensions either lower or higher than four. The lower- 
dimensional cases are usually analyzed as "toy models" for understanding the 
complicated problems of quantum gravity. The higher-dimensional models 
are motivated by efforts to find a theory which unifies gravity with the other 
forces. The most surprising and popular (2+l)-dimcnsional black hole is the 
BTZ (Baiiados-Teitelboim-Zanelli) black hole in the Einstein theory with a 
negative cosmological constant. Locally it is isometric to anti de Sitter space 
but its topology is different. In |121| the properties of (2+l)-dimcnsional 
black holes are reviewed. In (H-l)-dimcnsions one obtains black holes only if 
one includes at least a simple dilaton scalar field; the motivation for how to 
do this comes from string theory. In higher dimensions one can find general- 
izations of all basic black hole solutions in four dimensions ^122| . Interesting 
observations c once rning higher-dimensional black holes have been given a 



few years ago [123]. Perhaps one does not need to quantize gravity in order 



to remove the singularities of classical relativity. It may well be true that 
some new classical physics intervenes below Planckian energies. In |123| it is 
demonstrated that certain singularities of the four-dimensional extreme dila- 
ton black holes can be resolved by passing to a higher-dimensional theory of 
gravity in which usual spacetime is obtained only below some compactifica- 
tion scale. A useful, brief pedagogical introduction to black holes in unified 



theories is contained in |76 



One of the most admirable recent results of string theory, which un- 
doubtedly converted some relativists and stimulated many string theorists, 
has been the derivation of the exact value of the entropy of extreme and 
nearly extreme black holes. I shall just paraphrase a few statements from the 
March 1999 review for the centenary of the American Physical Society by 
Horowitz and Teukolsky |124]. There are very special states in string the- 
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ory called BPS (Bogomorny-Prasad-Sommerfield) states which saturate an 
equahty M > c|Q|, with M being the mass, Q the charge, and c is a fixed 
constant. The mass of these special states does not get any quantum cor- 
rections. The strength of the interactions in string theory is determined by 
a coupling constant g. One can count BPS states at large Q and small g. 
By increasing g one increases gravity, and then all of these states become 
black holes. (The BPS states are supersymmetric and one can thus follow 
the states from weak to strong coupling.) But they all become identical ex- 
treme Reissner-Nordstrom black holes, because there is only one black hole 
for given M = \Q\. When one counts the number N of BPS states in which 
an extreme hole can exist, and compares this with the entropy Sbh = \A of 
the hole as obtained in black hole thermodynamics ||8^ , |l25| , where A is the 
area of the event horizon [A — AttM'^ for the extreme Reissner-Nordstrom 
black hole), one finds exactly the "classical" result: Sbh = logA^!. The en- 
tropy of the classical black hole configuration is given in terms of the number 
of quantum microstates associated with that configuration, by the basic for- 



mula of statistical physics. For more detailed recent reviews, see [126,127| 



and references therein. Remarkably, the results for the black hole entropy 
have been obtained also within the canonical quantization of gravity |128|. 
A comprehensive review [129| of black holes and solitons in string theory 
appeared very recently. 

Allow me to finish this "all that" section with a personal remark. In 1980 
L. Dvorak and I found that in the Einstein-Maxwell theory, external mag- 
netic fiux lines are expelled from the black hole horizon as the hole becomes 
an extreme Reissner-Nordstrom black hole [104|. Hence, extreme black holes 
exhibit some sort of "Meissner effect" known from superconductivity. Last 
year it was demonstrated by Chamblin, Emparan and Gibbons |130] that this 
effect occurs also for black hole solutions in string theory and Kaluza-Klcin 
theory. Other extremal solitonic objects in string theory (like p-branes) can 
also have superconducting properties. Within the Einstein-Maxwell theory 
this effect was first studied to linear order in magnetic field - we analyzed 
Reissner-Nordstrom black holes in the presence of magnetic fields induc ed by 
current loops. However, we also used an exact solution due to Ernst |131|, 
describing a charged black hole in a background magnetic field, which asymp- 
totically goes over to a Melvin universe, and found the same effect (see also 
[ [132| for the case of the magnetized Kerr-Ncwman black hole). In |13C| the 
techniques of finding exact solutions of Einstein's field equations are employed 
within string theory and Kaluza-Klcin theory to demonstrate the "Meissner 
effect" in these theories. 
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4 The Kerr metric 

The discovery of the Kerr metric in 1963 and the proof of its unique role in 
the physics of black holes have made an immense impact on the development 
of general relativity and astrophysics. This can hardly be more eloquently 
demonstrated than by an emotional text from Chandrasekhar [Q: "In my 
entire scientific life, extending over forty-five years, the most shattering expe- 
rience has been the realization that an exact solution of Einstein's equations 
of general relativity, discovered by the New Zealand mathematician Roy Kerr, 
provides the absolutely exact representation of untold numbers of massive 

black holes that populate the Universe..." 

In Boyer-Lindquist coordinates the Kerr metric |l33| looks as follows (see 
e.g. |8||0|): 



= - ( 1 — j dr - 2 dt dip + 

+ ^dr^ + i:de^+ ^sm^edif^, (20) 



where 



= + cos^ 6*, A = r"^ - 2Mr + , 

S{r^ +a^) + 2Mra^sm^ e, (21) 

where M and a are constants. 



4.1 Basic features 

The Boyer-Lindquist coordinates follow naturally from the symmetries of the 
Kerr spacetime. The scalars t and can be fixed uniquely (up to additive 
constants) as parameters varying along the integral curves of (unique) sta- 
tionary and axial Killing vector fields k and m; and the scalars r and can 
be fixed (up to constant factors) as parameters related as closely as possi- 
ble to the (geometrically preferred) principal null congruences, which in the 
Kerr spacetime exist (see e.g. p^ , pO| ), and their projections on to the two- 



dimensional spacelike submanifolds orthogonal to both k and m (see [134 
for details). The Boyer-Lindquist coordinates represent the natural general- 
ization of Schwarzschild coordinates. With a = the metric (|o|) reduces to 
the Schwarzschild metric. 

By examining the Kerr metric in the asymptotic region r — > oo, one finds 
that M represents the mass and J — Ma the angular momentum pointing 
in the z-direction, so that a is the angular momentum per unit mass. One 
can arrive at these results by considering, for example, the weak field and 
slow motion limit, M/r <C 1 and a/r ^ 1. The Kerr metric ( pO[ ) can then be 
written in the form 
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dt^ 



2M 



dr^ 



AaM 



■ sin 6 dip dt, 



(22) 



which is the weak field metric generated by a central body with mass M 
and angular momentum J = Ala. A general, rigorous way of interpreting the 
parameters entering the Kerr metric starts from the definition of multipole 
moments of asymptotically flat, stationary vacuum spacetimes. This is given 
in physical space by Thorne 1 135 , using his "asymptotically Cartesian and 
mass centered" coordinate systems, and by Hansen |136 , who, generalizing 



the definition of Geroch for the static case, gives the coordinate independent 
definition based on the conformal completion of the 3-dimensional manifold 
of trajectories of a timelike Killing vector k. The exact Kerr solution has 
served as a convenient "test-bed" for such definitions.]^ The mass monopole 
moment - the mass - is M , the mass dipole moment vanishes in the "mass- 
centered" coordinates, the quadrupole moment components are ^Ma^ and 
— |Ma^. The current dipole moment - the angular momentum - is nonvan- 
ishing only along the axis of symmetry and is equal to J = Ma, while the 
current quadrupole moment vanishes. All other nonvanishing Z-pole moments 
are proportional to Afa' [135,136|. Because these specific values of the mul- 
tipole moments depend on only two parameters, the Kerr solution clearly 
cannot represent the gravitational field outside a general rotating body. In 
Section 6.2 we indicate how the Kerr metric with general values of M and a 
can be produced by special disk sources. The fundamental significance of the 
Kerr spacetime, however, lies in its role as the only vacuum rotating black 
hole solution. 

Many texts give excellent and thorough discussions of properties of Kerr 
black holes from various viewpoints [PJl9||2^,|30|,^0|j75l|76|,^,P j9^,p^ . The 



Kerr metric entered the new edition of "Landau and Lifshitz" [139|. A few 
years ago, a book devoted entirely to the Kerr geometry appeared |14C|. Here 
we can list only a few basic points. 

As with the Reissner-Nordstrom spacetime, one can make the maximal 
analytic extension of the Kerr geometry. This, in fact, has much in common 
with the Reissner-Nordstrom case. Loosely speaking, the "repulsive" charac- 
ters of both charge and rotation have somewhat similar manifestations. When 



< M^, the metric (Ed) has coordinate singularities at Z\ = 0, i.e. at (cf. 



r 



r± = M ± {M^ 



a 



(23) 



For the most complete, rigorous treatment of the asymptotic structure of station- 
ary spacetimes characterized uniquely by mul tipole moments defined at spatial 



infinity, see the work by Beig and Simon [ L37 ] , the article by Beig and Schmidt 
in this volume, and references therein. 
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Fig. 6. A schematic picture of a Kerr black hole with two horizons, ergosphere, and 
the ring singularity; local light wavefronts are also indicated. Particle A, entering 

the ergosphere from infinity, can split inside the ergosphere into particles B and C 
in such a manner that C arrives at infinity with a higher energy than particle A 
came in. 

The intrinsic three-dimensional geometry at r = r± reveals that these are 
null hypersurfaces - the (outer) event horizon and the inner horizon (Fig. 
6). As with the Reissncr- Nordstrom metric, the inner horizon - the Cauchy 
hypersurface - is unstable (see a more detailed discussion in Section 3.2 and 
references there). And as in the Penrose diagram in Fig. 5, one finds infinitely 
many asymptotically flat regions in the analogous Penrose diagram for the 
Kerr black hole spacetime. 

A crucial difference between the Reissner-Nordstrom and Kerr geome- 
try is the existence of the ergosphere (or, more precisely, ergoregion) in the 
Kerr case. This is caused by the dragging of inertial frames due to a non- 
vanishing angular momentum. The timelike Killing vector k, given in the 
Boyer-Lindquist coordinates by d/dt, becomes null "sooner", at r = ro, than 
at the event horizon, tq > r+ a.t 9 ^ 0, tt, as a consequence of this dragging: 



This is the location of the ergosurface, the ergoregion being between this sur- 
face and the (outer) horizon. In the ergosphere, schematically illustrated in 
Fig. 6, the "rotating geometry" drags the particles and light (with wavefronts 
indicated in the figure) so strongly that all physical observers must corotate 
with the hole, and so rotate with respect to distant observers - "fixed stars" 



= -gtt = 1 - 2Mr/E = 0, 
r = ro = M+ (M^ - cos^ 0)^ . 



(24) 
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- at rest in the Boyer-Lindquist coordinates.^ Static observers, whose world- 
lines (r, 9, if) = constant would have k as tangent vectors, cannot exist since 
k is spacelike in the ergosphere. Indeed, a non-spacelike worldline with r, 6 
fixed must satisfy the condition 

gttdt'^ + gtptpdip^ + 2g^td(pdt < 0, (25) 

in which gu = —k^ka, g^^ — m^ma, gipt = k^rria are invariants. In the 
ergosphere, the metric (20), ( |2l| ) yields gu > 0, g^p,^ > 0,g^pt < 0, so that 



dip/dt > - an observer moving along a non-spacelike worldline must corotate 
with the hole. The effect of dragging on the forms of photon escape cones 
in a general Kerr field (without restriction < M^) has been numerically 



studied and carefully illustrated in a number of figures only recently [142|. 

In order to "compensate" the dragging, the congruence of "locally nonro- 
tating frames" (LNRFs), or "zero-angular momentum observers" (ZAMOS), 
has been introduced. These frames have also commonly been used outside 
relativistic, rapidly rotating stars constructed numerically, but the Kerr met- 
ric played an inspiring role (as, after all, in several other issues, such as in 
understanding the ergoregions, etc.). The four-velocity of these (not freely 
falling!) observers, given in Boyer-Lindquist coordinates by 



{A/X:A)^,0,0,2aMr/{ASA)^ , (26) 



is orthogonal to the hypersurfaces t = constant. The particles falling from 
rest at infinity with zero total angular momentum fall exactly in the radial 
direction in the locally nonrotating frames with an orthogonal triad tied to 
the r, 9, ip coordinate directions (see [ [143| for the study of the shell of such 
particles falling on to a Kerr black hole). 

Now going down from the ergosphere to the outer horizon, we find that 
both Killing vectors k and m are tangent to the horizon, and are spacelike 
there (with k "rotating" with respect to infinity). The null geodesic gen- 
erators of the horizon are tangent to the null vectors Z = fc + f2m, where 
fi = a/2Mr_|_ — constant is called the angular velocity of the hole, fi is con- 
stant over the horizon so that the horizon rotates rigidly. Since i is a Killing 
vector, the horizon is a Killing horizon (cf. Section 2.4). 

Another notable difference from the Reissner-Nordstrom metric is the 
character of the singularity at 17 = 0, i.e. at r = 0, 6' = 7r/2. It is timelike, 
as in the Reissner-Nordstrom case, but it is a ring singularity (see Fig. 6). 
In the maximal analytic extension of the Kerr metric one can go through 
the ring to negative values of the coordinate r and discover closed timelike 



It is instructive to analyze the somewhat "inverse problem" of gravitational col- 
lapse of a slowly rotating dust shell which produces the Kerr metric, linearized 
in a/Af, outside (cf. Eq.(|22[)), and has flat space inside. Fixed distant stars seen 
from the centre of such shell appear to r otate due to the dragging of inertial 



frames, as was discussed in detail recently |141 
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lines since g^p^p < there. If the Kerr parameters are such that > A'P, the 
Kerr metric does not represent a black hole. It describes the gravitational 
field with a naked ring singularity. The Kerr ring singularity has a repulsive 
character near the rotation axis. It gives particles outward accelerations and 



collimates them along the rotation axis [I44|, which might be relevant in 
the context of the formation and precollimation of cosmic jets. However, the 
cosmic censorship conjecture is a very plausible, though difficult "to prove" 
hypothesis, and Kerr naked singularities are unlikely to form in nature. How- 
ever, the Kerr geometry with > M^, with a region containing the ring 
singularity "cut out" , can be produced by thin disks; though if they should 
be composed of physical matter, they cannot be very relativistic (see Section 
6.2 and references therein). 

If = Af^, the Kerr solution represents an extreme Kerr black hole, as is 
the analogous Reissner-Nordstrom case with — AP. The inner and outer 
horizons then coincide at r = M . The horizon is degenerate with infinite afHne 
length. Almost extreme Kerr black holes probably play the most important 
role in astrophysics (see below). In realistic astrophysical situations accreting 
matter will very likely have a sufficient amount of angular momentum to turn 
a Kerr hole to an almost extreme state. 

There exists a charged, electrovacuum generalization of the Kerr fam- 
ily found by Newman et al. ]145[ |. The Kerr-Newman metric in the Boyer- 
Lindquist coordinates can be obtained from the Kerr metric ( pO[ ) if all of 
the terms 2Mr explicitly appearing in (|20|), ( pT| ) are replaced by 2Mr — Q^, 
with Q being the charge. The metric describes charged, rotating black holes 
Q^, with two horizons located at r± — M ± (M^ — — Q^)^ ■ 

+ Q"^, and with A'P < + Q"^ one 
obtains naked (ring) singularities. The analytic extension, the presence of 
ergoregions and the structure of the singularity is similar to the Kerr case. 

In addition to the gravitational field, there exists a stationary electro- 
magnetic field which is completely determined by the charge Q and rotation 
parameter a. The vector potential of this field is given by the 1-form 

Ao,dx" = ~{Qr/S){dt - asin^ 9 dip), (27) 

so that if a 7^ the electric field is supplemented by a magnetic field. At large 
distances {r oo) the field corresponds to a monopole electric field with 
charge Q and a dipole magnetic field with magnetic moment n = Qa. Since 
the gyromagnetic ratio of a charged system with angular momentum J is 
defined by 7 = /i/ J, one finds the charged-rotating-black hole gyromagnetic 
ratio to satisfy the relation 7 — Q/M, i.e. it is twice as large as that of 
classical matter, and the same as that of an electron. By examining a black 
hole with a loop of rotating charged matter around it, the radius of the loop 
changing from large values to the size of the horizon, it is possible to gain 



some understanding of this result |146 
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4.2 The physics and astrophysics around rotating black holes 

In the introduction to their new 770 page monograph on black hole physics, 
Frolov and Novikov |7^ write: "... there are a lot of questions connected 
with black hole physics and its applications. It is now virtually impossible to 
write a book where all these problems and questions are discussed in detail. 
Every month new issues of Physical Review D, Astrophysical Journal, and 
other physical and astrophysical journals add scores of new publications on 
the subject of black holes..." Although Frolov and Novikov have also black 
hole-like solutions in superstring and other theories on their minds, we would 
not probably be much in error, in particular in the context of astrophysics, if 
we would claim the same just about Kerr black holes. Hence, first of all, we 
must refer to the same literature as in the previous section 4.1. A few more 
references will be given below. 

A remarkable fact which stands at the roots of these developments is that 
the wave equation is separable, and the geodesic equations are integrable in 
the Kerr geometry. Carter |7^ , who explicitly demonstrated the separability 
of the Hamilton- Jacobi equation governing the geodesic motion, has empha- 
sized that one can in fact derive the Kerr metric as the simplest nonstatic 
generalization of the Schwarzschild solution, by requiring the separability of 
the covariant Klein-Gordon wave equation.p| 

A thorough and comprehensive analysis of the behaviour of freely falling 
particles in the Kerr field would produce material for a book. We refer to e.g. 



[18 76,91,134 138 , 144 ,1471 for fairly detailed accounts and a number of further 
references. From the point of view of astrophysical applications the following 
items appear to be most essential: in contrast to the Schwarzschild case, 
where the stable circular orbits exist only up to r = 3r+ = 6M, in the field of 
rotating black holes, the stable direct (i.e. with a positive angular momentum) 
circular orbits in the equatorial plane can reach regions of "deeper potential 
well". With an extreme Kerr black hole the last stable direct circular orbit 



occurs at r = r_|_ = M. (See [147| for a clear discussion of the positions of the 
innermost stable, innermost bound, and photon orbits as the hole becomes 
extreme and a long cylindrical throat at the horizon develops.) A "spin-orbit- 
coupling" effect increases the binding energy of the direct orbits and decreases 
the binding energy of the retrograde (with a negative angular momentum) 
orbits relative to the Schwarzschild values. The binding energy of the last 
stable direct circular orbit is AE = 0.0572/i (^ is the particle's proper mass) 
in the Schwarzschild case, whereas AE — 0.4235/i for an extreme Kerr hole. 
A particle slowly spiralling inward due the emission of gravitational waves 



Although this is still not a "constructive, analytic derivation of the Kerr metri c 
which would fit its physical meaning", as required by Landau and Lifshitz |l39[ |, 
it is certainly more intuitive than the original derivation by Kerr. On the other 
hand, despite various hints like the existence of the Killi ng t ensor field (in addition 



to the Killing vectors) in the Kerr geometry (see e.g. |134]), it does not seem to 



be clear why the Kerr geometry makes it possible to separate these equations. 
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would radiate the total energy equal to this binding energy; hence much more 
- 42% of its rest energy - in the Kerr case. The second significant effect is 
the dragging of the particles moving on orbits outside of the equatorial plane. 
The draggin(jf^ will make the orbit of a star around a supermassive black hole 
to precess with angular velocity ^ 2J/r^. The star may go through a disk 



around the hole, subsequently crossing it at different places 150 . One can 
also show that as a result of the joint action of the gravomagnetic effect and 
the viscous forces in an accretion disk, the disk tends to be oriented in the 
equatorial plane of the central rotating black hole (the "Bardeen-Petterson 
effect"). 

The above examples demonstrate specific effects in the Kerr background 
which very likely play a significant role in astrophysics (see also Section 4.3 
below). The best known process in the field of a rotating black hole is probably 
astrophysically unimportant, but is of principal significance in the black hole 
physics. This is the Penrose process for extracting energy from rotating black 
holes. It is illustrated schematically in Fig. 6: particle A comes from infinity 
into the ergosphere, splits into two particles, B and C. Whereas C is ejected 
back to infinity, B falls inside the black hole. The process can be arranged 
in such a way that particle C comes back to infinity with higher energy than 
with which particle A was coming in. The gain in the energy is caused by the 
decrease of rotational energy of the hole. Such process is possible because the 
Killing vector k becomes spacelike in the ergosphere, so that the (conserved) 
energy of particle B "as measured at infinity" (see e.g. |Q), Eb = —kaP%, 
can be negative. Unfortunately, the "explosion" of particle B requires such a 
big internal energy that the process is not realistic astrophysically. 

More general considerations of the interaction of black holes with matter 
outside have led to the formulation of the four laws of black hole thermody- 
namics |p^ , [76t^2[|l25|| . These issues, in particular after the discovery of the 
Hawking effect that black holes emit particles thermally with temperature 
T = Kh/2'Kkc («;-surface gravity, /c-Boltzmann's constant), have been an in- 
spiration in various areas of theoretical physics, going from general relativity 



and statistical physics, to quantum gravity and string theory (see |125,126 
and some remarks and references in Section 3.2). The Kerr solution played 
indeed the most crucial role in these developments. I recall how during my 
visits to Moscow in the middle of the 1970s Zel'dovich and his colleagues were 
somewhat regretfully admitting that they were on the edge of discovering the 
Hawking effect. They realized that an analogue of the Penrose process occurs 
with the waves (in so called superradiant scattering) which get amplified if 
their energy per unit angular momentum is smaller than the angular velocity 



Relativists often consider the effects produced by moving mass currents as "the 
dragging of inertial frames", but the concept of the gravomagnetic field, or gravo- 



magnetism has some advantages, as has been stressed recently [146]. The gravo- 
magnetic viewpoint, however, has also been used in many works in the past - 
see, e.g. ||7l|,[75[ , and in particular |14£|] and references therein. 
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f2 of the horizon. Zel'dovich then suggested that there should be spontaneous 
emission of particles in the corresponding modes but did not study quantum 
fields on a nonrotating background (cf. an account of these developments, 
including the visit of Hawking to Moscow in 1973, by Israel Jzo]). 

Returning back to Earth or, rather, up to heavens, it is not so well known 
that an astrophysically more realistic example of the Penrose process exists: 
this is the Blandford-Znajek mechanism - see [ |75| , p8| , p"5l| - in which a mag- 
netic field threading a rotating hole (the field being maintained by external 
currents in an accretion disk, for example) can extract the hole's rotational 
energy and convert it into a Poynting flux and electron-positron pairs. A 
Kerr black hole with angular momentum parallel to an external magnetic 
field acts (by "unipolar induction") like a rotating conductor in an external 
field. There will be an induced electric field and a potential difference be- 
tween the pole and the equator. If these are connected, an electric current 
will fiow and power will be dissipated. In fact, this appears to be until now 
the most plausible process to explain gigantic relativistic jets emanating from 
the centres of some of the most active galaxies. The i?Z-mechanism has its 
problems: extremely rotating black holes expel magnetic flux ||75,152| - there 
probably exists a value of the angular momentum Jq < Jmax for which the 
power extracted will be g reatest. It is not clear whether the process can be ef- 
flciently maintained |153|; and perhaps more importantly, new astrophysical 
estimates of seed magnetic flelds seem to be too low to make the mechanism 
efficient |154]. The i?Z-mechanism will probably attract more attention in 
the coming years, in particular in view of the recent discovery of two "mi- 
croquasars" in our own Galaxy, which generate double radio structures and 
apparent superluminal jets similar to extragalactic strong radio sources |155|. 

A remarkable achievement of pure mathematical physics, with a great 
impact on astrophysics, has not only been the discovery of the Kerr solution 
itself but also the development of the theory of Kerr metric perturbations 
[[75y7q,|91,156|. By employing the Newman-Penrose null tetrad formalism, in- 



vented and extensively used in mathematical relativity, in particular in grav- 
itational radiation theory, it has been possible to separate completely all 
perturbation equations for non-zero spin fields. In particular, a single "mas- 
ter equation" - called the Teukolsky equation - governs scalar, electromag- 
netic and gravitational perturbations of a Kerr black hole.^^ If no sources are 
present on the right-hand side, the equation looks as follows: 



(r^ + 0^)2 
A 



2 • 2 

a sm I 



dr\ dr I sin 9 89 
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In the case of a Kerr-Newman black hole, the electromagnetic and gravitational 
perturbations necessarily couple. Until now, in contrast to the spherical Reissner- 
Nordstrom case, a way of how to decouple them has not been discovered. 
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-2s 



r — la cos f 



dtp 
'dt 



+ (s^ cot^ e ~ s)'ijj = 0. (28) 



The coordinates are the Boyer-Lindquist coordinates used in (pQ), A is de- 
fined in (^l|), and s is the spin weight of the perturbing field; s = 0, ±1, ±2. 
The variables in the Teukolsky equation can be separated by decomposing -0 
according to 



i'lrn = {l/V2T:)sRlmir,Uj)sSlrniO)e 



iniif —iLut 



(29) 



where sSim are so called spin weighted spheroidal harmonics. By solving the 
radial Teukolsky equation for sRim with appropriate boundary conditions one 
can find answers to a number of (astro)physical problems of interest like the 
structure of stationary electromagnetic or gravitational fields due to external 
sources around a Kerr black hole (e.g. [75 157|), the emission of gravitational 
waves from particles plunging into the hole (e.g. | 7G|,p7 15S]), or the scattering 
of the waves from a rotating black hole (e.g. ||76,15Gj and references therein). 
At present, the Teukolsky equation is being used to study the formation of a 
rotating black hole from a head-on collision of two holes of equal mass and 
spin, initially with small separation, to fi nd the wave forms of gravitational 
radiation produced in this process [159|. The first studies of second-order 
perturbations of a Kerr black hole are also appearing [ 160 1 . 

To find all gravitational (metric) perturbations by solving the complete 
system of equations in the Newman-Penrose formalism is in general a formi- 
dable task. As Chandrasekhar's "last observation" at the end of his chapter on 
gravitational perturbations of the Kerr black hole reads |Q : "The treatment 
of the perturbations of the Kerr spacetime in this chapter has been prolixious 
in its complexity. Perhaps, at a later time, the complexity will be unravelled 
by deeper insights. But mean time, the analysis has led us into a realm of 
the rococo: splendorous, joyful, and immensely ornate." 



4.3 Astrophysical evidence for a Kerr metric 

Very recently new observations seem to have opened up real possibilities of 
testing gravity in the strong-field regime. In particular, it appears feasible 
to distinguish the Kerr metric from the Schwarzschild, i.e. to measure a/M. 
Our following remarks on these developments are based on the review by M. 



Rees [ p8[ , and in particular, on the very recent survey by A. Fabian |161|, 
an authority on diagnosing relativistic rotation from the character of the 
emission lines of accretion disks around black holes. 

The interest here is not in optical lines since the optical band comes from 
a volume much larger than the hole. However, the X-rays are produced in 
the innermost parts of an accretion fiow, and should thus display substantial 
gravitational redshifts as well as Doppler shifts. This only became possible 
to observe quite recently, when the ASCA X-ray satellite started to operate. 
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and the energy resolution and sensitivity became sufficient to analyze line 
shapes. 

Typically the profile of a line emitted by a disk from gas orbiting around 
a compact object has a double-horned shape. The disk can be imagined to 
be composed of thin annuli of orbiting matter - the total line is then the 
sum of contributions from each annulus. If the disk is not perpendicular to 
our line of sight its approaching sides will - due to classical Doppler shifts - 
produce blue peaks, receding sides red peaks. The broadest parts of the total 
line come from the innermost annuli because the motion there is fastest. In 
addition, there are relativistic effects: they imply that the emission is beamed 
in the direction of motion, transverse Doppler shifted, and gravitationally 
redshifted. As a result, the total line is broad and skewed in a characteristic 
manner. Such lines are seen in the X-ray spectra of most Seyfert 1 galaxies. 
In the Seyfert galaxy MCG-6-30-15 the fluorescent iron line was observed to 
be (red)shifted further to lower energies]^ 
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Fig. 7. The broad iron line from MCG-6-30-15. The best-fitting, maximally spin- 
ning Kerr black hole model is shown (from Iwasawa, K. et al. (1996), Mon. Not. 
Roy. Astron. Soc. 282, 1038). 



This suggests that the emission took place below 3iZs = 6Af , i.e. below the 
innermost stable orbit for a Schwarzschild black hole. In 1996, the line shape 
was well fitted by the assumption that the line is produced in a close orbit 

In Seyfert 1 galaxies hard flares occur which irradiate the accretion disk, and 
produce a reflection component of continuum peaking at ~ 30 keV and the fluo- 
rescent iron line at about 6.5 keV. 
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around maximally rotating (extreme) Kerr black hole. In 1997, the param- 
eter a/M was quantified as exceeding 0.95. Hence, it has been "tentatively 
concluded that the line was the first spectroscopic evidence for a Kerr hole" 

There are alternative models for a broad skew iron line, including Comp- 
tonization by cold electrons, or the emission from irradiated matter falling 
from the inner edge of the disk around a nonrotating Schwarzschild black hole. 
It appears, however, that the data speak against these possibilities |161|. In 
any case, with future X-ray detectors, which will yield count rates orders 
of magnitude higher than ASCA, the line shapes should reveal in a much 
greater detail specific features of the Kerr metric. 

In addition, other possibilities to determine a/M exist. These include: 

(i) Observations of stars in relativistic orbits going through a disk around a 
supermassive rotating black hole ||88,150|. 

(ii) Characteristic frequencies of the vibrational modes in disks or tori around 
rotating black holes [ B8 162| . 

(iii) The precession of a disk which is tilted with respect to the hole's spin 
axis. This precession arises because of frame dragging and produces a periodic 
modulation of the X-ray luminosity. 

(iv) Astrophysically most important would be a discovery showing that the 
properties of cosmic jets depend on the value of a/M . This could indicate 
that the Blandford-Znajek mechanism (see section 4.2) is really going on. Its 
likelihood would increase if jets were found with Lorentz factors 7 significantly 
exceeding 10 (see for more details). 

(v) Last but not least, future observations of gravitational waves from black 
hole collisions offer great hopes of a clean observation of a black hole 
geometry, without astrophysical complications. 

It is hard to point out any other exact solution of Einstein's field equations 
(or of any kind of field equations?) discovered in the second half of the 20th 
century which has had so many impacts on so many diverse areas of physics, 
astrophysics, astronomy, and even space science as has had the Kerr metric. 



5 Black hole uniqueness and multi-black hole solutions 

Since black holes can be formed from the collapse of various matter configu- 
rations, it is natural to expect that there will be many solutions of Einstein's 
equations describing black holes. It is expected that the asymptotic final 
state of a collapse can be represented by a stationary spacetime, i.e. one 
which admits a 1-dimensional group of isometrics whose orbits are timelike 
near infinity. Strong arguments show that the event horizon of a station- 
ary black hole must be a Killing horizon. One of the most remarkable and 
surprising results of black hole theory are the sequence of theorems showing 
rigorously that the only stationary solution of the Einstein electrovacuum 
equations that is asymptotically flat and has a regular event horizon is the 
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Kerr-Newman solution. There is a number of papers on this issue - recent 
detailed reviews are given in pO| , ^ . The intuition gained from exact black 
hole solutions in proving the theorems has been essential. 

Roughly speaking, the uniqueness proof consists of the following three 
parts. First, one demonstrates the "rigidity theorem" , which claims that non- 
degenerate (k ^ 0) stationary electrovacuum analytic black holes are either 
static or axially symmetric. One then establishes that the Reissner- Nordstrom 
nondegenerate electrovacuum black holes are all static (nonrotating) nonde- 
generate black holes in electrovacuum. Finally, one separately proves that the 
nondegenerate Kerr-Newman black holes represent all nondegenerate axially 
symmetric stationary electrovacuum black holes. 

Although such results were proved more than 10 years ago, recently there 
has been new progress in the understanding of the global structure of station- 
ary black holes. Again, exact solutions have been inspiring: by gluing together 
two copies of the Kerr spacetime in a certain way, Chrusciel pC| ] constructed a 
black hole spacetime which is stationary but not axisymmetric, demonstrat- 
ing thus that the standard formulation and proof of the rigidity theorem [p6| 
is not correct. (The reason being essentially that when one extends the isome- 
tries from a neighbourhood of the horizon by analytic continuation one has no 
guarantee that the maximal analytic extension is unique.) Chrusciel proved 
"a corrected version of the black hole rigidity theorem" ; in the connected 
case one can prove a uniqueness theorem for static electrovacuum black holes 
with degenerate horizons. The uniqueness theorem for static degenerate black 
holes which demonstrates that the extreme Reissner-Nordstrom black hole is 
the only case, is of importance also in string theory. The most unsatisfactory 
feature of the rigidity theorem is the assumption of analyticity of the metric 
in a neighborhood of the event horizon. In this context, Chrusciel men- 
tions the case of Robinson- Trautman exact analytic metrics, which can be 



smoothly but not analytically extended through an event horizon [163 . We 
shall discuss this issue in somewhat greater detail in Section 10. 

The black hole uniqueness theorems indicated above are concerned with 
only single black holes. (Corresponding spacetimes contain an asymptotically 
flat hypersurface S with compact interior and compact connected boundary 
dS which is located on the event horizon.) Consequently a question naturally 
arises as to whether one can generalize the theorems to some multi-black 
hole solutions. In classical physics a solution exists in which a system of n 
arbitrarily located charged mass points with charges qi and masses vrii , such 
that \qi\ = \fGmi^ is in static equilibrium. In relativity the metric 



ds^ = -y-'^dt^ + V'^ {dx^ + dy^ + d??) , (30) 
with time-independent V satisfying Laplace's equation 

2 d'^V d'^V d'^V 
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is a solution of the Einstein-Maxwell equations with the electric field 

E = VV-\ (32) 

where V = {d/dx, d/dy, d/dz). (In standard units E — \/GWV~^ .) The sim- 
plest solution of this form is the Majumdar-Papapetrou metric, corresponding 
to a linear combination of n "monopole sources" with masses > and 
charges qi = rrii, located at arbitrary points xt: 

n 

y = l + V , , . (33) 



Hartle and Hawking |164| have shown that every such spacetime can be an- 
alytically extended to a spacetime representing n degenerate charged black 
holes in static equilibrium. The points x = Xi are actually event horizons of 



area Attiti^. For the case of one black hole, the metric (30) is just the extreme 
Reissner- Nordstrom black hole in isotropic coordinates. 

A uniqueness theorem for the Majumdar-Papapetrou metrics is not avail- 



able, although some partial answers are known (see [165 166 for more de- 
tails). It is believed that these are the only asymptotically flat, regular multi- 
black hole solutions of the Einstein-Maxwell equations. In fact, such a result 
would exclude an interesting possibility that a repulsive gravitational spin- 
spin interaction between two (or more) rotating, possibly charged, black holes 
can overcome their gravitational attraction and thus that there exists in Ein- 
stein's theory of gravitation ~ in contrast to Newton's theory - a stationary 
solution of the two-body problem. 

Among new solutions discovered by modern generating techniques, there 
are the solutions of Kramer and Neugebauer which represent a nonlinear su- 



perposition of Kerr black holes (see 167 for a review). These solutions have 



been the subject of a number of investigations which have shown that spin- 
spin repulsion is not strong enough to overcome attraction. In particular, 
two symmetrically arranged equal black holes cannot be in stationary equi- 
librium. The situation might change if one considers two Kerr-Newman black 



holes [168|. Here one has four forces to reckon with: gravitational and electro- 
magnetic Coulomb-type interactions, and gravitational and electromagnetic 
spin-spin interactions. One can then satisfy the conditions which render the 
system of two Kerr-Newman black hole free of singularities on the axis, and 
make the total mass of the system positive. However, there persists a singu- 



larity in the plane of symmetry away from the axis [ 168 1 . In view of this result 
we have conjectured that even with electromagnetic forces included one can- 
not achieve balance for two black holes, except for the exceptional case of two 
nonrotating extreme (degenerate) Reissner-Nordstrom black holes. Recently, 
some new rigorous results concerning the (non) existence of multi-black hole 



stationary, axisymmetric electrovacuum spacetimes have been obtained [169| 



(see also 80 ), but the "decisive theorem" is still missing. 
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In connection with the problem of the balance of gravity by a gravitational 
spin-spin interaction, we should mention that there exists the solution of Dietz 
and Hoenselaers [17C| in which balance of the two rotating "particles" is 



achieved. However, the "sources" are complicated naked singularities which 
become Curzon-Chazy "particles" (see Section 6.1) if the rotation goes to 
zero, and it is far from clear whether appropriate physical interior solutions 
can be constructed. 



In 1993, Kastor and Traschen |171] found an interesting family of solutions 
to the Einstein-Maxwell equations with a non-zero cosmological constant 
A. They describe an arbitrary number of charged black holes in a "back- 
ground" de Sitter universe. In the limit of yl = these solutions become 
Majumdar-Papapetrou static metrics. In contrast to these metrics, the cos- 
mological multi-black hole solutions with A > are dynamical. Remarkably, 
one can construct solutions which describe coalescing black holes. In some 
cases cosmic censorship is violated - a naked singularity is formed as a result 
of the collision |172| . Although these solutions do not have smooth horizons, 
the singularities are mild, and geodesies can be extended through them. The 
metric is always at least C^. Since the solutions are dynamical, one may in- 
terpret the non-smoothness of the horizons as a consequence of gravitational 
and electromagnetic radiation. In this sense, the situation is analogous to the 
case of the Robinson- Trautman spacetimes discussed in Section 10. In five or 
more dimensions, however, one can construct static multi-black hole solutions 



with A = 0, which do not have smooth horizons |173|. The solutions of Kas- 
tor and Traschen also inspired a new and careful analysis |174| of the global 
structure of the Reissner-Nordstrom-de Sitter spacetimes characterized by 
mass, charge, and cosmological constant. The structure is considerably richer 
than that with A = 0. Most recently, the hoop conjecture (giving the crite- 
rion as to whether a black hole forms from a collapsing system) was discussed 



[ |175| by analyzing the solution of Kastor and Traschen. 

6 On stationary axisymmetric fields and relativistic 
disks 

6.1 Static Weyl metrics 

The static axisymmetric vacuum metrics in Weyl's canonical coordinates p £ 
[0, oo), z,t £ lR,(p £ [0, 2n) have the form 

The bmction U{p, z) satisfies flat-space Laplace's equation 

d^U 1 dU d^U 



dp^ p dp dz"^ 



0. (35) 
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The function z) is determined from U by quadrature up to an addi- 
tive constant. The axis p = is free of conical singularities at places where 
limp^o k = 0. 

The mathematically simplest example is the Curzon-Chazy solution in 
which U = —m/ y'p^ + z^ is the Newtonian potential of a spherical point 
particle. The spacetime, however, is not spherically symmetric. In fact, one 
of the lessons which one has learned from this solution is the directional char- 
acter of the singularity at + = 0. For example, the limit of the invariant 
Ra/s-ysR"^''^ depends on the direction of approach to the singularity. The 
singularity has a character of a r ing through which some timelike geodesies 
may pass to a Minkowski region [ 176 1 . 

Various studies of the Weyl metrics indicated explicitly how important it 
is always to check whether a result is not just a consequence of the choice 
of coordinates. There is the subclass of Weyl metrics generated by the New- 
tonian potential of a constant density line mass ("rod") with total mass M 
and (coordinate) length I, which is located along the z-axis with the middle 
point at the origin. These are Darmois-Zipoy-Vorhees metrics, called also the 
7-metrics The Schwarzschild solution (a spherically symmetric metric!) 
is a special case in this subclass: it is given by the potential of the rod with 
/ — 2M. Clearly, in general there is no correspondence between the geometry 
of the physical source and the geometry of the Newtonian "source" from the 
potential of which a Weyl metric is generated. 

A survey of the best known Weyl metrics, including some specific solutions 
describing fields due to circular disks is contained in |64|. More recently, 
Bicak, Lynden-Bell and Katz |177| have shown that most vacuum static Weyl 
solutions, including the Curzon and the Darmois-Vorhees-Zip oy so lutions, 
can arise as the metrics of counterrotating relativistic disks (see 1 177 also for 
other references on relativistic disks). The simple idea which inspired their 
work is commonly used in Newtonian galactic dynamics |6|: imagine a point 
mass placed at a distance b below the centre p = of a plane z = 0. This gives 
a solution of Laplace's equation above the plane. Then consider the potential 
obtained by reflecting this z > potential in z = so that a symmetrical 
solution both above and below the plane is obtained. It is continuous but 
has a discontinuous normal derivative on 2 = 0, the jump in which gives 
a positive surface density on the plane. In galactic dynamics one considers 
general line distributions of mass along the negative z-axis and, employing 
the device described above, one finds the potential-density pairs for general 
axially symmetric disks. In |17S], an infinite number of new static solutions 
of Einstein's equations were found starting from realistic potentials used to 



describe flat galaxies, as given recently by Evans and de Zeeuw |179 



Although these disks are Newtonian at large distances, in their central 
regions interesting relativistic features arise, such as velocities close to the 
velocity of light, and large redshifts. In a more mathematical context, some 
particular cases are so far the only explicit examples of spacetimes with a 
"polyhomogeneous" null infinity (cf. [180| and Section 9), and spacetimes 
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with a meaningful, but infinite ADM mass |178]. New Weyl vacuum solutions 
generated by Newtonian potentials of flat galaxies correspond to both finite 
and semi- infinite rods, with the line mass densities decreasing according to 
general power laws. It is an open question what kinds of singularities rods 
with different density profiles represent. 

Very recently, new interesting examples of the static solutions describing 
self-gravitating disks or rings, and disks or rings around static black holes 
have been constructed [181 182|,|183|| and the effects of the fields on freely 
moving test particles studied il81[| . Exac t dis ks with electric currents and 
magnetic fields have also been considered |184|. 

Employing the Weyl formalism, one can describe nonrotating black holes 
strongly distorted by the surrounding matter. The influence of the matter can 
be so strong that it may even cause the horizon topology to be changed from 
spherical to toroidal (see |7^ and references therein). 

Finally, we have to mention two solutions in the Weyl class, which were 
found soon after the birth of general relativity, and have not lost their in- 
fluence even today. The flrst, discovered by Bach and Weyl, is assigned by 
Bonnor ||6^ as "probably the most perspicacious of all exact solutions in 
GR" . It refers to two Curzon-Chazy "monopoles" on the axis of symmetry. 
One flnds that the metric function k has the property that limp^o k ^ 0, so 
that there is a stress described by a conical singularity between the particles, 
which holds particles apart. A similar solution can be constructed for the 
Schwarzschild "particles" (black holes) held apart by a stress. These cases 
can serve as one of the simplest demonstrations of the difference between 
the Einstein theory and held theories like the Maxwell theory: it is only in 
general relativity in which fleld equations involve also equations of motion. 

The second "old" solution which has played a very significant role is the 
metric discovered by Levi-Civita. It belongs to the class of degenerate (type 
D) static vacuum solutions which form a subclass of the Weyl solutions. In the 
invariant classification of the degenerate solutions by Ehlers and Kundt |^ , 
this solution is contained in the last, third subclass. That is why Ehlers and 
Kundt called it the C-metric, and it is so well known today. We shall discuss 
the C-metric later (Section 11) in greater detail since, as it has been learned in 
the 1970s, it is actually a radiative solution representing uniformly accelerated 
black holes. What Levi-Civita found and Ehlers and Kundt analyzed is only 
a portion of spacetime in which the boost Killing vector is timelike, and 
the coordinates can thus be found there (analogous to the coordinates in 
a uniformly accelerated frame in special relativity) in which the metric is 
time- independent . 



6.2 Relativistic disks as sources of the Kerr metric and other 
stationary spacetimes 

Thanks to the black hole uniqueness theorems (Section 5), the Kerr metric 
represents the unique solution describing all rotating vacuum black holes. 
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Nevertheless, although the cosmic censorship conjecture, on which the phys- 
ical relevance of the Kerr metric rests, is a very plausible hypothesis, it re- 
mains, as was noted in several places above, one of the central unresolved 
issues in relativity. It would thus support the significance of the Kerr metric 
if a physical source were found which produces the Kerr field. The situation 
would then resemble the case of the spherically symmetric Schwarzschild 
metric which can represent both a black hole and the external field due to 
matter. 

This has been real ized by many workers. The review on the "Sources for 
the Kerr Metric" [185| written in 1978, contains 71 references, and concludes 
with: "Destructive statements denying the existence of a material source for 
the Kerr metric should be rejected until (if ever) they are reasonably justi- 
fied." The work from 1991 gives "a toroidal source", consisting of "a toroidal 
shell . . . , a disk . . .and an annulus of matter interior to the torus" |186|. 
The masses of the disk and annulus are negative. To summarize in Hermann 
Bondi's way, the sources suggested for the Kerr metric have not been the 
easiest materials to buy in the shops . . . 

The situation is somewhat different in the special case of the extreme 
Kerr metric, where there is a definite relationship between mass and angular 
momentum. The numerical study |187| of uniformly rotating disks indicated 
how the extreme Kerr geometry forms around disks in the "ultrarelativistic" 
limit. These numerical results have been supported by important analytical 
work (see Section 6.3). However, in the case of a general Kerr metric physical 
sources had not been found before 1993. 

A method similar to that of constructing disk sources of static Weyl space- 
times (described in Section 5.1) has been shown to work also for axisymmet- 
ric, reflection symmetric, and stationary spacetimes [ 188, 189 1. It is important 
to realize that although now no metric function solves Laplace's equation as 
in the static case, we may view the procedure described in Section 5.1 as the 
identification of the surface z = b with the surface z = —b. The field then 
remains continuous, but the jump of its normal derivatives induces a matter 
distribution in the disk which arises due to the identification of the surfaces. 
What remains to be seen, is whether the material can be "bought in the 
shops" . This idea can be employed for all known asymptotically fiat station- 
ary vacuum spacetimes, for example for the Tomimatsu-Sato solutions, for 
the "rotating" Curzon solution, or for other metrics (cf. |^ for references). 

Any stationary axisymmetric vacuum metric can be written in canonical 
coordinates {t, (p, p, z) in the form [ |6ll 

where [/, k, and A are functions of p, z. For the Kerr solution (mass M, 
specific angular momentum a > 0), the functions U, k,A are ratios of poly- 
nomials when expressed in spheroidal coordinates | |6l| . 

Now, identify the "planes" z ~ b = constant > and z — ~b (this iden- 
tification leads to disks with zero radial pressure). With the Kerr geometry 
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the matching is more complicated than in the static cases, and therefore, 
one has to turn to Israel's covariant formalism (see |19C] for its recent exposi- 
tion) . Using this formalism one is able to link the surface stress-energy tensor 
of the disk arising from this identification, to the jump of normal extrinsic 
curvature across the timelike hypersurface given hy z — b (with the jump be- 
ing determined by the discontinuities in the normal derivatives in functions 
U,k,A). 

The procedure leads to physically plausible disks made of two streams 
of coUisionless particles , that circulate in opposite directions with differen- 
tial velocities 1 188, 189 1. Although extending to infinity, the disks have fi- 
nite mass and exhibit interesting relativistic properties such as high veloci- 
ties, large redshifts, and dragging effects, including ergoregions. Physical disk 
sources of Kerr spacetimes with > AP can be constructed (though these 
are "less relativistic"). And the procedure works also for electrovacuum sta- 
tionary spacetimes. The disks with electric current producing Kerr-Newman 
spacetimes are described in |191|, where the conditions for the existence of 
(electro)geodesic streams are also discussed. 

The power and beauty of the Einstein field equations is again illustrated: 
the character of exact vacuum fields determines fully the physical character- 
istics of their sources. In a more sophisticated way, this is seen in the problem 
of relativistic rigidly (uniformly) rotating disks of dust. 



6.3 Uniformly rotating disks 

The structure of an infinitesimally thin, finite relativistic disk of dust par- 
ticles which rotate uniformly around a common centre was first explored 



by J. Bardeen and R. Wagoner 25 years ago |187]. By developing an effi- 
cient expansion technique in the quantity 5 = Zc/(1 + Zc), Zc denoting the 
central redshift, they obtained numerically a fairly complete picture of the 
behaviour of the disk, even in the ultrarelativistic regime ((5 — > 1). In their 
first letter from 1969 they noted that "there may be some hope of finding an 
analytic solution". Today such a hope has been substantiated, thanks to the 



work of G. Neugebauer and R. Meinel (see |192,193 and references therein). 
The solution had, in fact, to wait until the "soliton-type-solution generating 
techniques" for nonlinear partial differential equations had been brought over 
from applied mathematics and other branches of physics to general relativity, 
starting from the end of the 1970s. 

These techniques have been mainly applied only in the vacuum cases so 
far, but this is precisely what is in this case needed: the structure of the 
thin disk enters the field equations only through the boundary conditions at 
z = 0, < p < a (a is the radius of the disk). The specific procedures which 
enabled Neugebauer and Meinel to tackle the problem are sophisticated and 
lengthy. Nevertheless, we wish to mention them telegraphically at least, since 
they represent the first example of solving the boundary value problem for a 
rotating object in Einstein's theory by analytic methods. 



60 Jin Bicak 



In the stationary axisymmetric case, Einstein's vacuum field equations 
for the metric (|l8|) imply the well-known Ernst equation (see e.g. |^) - a 
nonlinear partial differential equation for a complex function f of p and z: 



f,pp ~l~ f,zz + f.p 

P . 



= fl + f% (37) 



(Re/) 

where the Ernst potential 

/(p,z) = e2^ + z6, (38) 

with U{p, z) being the function entering the metric (|3^), function b{p, z) is a 
"potential" for A{p,z) in (|^, 

4p = pe-4%, 4, ^-pe-^X, (39) 



and the last function k{p, z) in ( |37| ) can be determined from U and b by 
quadratures. 

The Ernst equation can be regarded as the integrability condition of a 
system of linear equations for a complex matrix <P, which is a function of 
p + iz, p — iz, and of a (new) complex parameter A. Knowing one can de- 
termine / from ^ at A = 1 . Now the problem of solving the linear system can 
be reformulated as the so called Riemann-Hilbert problem in complex func- 
tion theory. (This, very roughly, means the following: let K he a closed curve 
in the complex plane and F{K) a matrix function given on A'; find a matrix 
function which is analytic inside i, and <Pout analytic outside K such 
that ^in^out = F on K.) The Riemann-Hilbert problem can be formulated 
as an integral equation. The hardest problem with which Neugebauer and 
Meinel were faced was in connecting the specific physical boundary values of 
/ on the disk with the functions entering the Riemann-Hilbert problem (with 
contour K being determined by the position of the disk in the p, z plane) , 
and with the corresponding integral equation. The fact that they succeeded 
and found the solution of their integral equation is a remarkable achievement 
in mathematical physics. The gravitational field and various physical char- 
acteristics of the disk (e.g. the surface density) are given up to quadratures 



in terms of ultraelliptic functions |192|, which can be numerically evaluated 



without difficulties. This result, however, may appear as a "lucky case": it 
does not imply that one will be able to tackle similarly more complicated sit- 
uations as, for example, thin disks with pressure, with non-uniform rotation, 
or 3-dimensional rotating bodies such as neutron stars. 

Many physical characteristics of uniformly rotating relativistic disks such 
as their surprisingly high binding energies, the high redshifts of photons emit- 
ted from the disks, or the dragging of inertial frames in the vicinity of the 



disks, were already obtained with remarkable accuracy in |187 , as the exact 
solution now verifies. Here we only wish to demonstrate the fundamental dif- 
ference between the Newtonian and relativistic case, as it is illustrated in Fig. 
8. The rigidly rotating disk of dust of Neugebauer and Meinel represents the 
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Fig. 8. The general relativistic ("Einsteinian") thin disk of rigidly rotating dust 
constructed by Neugebauer and Meinel, compared with the analogous disk in New- 
tonian theory. If th e an gular momentum is too low, the disk forms a rotating (Kerr) 
black hole. (From ]l93t.) 



relativistic analogue of a classical Maclaurin disk. For the Maclaurin disk, 
it is easy to show that the (dimensionless) quantities y — 2GfiM / and 
X = GM^/cJ (M and J are the total mass and angular momentum respec- 
tively, and f2 is the angular velocity) are related hy y = (97r^/125)a;'^. For a 
fixed M the angular velocity f2 ^ y can be increased arbitrarily, with J being 
correspondingly decreased. For relativistic disks, however, there is an upper 
bound on Q given by fimax = C^/2GM, whereas J is restricted by the lower 
bound Jmin = GAP /c. With an angular momentum too low, a rigidly rotat- 
ing disk cannot exist. If we "prepare" such a disk, it immediately begins to 
collapse and forms - assuming the cosmic censorship - a rotating Kerr black 
hole with x = GM'^/cJ > 1. (Notice that the assumption of rigid rotation 
is here crucial: the differentially rotating disks considered in the preceding 
section can have an arbitrary value of x.) Since one can define the angular 
velocity fl{M, J) of the horizon of a Kerr hole, one may consider y{x) for 
black hole states with x> 1 (cf. Fig 8). The rigidly rotating disk states and 
the black hole states just "meet" at y = a; = 1. In the ultrarelativistic limit 
the gravitational field outside the disk starts to be unaffected by the detailed 
structure of the disk - it approaches the field of an extremely rotating Kerr 
black hole with x = 1. Such a result had already been obtained by Bardeen 
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and Wagoner. However, it is only now, with the exact solution available, that 
it can be investigated with full rigor. It gives indirect evidence that Kerr 
black holes are really formed in the gravitational collapse of rotating bodies. 

As noticed also by Bardcen and Wagoner, in the ultrarelativistic limit 
the disk itself "becomes buried in the horizon of the extreme Kerr metric, 
surrounded by its own infinite, non asymptotically flat universe" (see [194] for 
a recent detailed analysis of the ultrarelativistic limit). Similar phenomena 
arise also in th e cas e of some spherical solutions of Einstein- Yang-Mills-Higgs 
equations (cf. [195| and Section 13). 



7 Taub-NUT space 

The name of this solution of vacuum Einstein's equations fits both to the 
names of its discoverers (Taub-Newman-Unti-Tamburino) a nd to its curious 



properties. Owing to these properties (which induced Misner |196] to consider 
the solution "as a counterexample to almost anything"), this spacetime has 
played a significant role in exhibiting the type of effects that can arise in 
strong gravitational fields. 



Taub 1 197] discovered an empty universe with four global Killing vectors 
almost half a century ago, during his pioneering study of metrics with several 
symmetries. By continuing the Taub universe through its horizon one arrives 
in NUT space. NUT space, however, was only discovered in 1963 by a different 



method ]198]. In fact, it could have been obtained earlier by applying the 



transformation given in Jiirgen Ehlers' dissertation ]54 and his talk at the 
GR2 conference in Royaumont in 1959 ]Q. This transformation gives the 
recipe for obtaining stationary solutions from static ones. How the NUT 
space can be obtained by applying this transformation to the Schwarzschild 
metric was demonstrated explicitly by Ehlers at GR4 in London in 1965 p&. 



7.1 A nevf way to the NUT metric 

Here we shall briefly mention a simple, physically appealing new derivation 
of the NUT metric given recently by Lyndcn-BcU and Nouri-Zonoz (LBNZ) 
[ [199| . Their work also shows how even uncomplicated solutions may still be 
of interest in unexpected contexts. LBNZ's inspiration to study the NUT 
space has in fact come from Newton's Principia! In one of his scholia Newton 
discusses motion under the standard central force plus a force which is normal 
to the surface swept out by the radius vector to the body which is describing 
the non-coplanar path. A simple interesting case is the motion of mass tuq 
satisfying the equation 

™o 5 = -^'(^) f + ^ ^ X -Bg, (40) 

where f = r/r, 

Bg = -Q f/r2, Q = Q c/mo, (41) 
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Q and Q are constants, and c is the velocity of light. Here we write c explicitly 
though c = 1, to make the analogy with magnetism. Indeed, Bg is the field of 
a "gravomagnetic" monopole of strength Q. The classical orbits of particles 



lie on cones which, if the monopole is absent, flatten into a plane [199|. 

It was known that NUT space corresponds to the mass with a gravomag- 
netic monopole, but this was never used in such a physical way as by LBNZ 
for its derivation. The main point is to start from the well-known split of the 



stationary metrics as described in Landau and Lifshitz [139| (see [200| for a 
covariant approach, and the contribution of Beig and Schmidt in the present 
volume) 

ds^ = -e-^^idt - A,dx'f + -iijdxUx\ (42) 

where v^Ai,^ij are independent of t. This form is unique up to the choice 
of time zero: t = t + x(x*) implies again the metric in the form ( ^2| ) in 
{t , a;'), with the "vector potential" undergoing a gauge transformation A[ = 
■^i + ViX- Writing down the equation of motion of a test particle in metric 
(^2[), in analogy with the equation of motion of a charged particle in an 
electromagnetic field, one is naturally led to define the "gravoelectric" and 
"gravomagnetic" fields by 

Eg ^Vv, Bg=V X A, (43) 



where "Vx" is with respect to 7ij. Following the problem of §95 in [139| 
one then rewrites all Einstein's equations in terms of the fields (^), the 
metric 7ij, and their derivatives. To find the vacuum spherically symmetric 
spatial 7-metric one takes "fijdx^dx^ = e^^dr^ +r'^{d9'^ -|-sin^ 6 dip"^), and one 
assumes v = i^(r), and Bg — —Qe^^/r^. The Einstein equations then imply 
the spacetime metric, which is not spherically symmetric, in the form 

ds^ = -e-'^''(^dt-2q{l + C0Be)d<^ + (l - / r"^)'^ e^" dr'^ 

where q — Q /2 ^ constant, 

= 1 - 2r-2 [q^ + mv/r2 - , (45) 
and the vector potential = 2q{\ + cos 9) satisfies (p3|). (The factor (1 — 



q /r ) should be raised to the power —1 in equation (3.22) in [199|, as it 
is clear from (3.20).) Equation (Q) is the NUT metric, with r being the 
curvature coordinate of spheres r = constant. With g = the metric ( p[ ) 
becomes the Schwarzschild metric in the standard Schwarzschild coordinates. 
More commonly the metric (E4) is written in the form 



ds"^ = -V[di+Aq sin^ - d(p ) + V "^dr^ -f (f ^ + q^) [dO^ + sin^ dtp^ 



(46) 



64 Jin Bicak 



which can be obtained from (Bj) by putting 



v/r2 - q^, i=t- Aqip. (48) 



Recently, Ehlers [201| considered a Newtonian hmit of NUT space within 
his frame theory which encompasses general relativity and Newton-Cartan 
theory (a slight generalization of Newton's theory). The main purpose of 
Ehlers' frame theory is to define rigorously what is meant by the statement 
that a one parameter family of relativistic spacetime models converges to a 
Newton-Cartan model or, in particular, to a strictly Newtonian model. 

The strictly Newtonian limit occurs when the Coriolis angular velocity 
field w, related to the connection coefficients F^^ in the Newton-Cartan the- 
ory, depends on time only. NUT spacetimes approach a truly Newton-Cartan 
limit with spatially non-constant radial Coriolis field w'" ~ —q/r'^, which in 
this limit coincides with the Newtonian gravomagnetic field. As in the anal- 
ogous classical problem with the equation of motion (^0[), the geodesies in 
NUT space lie on cones. This result has been used to study gravitational lens- 
ing by gravomagnetic monopoles ]199| | : they twist the rays that pass them in 
a characteristic manner, different from that due to rotating objects. 

The metrics and (^ ) appear to have a preferred axis of fixed points 
of symmetry. This is a false impression since we can switch the axis into 
any direction by a gauge transformation. For example, the metric (|4^ ) has 
a conical singularity at = but is regular at = tt, whereas the metric 
( ^6|) has a conical singularity at 6* = tt but is regular at = 0. The metrics 
are connected by the simple gauge transformation, i.e. t ^ t — t — Aqip. 
A mass endowed with a gravomagnetic monopole appears as a spherically 
symmteric object but the spacetime is not spherically symmetric according 
to the definition given in Section 2.1. Nevertheless, there exist equivalent 
coordinate systems in which the axis can be made to point in any direction 
- just as the axis of the vector potential of a magnetic monopole can be 
chosen arbitrarily. For further references on interpreting the NUT metric as 
a gravomagnetic monopole, see |199{ and the review by Bonnor [p3. 



7.2 Taub-NUT pathologies and applications 

By introducing two coordinate patches, namely the coordinates of metric ( ^4[ ) 
to cover the south pole {9 = tt) and those of ( p^ ) the north pole (9 = 0), the 
rotation axis can be made regular. However since ip is identified with period 
27r, equation ( ^ ) implies that t and i have to be identified with the period 
87rg. Then observers with (f, 9, ip) = constant follow closed timelike lines if V 
in ( p7| ) is positive, i.e. if f > fp = m + (m^ + q^) 2 . The hypersurface f = fo is 
the null hypersurface - horizon, below which lines {t, 9, ip) — constant become 
spacelike. Because of the periodic identification of t and t, the hypersurfaces 
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of constant f cliange the topology from 5*^ x to S^, on which t/2q, 6, ip 
become Euler angle coordinates. 

The region with y < is the Taub universe: it has homogeneous but 
non-isotropic space sections f = constant. The coordinate f , allowed to run 
from — oo to +oo, is a timelike coordinate, and is naturally denoted by t in 
the Taub region. 

In addition to the closed timelike lines in the NUT region there are fur- 
ther intriguing pathologies exhibited by the Taub-NUT solutions . He re we 
just list some of them and refer to the relevant literature |^ , ^ ,196|. The 
Taub region is globally hyperbolic: its entire future and past history can 
be determined from conditions given on a spacelike Cauchy hypersurface. 
However, this is not the case with the whole Taub-NUT spacetime. As in 
the Reissner-Nordstrom spacetimes (Section 3), there are Cauchy horizons 
H^[E) of a particular spacelike section S of maximal proper volume ly- 
ing between the globally hyperbolic Taub regions and the causality violating 
NUT regions. H^{S) are smooth, compact null hypersurfaces diffeomorphic 
to - the generators of such null surfaces are closed null geodesies. The 
Taub region is limited between <_ < ^ < where t± are roots of V in 
equation ( ^7| ) (with the interchange i <~> f) . This region is compact but there 
are timelike and null geodesies which remain within it and are not complete. 
(See [ p7| for a nice picture of these geodesies spiralling around and approach- 
ing H^{S) asymptotically.) This pathological behaviour of "the incomplete 
geodesies imprisoned in a compact neighbourhood of the horizon" was inspi- 
rational in the definition of singularities flT^ ~ one meets here the example in 
which the geodesic incompleteness is not necessarily connected with strong 
gravitational fields. It can be shown, however, that after an addition of even 
the slightest amount of matter this pathological behaviour will not take place 
- true singularities arise. 

This enables one to consider the time between t_ and as the lifetime 
of the Taub universe. Wheeler |202| constructed a specific case of the Taub 
universe which will live as long as a typical Friedmann closed dust model 
(~ 10^*^ years) but will have a volume at maximum expansion smaller by 
a factor of 5 x 10^°. This example thus appears to be a difficulty for the 
anthropic principle. 

Taub space seems also to be the only known example giving the possibility 
of making inequivalent NUT-like extensions which lead to a non-Hausdorff 
spacetime manifold ^6 200 202 1. 

The Taub-NUT solution plays an important role in cosmology and quan- 
tum gravity. Here we wish to note yet two other recent applications of this 
space. About ten years ago, interest was revived in closed timelike lines, time 
machines, and wormholes. One of the leaders in this activity. Kip Thorne, 
explains in 1 204 1 in a pedagogical way the main recent results on closed time- 
like curves and wormholes by using "Misner space" - Minkowski spacetime 
with identification under a boost, which Misner introduced as a simplified 
version of Taub-NUT space. 
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The second application of the Taub-NUT space is still more remarkable - 
it plays an important role outside general relativity. The asymptotic motion of 
monopoles in (super-)Yang-Mills theories corresponds to the geodesic motion 



in Euclidean Taub-NUT space [205]. Euclidean Taub-NUT spaces have been 



discussed in many further works on monopoles in gauge theories. One of 



the latest of these works |206|, on the exact T-duality (which relates string 
theories compactified on large and small tori) between Taub-NUT spaces and 
so called "calorons" (instantons at finite temperature defined on x S^), 
gives also references to previous contributions. 

8 Plane waves and their collisions 
8.1 Plane-fronted waves 

The history of gravitational plane waves had began already by 1923 with 
the paper on spaces conformal to flat space by Brinkmann. Interest in these 
waves was revived in 1937 by Rosen, and in the late 1950s by Bondi, Pirani 
and Robinson, Holy, and Peres (see [ p3|j6l| ] for references). A comprehensive 
geometrical approach to these spacetimes soon followed in the classical trea- 
tise by Jordan, Ehlers and Kundt jSTj , and in the subsequent well-known 
chapter by Ehlers and Kundt [53| . As an application of various newly devel- 
oped methods to analyze gravitational radiation, and as a simple background 
to test various physical theories, plane waves have proved to be a useful and 
stimulating arena which offers interesting contests even today, as we shall 
indicate by a few examples in Section 8.2. 

Consider a congruence of null geodesies (rays) x°'{v) such that dx'^ /dv = 
k", kak" — 0, ka-fjk^ — 0, V being an affine parameter. In general a geodesic 
congruence is characterized by its expansion 9, shear \a\ and twist uj given 
by (see e.g. @) 



lk%^ Wl^Jlho^-.mk'^'^-d', (49) 



(50) 

According to the definition given by Ehlers and Kundt a vacuum 
spacetime is a "plane-fronted gravitational wave" if it contains a shearfree 
\a\ — geodesic null congruence, and if it admits "plane wave surfaces" 
(spacelike 2-surfaces orthogonal to /c"). This definition is inspired by plane 
electromagnetic waves in Maxwell's theory. Electromagnetic plane waves are 
null fields ("pure radiation fields"): there exists a null vector k", tangent to 
the rays, which is transverse to the electromagnetic field -Fk/j, i.e. F^pk^ = 0, 
F*^k^ = 0, and the quadratic invariants of which vanish, FapF"^ = — 
Fa/sF*"^, where F* p is dual to Fa/s- Analogously, Petrov type N gravita- 
tional fields (see pn) are null fields with rays tangent to fc" (the "quadru- 
ple Debever-Penrose null vector"), and with the Riemann tensor satisfying 
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Ra^^sk^ = 0, Rafj^sR"'^'^^ = 0, and R"''^^^ R^p^^s = 0.0 Then the Bianchi 
identities and the Kundt- Thompson theorem for type N sohitions in vacuum 
spacetimes (also more generally, under the presence of a nonvanishing cos- 
mological constant) imply that the shear of k" must necessarily vanish (see 
]). Because of the existence of plane wave surfaces, the expansion ( p9| ) 
and twist ( |50| ) must vanish as well, = lo = 0. In this way we arrive at 
the Kundt class of nonexpanding, shearfree and twistfree gravitational waves 
[ |6T| . The best known subclass of these waves are "plane-fronted gravitational 
waves with parallel rays" (pp-waves) which are defined by the condition that 
the null vector k" is covariantly constant, ka-f} = 0. Thus, automatically fc" 
is the Killing vector, and 9 ^ \a\ ^ oj ~ 0. 

Ehlers and Kundt give several equivalent characterizations of the 
pp-waves and show, following their previous work \ ^7\ , that in suitable null 
coordinates with a null coordinate u such that fc^ — u^a and fc" — (d/dv)", 
the metric has the form 

ds'^ = 2dCdC - 2dudv - 2H{u, C, C)du^, (51) 

where H is a real function dependent on u, and on the complex coordinate C 
which spans the wave 2-surfaces u = constant, v = constant. These 2-surfaces 
with Euclidean geometry are thus contained in the wave hypersurfaces u = 
constant and cut the rays given by (w, Q = constant, v changing. The vacuum 
field equations imply 2-dimensional Laplace's equation 

H,cc = 0, (52) 

so that we can write 

2H^fiu,C) + f{u,C), (53) 

where /(m, ^) is an arbitrary function of u, analytic in ^. To characterize the 
curvature in the waves and their effect on test particles it is convenient to 
introduce the null complex tetrad, such that at each spacetime point, together 
with the preferred null vector fc", we have a null vector l^J^ka = —1, and 
complex spacelike vector m" satisfying niafh"' = l,mak" = nial" = 0. For 
the metric (|5l| ) the only nonvanishing projection of the Weyl (in the vacuum 
case, the Riemann) tensor onto this tetrad is the (Newman-Penrose) scalar 

which denotes a transverse component of the wave propagating in the fc" 
direction. As shown by Ehlers and Kundt [|3| (see also e.g. ||6^), though in 
a somewhat different notation, we can use again an analogy with the elec- 
tromagnetic field - described for an analogous plane wave by the transverse 

This algebraic (local) analogy between null fields exists also between electromag- 
netic and gravitational shocks (possible discontinuities across null hypersurfaces) , 
and in the asymptotic behaviour of fields at large distances from sources (the 
"peeling property" - see e.g. p£|p7| ). 
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component 02 = Fapfhfl^ - and write !f4 = Ae'^, where real A > is con- 
sidered as the amplitude of the wave, and at each spacetime point associate 
with the plane of polarization. Vacuum pp-waves with — constant are 
called linearly polarized. 

Consider a free test particle (observer) with 4- velocity u and a neighbour- 
ing free test particle displaced by a "connecting" vector Z°'(t). Introducing 
then the physical frame e(j-| which is connected with the observer such that 
6(0) ~ u and e(j) are connected with the null tetrad vectors by 

™= 75 (*'(i)+^^(2))' m - ^ (e(i) - ie(2)) , 
we find that the equation of geodesic deviation in spacetime with only S'4 7^ 



implies (see |207l) 



- + ^ + AxZ(i), = 0, (56) 

where A+ = ^ Re tf4, Ax = ^ Im are amplitudes of and "x" polar- 
ization modes, and Z^'^ are the frame components of the connecting vector 
Z. Since the frame vector 6(3) is chosen in the longitudinal direction (the di- 
rection of the rays), equation ( p6[ ) clearly exhibits the transverse character of 
the wave. If particles, initially at rest, lie in the (e(i), 6(2)) plane, there is no 
motion in the longitudinal direction of 6(3). The ring of particles is deformed 
into an ellipse, the axes of different polarizations are shifted one with respect 
to the other by ^ (such behaviour is typical for linearized gravitational waves 

- cf. e.g. [Q). Making a rotation in the transverse plane by an angle 

e'(i) = cosz?e(i) + sint9e(2) , e'(2) = — sini^eji) + cosi?e(2) , (57) 

and taking d = d+ir) = -\ Arg = -\0, then A'^ = A'^ = 

- the wave is purely polarized, li = constant, the rotation angle is 
independent of time - the wave is rightly considered as linearly polarized. 

Hence, with the discovery of pp-waves, the understanding of the proper- 
ties of gravitational radiation has become deeper and closer to physics. In 
addition, the pp-waves can easily be "linearized" by taking the function H 
in the metric (|5l|) to be so small that the spacetime can be considered as a 
perturbation of Minkowski space within the linearized theory. Such an "easy 
way" from the linear to fully nonlinear spacetimes is of course paid by their 
simplicity. 

In general, in fact, the pp-waves have only the single isometry generated by 
the Killing vector k" = (d/dv)" . However, a much larger group of symmetries 
may exist for various particular choices of the function H(u,(^,(^). Jordan, 
Ehlers and Kundt Q (see also [ p3pl| ]) gave a complete classification of the 
pp-waves in terms of their symmetries and corresponding special forms of H. 
For example, in the best known case of plane waves to which we shall turn 
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in greater detail below, W/^ is independent of ^, so that after removing linear 
terms in C by a coordinate transformation, we have 

H{u,C,C)^A{u)e+A{u)(\ (58) 

with A(u) being an arbitrary function of u. This spacetime admits five Killing 
vectors. 



Recently, Aichelburg and Balasin [P08 20£ ] generalized the classification 



given in 1 57 by admitting distribution- valued profile functions and allowing 
for non- vacuum spacetimes with metric ([si]), but with H which in general 
does not satisfy (|5^). They have shown that with H in the form of delta- like 
pulses, 

HiuXX)^ fi(:X)6iu), (59) 

new symmetry classes arise even in the vacuum case. 

The main motivation to consider impulsive pp-waves stems from the met- 
rics describing a black hole or a "particle" boosted to the speed of light. The 



simplest metric of this type, given by Aichelburg and Sexl [210|, is a Schwarz- 
schild black hole with mass m boosted in such a way that fi = mj \J\ — 
is held constant as ui ^ 1. It reads 

ds^ = 2dC,dl - 2dudv ~ 4/ilog(CC)'5(M)d^i^ (60) 

with H clearly in the form (|5^). This is not a vacuum metric: the energy- 
momentum tensor Ta/s = iiS{u)6{Qkakp indicates that there is a "point-like 
particle" moving with the speed of light along u = 0. The Aichelburg-Sexl 
metric and its more recent generalizations have found interesting applications 
even outside of general relativity. Some of them will be briefly mentioned in 
Section 8.2. 

Let us now turn to the simplest class of pp-waves, which comprises of the 
best known and illuminating examples of exact gravitational waves. These are 
the plane waves. They are defined as homogeneous pp-waves in the sense that 
the curvature component (se e (|54| )) is constant along the wave surfaces so 
that function H is in the form (|58|). One can write H as in ( p3| ) where 

/(u,C)-i^(^)e'^(")C', (61) 

with linear terms being removed by a coordinate transformation. Just as a 
plane electromagnetic wave, a plane gravitational wave is thus completely 
represented by its amplitude A{u) and polarization angle 0{u) as functions 
of the phase u. 

The plane waves, including their generalization into the Einstein-Maxwell 
theory (an additional term B{u)(^^ then appears in H, both !f4 and the elec- 
tromagnetic quantity $2 being independent of C), were already studied in 
1926 (see |^). A real understanding however came only in the late 1950s. 
Ehlers and Kundt j5^ give various characterizations of this class. For ex- 
ample, they prove that a non-flat vacuum fleld is a pp-wave if and only if 
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the curvature tensor is complex recurrent, i.e. if PaP'yS,ii = Pap-ySQ^n where 
PaP'yS = RaP'yS + * Rap-yS', and it is a plane wave if and only if the re- 
currence vector is collinear with a real null vector. They also state a 
nice theorem showing that the plane wave spacetimes defined by the met- 
ric H and / given by (p3[), (|69|), C, = x + iy, and with coordinate 
ranges — oo < x,y,u,v < oo, are geodesically complete if functions A{u) and 
0{u) are C^-functions. Quoting directly from | |53[ |, "there exist ... complete 
solutions free of sources (singularities), proving to think of a graviton field 
independent of any matter by which it be generated. This corresponds to the 
existence of source- free photon fields in electrodynamics". Ehlers and Kundt 
also state an open problem which, as far as I am aware, has not yet 
been solved: to prove that plane waves are the only geodesically complete 
pp- waves. 

The most telling examples of plane waves are sandwich waves. The ampli- 
tude A{u) in ( |6l|) need not be smooth: either it can only be continuous and 
nonvanishing on a finite interval of u (sandwich), or a step function (shock), 
or a delta function (impulse) . A physical interpretation of such waves is bet- 
ter achieved in other coordinate systems, in which the metric "before" and 
"after" the wave is not Minkowskian but has a higher degree of smooth- 
ness. For linearly polarized waves (6* equal to zero), a convenient coordinate 
system can be introduced by setting (see e.g. [ ^ll[ ) ( ~ (l/\/2)(pa; -I- i^y), 
V = (l/2)(t -f z + pp' x^ qq'y^), u = t — z, where ' = d/du, and functions 
p — p{u) and q — q{u) solve equations p" + A{u)p = and q" — A{u)q = 0. 
In these coordinates the metric turns out to be 

ds^ = -dt^ + p^dx^ + q^dy^ + dz^ . (62) 

In double-null coordinates w, v, with u = u^t — z, v — t + z, and with a 
general polarization, the metric can be cast into the form (see e.g. |65| , |212|] ) 

ds^ ^ -dudd + (e^ coahWdx^ + e-^ coshWdy^ -2amhWdxdy), (63) 

where U, V, W depend on u only. This so called Rosen form was used in t he 
classical paper on exact plane waves by Bondi, Pirani and Robinson ||213| . 

A simple, textbook example [pl4| of a sandwich wave is the wave with a 
"square profile": A{u) = for u < and u > a'^,A{u) = = constant for 



< u < a^. The functions p and q which enter (62) are then p = q = 1 at 
u < 0,p = coa{u/a),q — cosh(M/a) at < u < a^, and p = —(u/a)sma + 
constant, g — (w/a)sinha + constant at < u. This example can be used 
to demonstrate explicitly various typical features of plane sandwich gravita- 
tional waves within the exact theory: (i) the wave fronts travel with the speed 
of light; (ii) the discontinuities of the second derivatives of the metric tensor 
are permitted along a null hypersurface, but must have a special structure; 
(iii) the waves have a transverse character and produce relative accelera- 
tions in test particles; (iv) the waves focus astigmatically initially parallel 
null congruences (rays) that are pointing in other directions than the waves 
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themselves; (v) as a consequence of the focusing, Rosen-type line elements 
contain coordinate singularities on a hypersurface behind the waves, and in 
general caustics will develop there |214]. 

The focusing effects imply a remarkable property of plane wave space- 
times: no spacelike global hypersurface exists on which initial data can be 
specified, i.e. plane wave spacetimes contain no global Cauchy hypersurface. 
This can be understood from Fig. 9. Considering a point Q in flat space in 
front of the wave, Penrose | |215{ has shown that its future null cone is dis- 
torted as it passes through the wave in such a manner that it is refocused to 
either a point i? or a line passing through R parallel to the wave front. Any 
possible Cauchy hypersurface going through Q must lie below the future null 
cone through Q, i.e. below the past null cone of R. Hence, it cannot extend 
as a spacelike hypersurface to spatial infinity. 




8.2 Plane-fronted waves: new developments and applications 



The interest in impulsive waves generated by boosting a "particle" at rest to 
the velocity of light by means of an appropriate limiting procedure persists 
up to the present. The ultrarelativistic limits of Kerr and Kerr-Newman black 
holes were obtained in | 216| ,217, 218|, and recently, boosted static multipole 
(Weyl) particles were studied] [219| . Impulsive gravitational waves were also 
generated by boosting the Schwarzschild-d e Sitter a nd Schwarzschild-anti de 
Sitter metrics to the ultrarelativistic limit 220 ,221]. 

These types of spacetimes, especially the simple Aichelburg-Sexl metrics, 
have been employed in current problems of the generation of gravitational 
radiation from axisymmetric black hole collisions and black hole encounters. 
The recent monograph by d'Eath [ p22[ gives a comprehensive survey, includ- 
ing the author's new results. There is good reason to believe that spacetime 
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metrics produced in high speed coUisions wiU be simpler than those corre- 
sponding to (more reahstic) situations in which black holes start to collide 
with low relative velocities. The spacetimes corresponding to the collisions at 
exactly the speed of light is an interesting limit which can be treated most 
easily. Aichelburg-Sexl metrics are used to describe limiting "incoming states" 
of two black holes, moving one against the other with the speed of light. An 
approximation method has been developed in which a large Lorentz boost is 
applied so that one has a weak shock propagating on a strong shock. One 
finds an estimate of 16.8 % for the ef ficien cy of gravitational wave generation 
in a head-on speed-of- light collision |222|. 

Great interest has been stimulated by 't Hooft's 223 work on the quan- 
tum scattering of two pointlike particles at centre-of-mass energies higher or 
equal to the Planck energy. This quantum process has been shown to have 
close con nection with classical black hole collisions at the speed of light (see 
[222 224] and references therein). 

Recently, the Colombeau algebra of generalized functions, which enables 
one to deal with singular products of distributions, has been brought to gen- 
eral relativity and used in the description of impulsive pp-waves in various 
coordinate systems |225|, and also for a rigorous solution of the geodesic and 
geodesic deviation equations for impulsive waves [226|. The investigation of 
the equations of geodesies in non- homogeneous pp-waves (with f C^) has 
shown that the motion of test particles in these spacetimes is described by 
the Henon-Heiles Hamiltonian which implies that the motion is chaotic |227|. 

Plane-fronted waves have been used as simple metrics in various other 
contexts, for example, in quantum field theory on a given background (see 
[ ^28| for recent work), and in string theory [229|. As emphasized very recently 
by Gibbons [230|, since for pp-waves and type Kundt's class (see the begin- 
ning of Section 8.1) all possible invariants formed from the Weyl tensor and 
its covariant derivatives vanish | p3l| , these metrics suffer no quantum correc- 
tions to all loop orders. Thus they may offer insights into the behaviour of a 
full quantum theory. The invariants vanish also in type /// spacetimes with 
nonexpanding and nontwisting rays [232|. 



8.3 Colliding plane waves 

As with a number of other issues in gravitational (radiation) theory, the 
pioneering ideas on colliding plane gravitational waves are connected with 
Roger Penrose. It does not seem to be generally recognized that the basic 
idea appeared six years before the well-known paper by Khan and Penrose 
[ ^33| in which the metric describing the general spacetime representing a 
collision of two parallel-polarized impulsive gravitational waves was obtained. 
Having demonstrated the surprising fact that general relativistic plane wave 
spacetimes admit no Cauchy hypcrsurface due to the focusing effect the waves 
exert on null cones, Penrose ]215[ (in footnote 12) remarks: "This fact has 
relevance to the question of two colliding weak plane sandwich waves. Each 
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wave warps the other until singularities in the wave fronts ultimately appear. 
This, in fact, causes the spacetime to acquire genuine physical singularities in 
this case. The warping also produces a scattering of each wave after collision 
so that they cease to be sandwich waves when they separate (and they are 
no longer plane - although they have a two-parameter symmetry group)." 
The first detailed study of colliding plane waves, independently of Khan 



and Penrose, was also undertaken by Szekeres (see [234,235|). He formulated 
the problem as a characteristic initial value problem for a system of hyper- 
bolic equations in two variables (null coordinates) u, v with data specified on 
the pair of null hypersurfaces u = 0,v = intersecting in a spacelike 2-surface 
(Fig. 10). In the particular case of spacetimes representing plane waves prop- 
agating before the collision in a flat background, Szekeres has shown that 
coordinates (of the "Rosen type" , as known from the case of one wave - see 
Eq. (|63|)) exist in which the metric reads 

ds'^ = - e^^'du dv + 

+ e-^ [ cosh Wdx^ + 6'^ cosh Wdy^ - 2 sinh Wdx dy] , (64) 

where M , U, V and W are functions of u and v. Coordinates x and y are 
aligned along the two commuting Killing vectors d/dx and d/dy, which are 
assumed to exist in the whole spacetime representing the colliding waves (cf. 
the note by Penrose above). In almost all recent work on colliding waves, 
region IV in Fig. 10, where u < 0, w < 0, is assumed to be fiat. The null lines 
u = 0,v < and v — 0,u < are wavefronts, and in regions II {u < 0,v > 0) 
and /// (u > 0, w < 0) one has the standard plane wave metric corresponding 
to two approaching plane waves from opposite directions. In region IT, func- 
tions M, U, V, W depend on v only, and in region III only on u. The waves 
collide at the 2-surface u = u = 0, in region / they interact. The spacetime 
here can be determined by the initial data posed on the v > portion of the 
hypersurface u = (which in Fig. 10 are "supplied" by the wave propagating 
to the right) and by the data on the u > portion of the hypersurface v — 
(given by the wave propagating to the left). Unfortunately, the integration of 
such an initial value problem does not seem to be possible for general incom- 
ing wave forms and polarizations. If, however, the approaching waves have 
constant and aligned (parallel) polarizations, one may set the function W — 
globally. The solution of the initial value problem then reduces to a one di- 
mensional integral for the function V, and two quadratures for the function 
M. (The function exp(— C/) must have the form f{u) + g{v) as a consequence 
of the field equations everywhere, and it can be determined easily from the 
initial data.) Despite these simplifications it is very difficult to obtain exact 



solutions in closed analytic form. Szekeres |235| found a solution (as he puts it 



"more or less by trial and error") which, as special cases, includes the solution 



given by himself earlier |234| and the solution obtained independently and 



simultaneously by Khan and Penrose [233|. Although Szekeres' formulation 



of a general solution for the problem of colliding parallel-polarized waves is 



74 Jiff Bicak 



difficult to use for constructing other specific explicit examples, it has been 
employed in a general analysis of the structure of the singularities produced 
by the collision p36[, which will be discussed in the following. 




Fig. 10. The spacetime diagram indicating the collision of two plane-fronted grav- 
itational waves which come from regions //and ///, collide in region /, and produce 
a spacelike singularity. Region IV is flat. 



It has also inspired an important, difficult piece of mathematical physics 
which was develop ed a t the beginning of the 1990s in the series of papers by 
Hauser and Ernst |237]. Their new method of analyzing the initial value prob- 
lem can be used also for the case when the polarization of the approaching 
waves is not aligned. They formulated the initial value problem in terms of 
the equivalent matrix Riemann-Hilbert problem in the complex plane. Their 
techniques are related to those used by Neugebauer and Meinel to analyze 
and construct the rotating disk solution as a boundary value problem (Sec- 
tion 6.3). No analogous solution for colliding waves in the noncollinear case is 
available at present, but investigations in this direction are still in progress. 
Most recently, Hauser and Ernst prepared an extensive treatise [238| in which 
they give a general description and detailed mathematical proofs of their 
study of the solutions of the hyperbolic Ernst equation. 

The approach of Khan and Penrose for obtaining exact solutions describ- 
ing colliding plane waves starts in the region / where the waves interact: (i) 
find a solution with two commuting spacelike Killing vectors d/dx and d/dy, 
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transform to null coordinates, and look back in time whether this solution 
can be extended across the null hypersurface u = 0, w = so that it describes 
a plane wave propagating in the u-direction in region // and another plane 
wave propagating in the w-direction in region /// ; (ii) satisfy boundary con- 
ditions not only across boundaries between regions / and //, and regions / 
and ///, but also across the boundaries between //and IV, and ///and IV in 
such a manner that IV is flat. The original prescription of Khan and Penrose 
for extending the solution from region / to regions // and /// consists in the 
substitutions uH{u) and vll{v) in place of u and v everywhere in the metric 
coefficients; here //(m) = 1 for w > 0, // = for u < is the usual Heaviside 
function. We then get the metric as a function of v (respectively u) in region 
// (respectively III) corresponding to the wave propagating to the right (re- 
spectively to the left) in Fig. 10. Finally, it remains to investigate carefully 
the structure of discontinuities and possible singularities on the null bound- 
aries between these regions. In the original Khan and Penrose solutions the 
Riemann tensor has a (5-function character on the boundaries between // and 
IV, and /// and IV; but inside regions // and /// themselves the spacetime 
is flat (the collision of impulsive plane waves). In the solution obtained by 
Szekeres [[235[ , regions //and ///are not flat, and the Riemann tensor at the 
boundaries between // (respectively /// ) and IV is just discontinuous (the 
collision of shock waves). 

Nutku and Halil |239 constructed an exact solution describing the colli- 
sion of two impulsive plane waves with non-aligned polarizations. In the limit 
of collinear polarizations their solution reduces to the solution of Khan and 
Penrose. All of these solutions reveal that the spacelike singularity always 
develops in region / (given by + = 1 in Fig. 10) ~ in agreement with the 
original suggestion of Penrose. Moreover, the singularity "propagates back- 
ward" and s o ca lled fold singularities, analyzed in detail in 1984 by Matzner 
and Tipler [240|, appear also at u = 1 and u = 1 in regions // and ///. 
This new type of singularity provides evidence of how even relatively recent 
studies of explicit exact solutions may reveal unexpected global features of 
relativistic spacetimes. 

The remarkable growth of interest in colliding plane waves owes much 
to the systematic (and symptomatic) effort of S. Chandrasekhar who, since 
1984, together with V. Ferrari, and with B. Xanthopoulos, pub lished a num- 
ber of papers on colliding plane vacuum gravitational waves |241, 242| , and 
on gravitational waves coupled with electromagnetic waves, with null dust, 
and with perfect fluid (see ]212| | for references). The basic strategy of their 
approach follows that of Khan and Penrose: first a careful analysis of the pos- 
sible solution is done in the interaction region /, and then one works backward 
in time, extending the solutions to regions //, /// and IV. 

The main new input consists in carrying over the techniques known from 
stationary, axisymmetric spacetimes with one timelike and one spacelike 
Killing vector to the case of two spacelike Killing vectors, d/dx,d/dy, and 
exploring new features. 
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Taking a simple linear solution of the Ernst equation, E = Prj + iQfi, 
where P and Q are real constants which satisfy + = 1, and 77, fj, are 
suitable time and space coordinates, Chandrasekhar and Ferrari show 
that one arrives at the Nutku-Halil solution. In particular, if Q = 0, the Khan- 
Penrose solution emerges. Since by starting from the same simplest form of 
the Ernst function in the axisymmetric stationary case one arrives at the Kerr 
solution (or at the Schwarzschild solution for the real Ernst function), we may 
conclude that in region / the solutions of Khan and Penrose and of Nutku 
and Halil are, for spacetimes with two spacelike Killing vectors, the analogues 
of the Schwarzschild and Kerr solutions. This mathematical analogy can be 
generalized to colliding electromagnetic and gravitational wa ves within the 
Einstein-Maxwell theory - Chandrasekhar and Xanthopoulos [ ^43[ found the 
analogue of the charged Kerr-Newman solution. Such a generalization is also 
of interest from a conceptual viewpoint: the (5-function singularity in the Wcyl 
tensor of an impulsive gravitational wave might imply a similar singularity in 
the Maxwell stress tensor, which would seem to suggest that the field itself 
would contain "square roots of the J-function" . 

In the most important paper [242| of the series, Chandrasekhar and Xan- 
thopoulos, starting from the simplest linear solution for the Ernst conjugate 
function _E+ = Prj + iQfi, P^ + = 1, obtained a new exact solution for col- 
liding plane impulsive gravitational waves accompanied by shock waves. This 
solution results in the development of a nonsingular Killing-Cauchy horizon 
instead of a spacelike curvature singularity. The metric can be analytically 
extended across this horizon to produce a maximal spacetime which con- 
tains timelike singularities. (The spacelikc singularity in region / in Fig. 10 is 
changed into the horizon, to the future of which timelike singularities occur.) 
In the region of interaction of the colliding waves, the spacetime is isometric 
to the spacetime in a region interior to the ergosphere. 

Many new interesting solutions were discovered by using the Khan and 
Penrose approach. In addition, inverse scattering (soliton) methods and other 
tools from the s olution g eneration techniques were applied. They are reviewed 
in detail in (6|,|l|,|4|. 

Although very attractive mathematical methods are contained in these 
works, one feels that physical interpretation has receded into the background 
- as seemed to be the case when the new solution generating techniques were 
exploited in all possible directions for stationary axisymmetric spacetimes. 
It is therefore encouraging that a more physical and original approach to 
the problem has been initiated by Yurtsever. In a couple of papers he dis- 
cusses Killing-Cauchy horizons |245|] an d the structure of the singularities 
produced by colliding plane waves |236]. Similar to the Cauchy horizons in 
black hole physics, one finds that the Killing-Cauchy horizons are unstable. 
We thus expect that the horizon will be converted to a spacelike singularity. 
By using the approach of Szekeres described at the beginning of this section, 
it is possible to relate the asymptotic form of the metric near the singu- 
larity - which approaches an inhomogeneous Kasner solution (see Section 
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12.1) - to the initial data given along the wavefronts of the incoming waves. 
For specific choices of initial data the singularity degenerates into a coordi- 
nate singularity and a Killing-Cauchy horizon arises. However, Yurtsever's 
analysis |236| shows that such horizons are unstable (within full nonlinear 
theory) against small but generic perturbations of the initial data. These re- 
sults are stronger than those on the instability of the inner horizons of the 
Reissner-Nordstrom or Kerr black holes. In particular, Yurtsever constructs 
an interesting (though unstable) solution which, when analytically extended 
across its Killing-Cauchy horizon, represents a Schwarzschild black hole cre- 
ated out of the collis ion between two plane sandwich waves propagating in a 
cylindrical universe |^36| . 

Yurtsever also introduced "almost plane wave spacetimes" and analyzed 
collisions of almost plane waves [246|. These waves have a finite but very 
large transverse sizes. Some general results can be proved (for example, that 
almost plane waves cannot have a sandwich character, but always leave tails 
behind), and an order-of-magnitude analysis can be used in the discussion of 
the outcome of the collision of two almost plane waves; i.e. whether they will 
focus to a finite minimum size and then disperse, or whether a black hole 
will be created. Although in the case of almost plane waves one can hardly 
hope to find an exact spacetime in an explicit form, this is a field which was 
inspired by exact explicit solutions, and may play a significant role in other 
parts of general relativity. 



9 Cylindrical waves 

In 1913, before the final formulation of general relativity, Einstein remarked 
in a discussion with Max Born that, in the weak-field limit, gravitational 
waves exist and propagate with the velocity of light (Poincare pioneered the 
idea of gravitational waves propagating with the velocity of light in 1905 - 
see [|l^). Yet, in 1936 Einstein wrote to Born |247| : "... gravitational waves 
do not exist, though they had been assumed a certainty to the first approxi- 
mation. This shows the nonlinear general relativistic field equations can tell 
us more, or, rather, limit us more than we have believed up to now. If only 
it were not so damnably difficult to find rigorous solutions" . However, after 
finding a mistake in his argumentation (with the help of H. Robertson) and 



discovering with Nathan Rosen cylindrical gravitational waves [248| as the 
first exact radiative solutions to his vacuum field equations, Einstein changed 
his mind. In fact, cylindrical waves were found more than 10 years before Ein- 
stein and Rosen by Guido Beck in Vienna [ |249[ . Beck was mainly interested 
in time-independent axisymmetric Weyl fields, but he realized that through 
a complex transformation of coordinates (z — > li, t — > iz) one obtains cylin- 
drically symmetric time-dependent fields which represent cylindrical gravi- 
tational waves, and wrote down equations ( |7l| ) and (^ below. The work 
of Einstein and Rosen is devoted explicitly to gravitational waves. It inves- 
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tigates conditions for the existence of standing and progressive waves, and 
even notices that the waves carry away energy from the mass located at the 
axis of symmetry. We shah thus not modify the tradition and will call this 
type of waves Einstein- Rosen waves (which some readers may wish to shorten 
to EROS-waves). 

This type of waves, symmetric with respect to the transformation z —z 
(z - the axis of symmetry), contains one degree of freedom of the radiation 
field and corresponds to a fixed state of polarization. The metric can be 
written in the form 



(65) 



where p and t are invariants ("Weyl-type canonical coordinates"), and ip = 
ip{t,p), 7 = j{t,p). The Killing vectors d/dip and d/dz are both spacelike 
and hypersurface orthogonal. 

The metric containing a second degree of freedom was discovered by 
Jiirgen Ehlers (working in the group of Pascual Jordan), who used a trick 
similar to Beck's on the generalized (stationary ) We yl metrics, and indepen- 
dently b y Kompaneets (see the discussion in |25(| ). In the literature (e.g. 
[ 251 252|] ) one refers to the Jordan-Ehlcrs-Kompanects form of the metric: 



{~dt^ + dp^) + e^"^ {dz + ujdipf + p^e-^'f'dip^ 



(66) 



Here, the additional function Ld{t,p) represents the second polarization. 

Despite the fact that cylindrically symmetric waves cannot describe ex- 
actly the radiation from bounded sources, both the Einstein-Rosen waves and 
their generalization ( |66[ ) have played an important role in clarifying a nu mbe r 
of complicated issues, such as the energy loss due to gr avitational waves 1 253 1 , 
the interaction of waves with cosmic strings [254 255 [, the asym ptotic struc- 
ture of radiative spacetimes [25C[, the dispersion of waves [256[, testing the 
quasilocal mass-energy [257|, testi ng co des in numerical relativity [251|, inves- 
tigation of the cosmic censorship [P58[ , and quantum gravity in a simplifi ed 



but field theoretically interesting context of midisuperspaces [259,260 261 



In the following we shall discuss in some detail the asymptotic structure 
and midisuperspace quantization since in these two issues cylindrical waves 
have played the pioneering role. Some other applications of cylindrical waves 
will be briefly mentioned at the end of the section. 



9.1 Cylindrical waves and the asymptotic structure 
of 3-dimensional general relativity 



In recent work with Ashtekar and Schmidt |262 263 1, which started thanks 
to the hospitality of Jiirgen Ehlers' group, we considered gravitational waves 
with a space-translation Killing field ("generalized Einstein- Rosen waves"). 
In (2-fl)-dimensional framework the Einstein-Rosen subclass forms a sim- 
ple instructive example of explicitly given spacetimes which admit a smooth 
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global null (and timelike) infinity even for strong initial data. Because of the 
symmetry, the 4-dimensional Einstein vacuum equations are equivalent to 
the 3-diniensional Einstein equations with certain matter sources. This re- 
sult has roots in the classical paper by Jordan, Ehlers and Kundt ||5^ which 
includes "reduction formulas" for the calculation of the Riemann tensor of 
spaces which admit an Abelian isometry group. 

Vacuum spacetimes which admit a spacelike, hypersurface orthogonal 
Killing vector d/dz can be described conveniently in coordinates adapted 
to the symmetry: 

ds'^ = V'^{x)dz'^ +gabix)dx''dx\ a, 5, . . . = 0, 1, 2, (67) 

where x = x'^ and (jat is a metric with Lorentz signature. The field equations 
can be simplified if one uses a metric in the 3-space which is rescaled by the 
norm of the Killing vector, and writes the norm of the Killing vector as an 
exponential. Then (^^ becomes 

ds2 ^ e2^(^)d22 ^ e-^'f''-''^gab{x)dx''dx\ (68) 
and the field equations, 

Rab - 2VaV'VfcV = 0, 5°'VaVfcV = 0, (69) 

where V denotes the derivative with respect to the metric gab, can be rein- 
terpreted as Einstein's equations in 3 dimensions with a scalar field (P = v^V" 
as source. Thus, 4-dimensional vacuum gravity is equivalent to 3-dimensional 
gravity coupled to a scalar field. In 3 dimensions, there is no gravitational 
radiation. Hence, the local degrees of freedom are all contained in the scalar 
field. One therefore expects that Cauchy data for the scalar field will suffice to 
determine the solution. For data which fall off appropriately, we thus expect 
the 3-dimensional Lorentzian geometry to be asymptotically flat in the sense 
of Penrose |27,264|, i.e. that there should exist a 2-dimensional boundary 



representing null infinity. 

In general cases, this is analyzed in |262 . Here we shall restrict ourselves 
to the Einstein-Rosen waves by assuming that there is a further spacelike, 
hypersurface orthogonal Killing vector d/dip which commutes with d/dz. 
Then, as is well known, the equations simplify drastically. The 3-metric is 
given by 

da^ = gabdx'^dx'' = e'^^i-dt^ + dp^) + p^dif^, (70) 
the field equations (|69| ) become 

y = p(V>2 + V''2), i^2p^^p', (71) 



and 



(72) 
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where the dot and the prime denote derivatives with respect to t and p 
respectively. The last equation is the wave equation for the non-flat 3-metric 
( [to] ) as well as for the flat metric obtained by setting 7 = 0. 

Thus, we can first solve the axisymmetric wave equation ( [72[ ) for ip on 
Minkowski space and then solve ( |7l| ) for 7 - the only unknown metric coeffi- 
cient - by quadratures. The "method of descent" from the Kirchhoff formula 
in 4 dimensions gives the representation of the solution of the wave equation in 
3 dimensions in terms of Cauchy data <fo — iplt — 0, x, y),^i — ip.tii — 0, x, y) 



(see |262|). We assume that the Cauchy data are axially symmetric and of 
compact support. 

Let us look at the behaviour of the solution at future null infinity J^^ . 
Let p, if be polar coordinates in the plane, and introduce the retarded time 
coordinate u = t—p to explore the fall-off along constant u null hypersurfaces. 
For large p, the function ip at u ^ constant admits a power series expansion 
in p~^: 

The coefficients in this expansion are determined by integrals over the Cauchy 
data. At M 3> PO) Po being the radius of the disk in the initial Cauchy surface 
in which the data are non-zero, we obtain 

Mu) = ^ + ^ + ..., (74) 

where fco , ki are constants which are determined by the data. If the solution 
happens to be time-symmetric, so that vanishes, we find /o ^ u^^ for 
large u. Similarly, we can also study the behaviour of the solution near the 
timelike infinity of 3-dimensional Minkowski space by setting t — U + up, 



K > 1, and investigating for p ^ 00 with U and k fixed. We refer to [262| 
for details. 

In Bondi-type coordinates {u = t — p, p,ip), equation ( ^0[ ) yields 

da^ = e'^^{~du^ ~2dudp) + p^dip^. (75) 
The Einstein equations take the form 

7,„ = 2pV,«(V',p-V'.«), 7,p = pV'!^p, (76) 
and the wave equation on ip becomes 

-2i;^up + ip,pp + P~\i',p-'ip,u) = 0- (77) 



The asymptotic form of ip(t, p) is given by the expansion ([73[). Since we can 
differentiate (|73| ) term by term, the field equations ( [76|) and (|77| ) imply 

9k{u) 

k=l 



7.„ = ~2[/o(u)p+5:^, (78) 
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where the functions gk,hk are products of the functions fa, fk, fo, fk- Inte- 
grating (^9|) and fixing the arbitrary function of u in the resuh using ([ts]), 
we obtain 

Thus, 7 also admits an expansion in p^^, where the coefficients depend 
smoothly on u. It is now straightforward to show that the spacetime ad- 
mits a smooth future null infinity, J^"*". Setting p = p~^,u = u,ip = tp and 
rescaling gab by a conformal factor f] = p, we obtain 

da^ = Q^da^ = e^''{~p^du^ + 2dudp) + #^ (81) 

where ^(u, p) = 7(u, p~^)- Because of (^o|), 7 has a smooth extension through 
p = 0. Therefore, gab is smooth across the surface p — Q. This surface is the 
future null infinity, . Hence, the (2-|-l)-dimensional curved spacetime has 
a smooth (2-dimensional) null infinity. Penrose's picture works for arbitrarily 
strong initial data 'Fi. 

Using (^ij), we find that at we have: 

7(«, 00) = 70 - 2 / f^du. (82) 



Since one can make sure that 7 = at i+ [263|, one finds the simple result 
that 

/+00 
fodu. (83) 
-00 

At spatial infinity {t = constant, p — > 00), the metric is given by 

da"^ = e^^">{^dt^ + dp^) + p^dif^. (84) 

For a non-zero data, constant 70 is positive, whence the metric has a "coni- 
cal singularity" at spatial infinity. This conical singularity, present at spatial 
infinity, is "radiated out" according to equation ( |8^ . The future timelike in- 
finity, i^, is smooth. In (2-|-l)-diniensions, modulo some subtleties [ |262| |, equa- 
tion ( ^^ plays the role of the Bondi mass loss formula in (3-|-l)-dimensions, 
relating the decrease of the total (Bondi) mass-energy at null infinity to the 
flux of gravitational radiation. We can thus conclude that cylindrical waves 
in (2+l)-dimensions give an explicit model of the Bondi-Penrose radiation 
theory which admits smooth null and timelike infinity for arbitrarily strong 
initial data. There is no other such model available. The general results on the 
existence of in 4 dimensions, mentioned at the end of Section 1.3, assume 
weak data. 
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It is of interest to investigate cylindrical waves also in a (3+l)-diniensional 
context. The asymptotic behaviour of these waves was discussed by Stachel 
[ ^50| many years ago. However, his work deals solely with asymptotic direc- 
tions, which are perpendicular to the axis of symmetry, i.e. to the d/dz - 
Killing vector. Detailed calculations show that, in contrast to the perpendic- 
ular directions, where null infinity in the (3+l)-dimensional framework does 
not exist, it does exist in other directions for data of compact support. If the 
data are not time-symmetric, the fall-off is so s low that (l ocal ) null infinity 
has a poly homogeneous (logarithmic) character |18C| ~ see |263tl for details. 

We have concentrated on the simplest case of Einstein- Rosen waves. They 
served as a prototype for developing a general framework to analyze the 
asymptotic structure of spacetime at null infinity in three spacetime dimen- 
sions. This structure has a number of quite surprising featu res w hich do not 
arise in the Bondi-Penrose description in four dimensions [262|. One of the 
motivations for developing such a framework is to provide a natural point 
of departure for constructing the stage for asymptotic quantization and the 
S-matrix theory of an interesting midisuperspace in quantum gravity. 



9.2 Cylindrical waves and quantum gravity 

As the editors of the Proceedings of the 117th WE Heraeus Seminar on 
canonical gravity in 1993 [ p65[ , Jiirgen Ehlers and Helmut Fricdrich start 
their Introduction realistically: "When asking a worker in the field about the 
progress in quantum general relativity in the last decade, one shouldn't be 
surprised to hear: 'We understand the problems better'. If it referred to a 
lesser task, such an answer would sound ironic. But the search for quantum 
gravity... has been going on now for more than half a century and in spite 
of a number of ingenious proposals, a satisfactory theory is still lacking..." 
Although I am following the subject from afar, I believe that one would 
not be too wrong if one repeated the same words in 1999. However, apart 
from general theoretical developments, many interesting quantum gravity 
models have been studied, and exact solutions have played a basic role in 
them. In particular, in the investigations of (spherical) gravitational collapse 
and in quantum cosmology based typically on homogeneous cosmological 
models (cf. S ection 12.1), one star ts from simple classical solutions - see e.g. 
[266 267l , 268| for reviews and |269| for a bibliography up to 1990. A common 



feature of such models is the reduction of infinitely many degrees of freedom 
of the gravitational field to a finite number. In quantum field theory (such 
as quantum electrodynamics) a typical object to be quantized is a wave with 
an infinite number of degrees of freedom. The first radiative solutions of 
the gravitational field equations which were subject to quantization were 



the Einstein-Rosen waves. Kuchaf [259| applied the methods of canonical 



quantization of gravity to these waves, using the methods employed earlier in 
the minisuperspace models, i.e. restricting himself only to geometries (fields) 
preserving the symmetries. 
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The Einstein-Rosen cylindrical waves have an infinite number oo^ of de- 
grees of freedom contained in one polarization, one degree of freedom for each 
cylindrical surface drawn around the axis of symmetry. Moreover, the slicing 
of spacetime by spacelike (cylindrically symmetric) hypersurfaces is not fixed 
completely by the symmetry - an arbitrary cylindrically symmetric defor- 
mation of a given slice leads again to an allowed slice. Such a deformation 
represents an cx3^ "fingered time". Hence, the resulting space of 3-geometries 
on cylindrically symmetric slices is infinitely richer than the minisuperspaces 
of quantum cosmology. The exact Ein stein-Rosen waves thus inspired the 



first construction of what Kuchaf [259| called the "midisuperspace^\ 



Let us briefly look at the main steps in Kuchaf 's procedure. The sym- 
metry of the problem implies that the spatial metric has the form 

511 = e^"*, 522 = i?^e"*, ^33 = e*, (85) 

where 7, <P, and R are functions of a single cylindrical coordinate = r {x^ — 
(fi, — z). Similarly the lapse function N — N(r) depends only on r, and the 
shift vector has the only nonvanishing radial component — N^{r),N'^ = 
= 0. We have adopted here Kuchaf's notation. When we put R — r = 
<P = 2iIj, 7 — > 27, N — e^~'^, and = 0, we recover the standard 
Einstein-Rosen line element (^5|); however, in general the radial and time 
coordinates t and r differ from the canonical Einstein- Rosen radial and time 
coordinates in which the metric has the standard form (^). The symmetry 
implies that the canonical momentum tt**^ is diagonal and expressible by three 
functions 7r-y, tt^, 7r<g of r; for example, tt^^ = 'K~^e^~^ , and similarly for the 
other components. After the reduction to cylindrical symmetry, the action 
functional assumes the canonical form 



/oo 1*00 
dt I dr(7r^7 + ttrR + TT.p^ - NH - N^Hi) , 
'OO Jq 



(86) 



in which 7,i?, ^ are the canonical coordinates and 7r-y, tt^j, tt^ the conjugate 
momenta (the integration over z has been limited hy z = zq and z = zq + 1). 
The superhamiltonian Ti. and supermomentum Tii arc rather complicated 
functions of the canonical variables: 

n = ei(*-^) (-^^^fl + + 2R" - 7'i?' + ii?<?'2) , (87) 

Hi = -27r; + j't:^ + R'ttr + ^'t:^. (88) 

The most important step now is the replacement of the old canonical 
variables 7, tt^, i?, tt^ by a new canonical set T, Ut, i?, IIji through a suitable 
canonical transformation. We shall write here only one of its components (see 



■^^ For the basic concepts and ide as o f canonical gravity, we refer to e.g. and 



especially to Kuchai^'s review |27C], where the canonical quantization of cylindri- 
cal waves is also analyzed. 
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[259 27C] for the complete transformatfon): 

T(r) = T{oo) + f [~TT^{r)] dr. (89) 

J oo 

By integrating the Hamilton equations following from the action (^6|), rewrit- 
ten in the new canonical coordinates, one finds that T and R are the Einstein- 
Rosen privileged time and radial coordinates, i.e. those appearing in the 
canonical form ( |65| ) of the Einstein- Rosen metric (with T = t, R = p). Ac- 
cording to (|8S|), the Einstein- Rosen time can be reconstructed, in a non-local 
way, from the momentum tt^, which characterizes the extrinsic curvature of 
a given hypersurface. In this way, the concept of the "extrinsic time repre- 
sentation" entered canonical gravity with cylindrical gravitational waves. 

In terms of the new canonical variables, the super hamiltonian and super- 
momentum become 

ji = ei(*-7) (^'Tjj, + T'nR + i {R-\l + i?<P'2)) , (90) 
Hi = T'Ut + R'Ur + (^'TTj,. (91) 

Since TL and TLi are linear in Ut and i7jj, the classical constraints Ti — 
0, Til — can immediately be resolved with respect to these momenta, con- 
jugate to the "embedding" canonical variables T{r) and R{r): 

- TIt - (i?'' - T'Y^ [\{R'^^l + m'^)R' - ^'tt^T'] = 0, (92) 



and similarly for 11 ji. It is easy to see [p59y27Cl(| that the constraints have the 
same form as the constraints for a massless scalar field ^ propagating on a flat 
background foliated by arbitrary spacelike hypersurfaces T = T{r), R = R{r). 
The canonical variables tt^ represent the true degrees of freedom, and 
the remaining canonical variables play the role of spacelike embeddings of a 
Cauchy hypersurface into spacetime. 

After turning the canonical momenta Ut, Uj^, tt^, into variational deriva- 
tives, e.g. Ut — —iS/6T{r), one can impose the classical constraints H = 
0,^1 = as restrictions on the state functional ^{T,R,(l>): Wl/ = 0, Hi^ = 
0. In particular, the Wheeler-DeWitt equation T-LW = in the extrinsic time 
representation assumes the form of a many-fingered time counterpart of an or- 
dinary Schrodinger equation. This reduces to the ordinary Schrodinger equa- 
tion for a single massless scalar field in Minkowski space if we adopt the 



standard foliation T — constant (see [259,270 for details) 



The described procedure, first realized in the case of the Einstein-Rosen 
waves, has opened a new route in canonical and quantum gravity. In contrast 
to the Arnowitt-Deser-Misner approach, in which the gravitational dynamics 
is described relative to a fixed foliation of spacetime, in this new approach 
(called "bubble time" dynamics of the gravitational field or the "internal time 



formalism" [271[) one tries to extract the many-fingered time (i.e. embeddings 



of Cauchy hypersurfaces) from the gravitational phase space, but does not fix 
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the foliation in the "target manifold" by coordinate conditions. However, the 
definition of the target manifold by a gauge (coordinate) condition is needed. 

This new approach has been so far successfully applied to a few other 
models (based on exact solutions) with infinite degrees of freedom, for ex- 
ample, plane gravitational waves, bosonic string, and as late as 1994, to 



spherically symmetric vacuum gravitational fields |272 . The internal time 



formalism for spacetimes with two Killing vectors was developed in [273| 
(therein references to previous works can also be found). Recently, canonical 
transformation techniques have been appli ed to Hamiltonian spacetime dy- 
namics with a thin spherical null-dust shell [ 274 1 . One would like to construct 
a midisuperspace model of spherical gravitational collapse, or more specifi- 
cally, a model for Hawking radiation with backreaction. The extensive past 
work on Hamiltonian approaches to spherically symmetric geometries (see 
[ ^74| for more than 40 references in this context) have not yet led to con- 
vincing insights. The very basic question of existence of the "internal time" 
formalism in a general situation has been most recently addressed by Hajicek 
; the existence has been proven, and shown to be related to the choice of 
gauge. 



9.3 Cylindrical waves: a miscellany 



Chandrasekhar [ p47[ constructed a formalism for cylindrical waves with two 
polarizations (cf. the metric (]6^)), similar to that used for the discussion of 
the collision of plane-fronted waves (Section 8.3). He obtained the "cylindri- 
cal" Ernst equation and corroborated (following the suggestion of O. Reula) 
the physical meaning of Thome's C- energy [p53| - the expression for energy 
suggested for cylindrical fields - by defining a Hamiltonian density corre- 
sponding to the Lagrangian density from which the Ernst equation can be 
derived. A brief summary of older work on the mass loss of a cylindrical source 
radiating out cylindrical waves and its relation to the C-energy is given in 
[ |65| . It should be pointed out, however, that although C-energy is a useful 
quantity, it was constructed by exploiting the local field equations, without 
direct reference to asymptotics. The physical energy (per unit z length) at 
both spatial and null infinity, which is the generator of the time translation, 
is in fact a non-polynomial function of the C-energy. In the weak field limit 



the two agree, but in strong fields they are quite different |262|. 

In | p56t , an exact solution was constructed with which one can study the 
dispersion of waves: a cylindrical wave packet, which though initially impul- 
sive, after reflection at the axis disperses, and develops shock wave fronts 
when the original wave meets the waves that are still ingoing. Cylindrical 
waves have been also analyzed in the context of phase shifts occurring in 



gravitational soliton interactions (see 1 275 and references therein) 



An exact explicit solution for cylindrical waves with two degrees of polar- 



ization has been obtained [252| from the Kerr solution after transforming the 



metric into "cylindrical" coordinates and using the substitution t iz, z 
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it, a — > ia. Both this solution and the well-known Weber- Wheeler-Bonnor 



pulse I pql have been employed as test beds in numerical relativity |276 , in 
particular in the approach which combines a Cauchy code for determining 
the dynamics of the central source with a characteristic code for determining 
the behaviour of radiation lESlI . 



In a number of works cylindrical waves have been considered in interaction 



with cosmic strings 1 254, 255 1. The strings are usually modelled as infinitely 
thin conical singularities. Recently Colombeau's theory of generalized func- 
tions was used to calculate the distributional curvature at the axis for a 



time-dependent cosmic string [ 277 | 



A somewhat surprising result concerning cosmic strings and radiation the- 
ory should also be noted: although an infinite, static cylindrically symmet- 
ric string does not, of course, radiate, it generates a no nvanishin g (though 



"non-radiative") contribution to the Bondi news function [ 278, 279 1. Recently, 
the asymptotics at null infinity of cylindrical waves with both polarizations 
(and, in general, an infinite cosmic string along the axis) has been analyzed 
in the context of axisymmetric electrovacuum spacetimes with a translational 



Killing vector at null infinity 1 28C | 



Finally, the cylindrically symmetric electrovacuum spacetimes with both 
polarizations, satisfying certain completeness and asymptotic flatness con- 
ditions in spacelike directions have been shown rigorously to imply that 
strong cosmic censorship holds ]25^ . This means that for generic (smooth) 
initial data the maximal globally hyperbolic development of the data is inex- 
tendible (no Cauchy horizons as for example, those discussed in Section 3.1 
for the Reissner-Nordstrom spacetime arise). This global existence result is 
non-trivial since with two polarizations and electromagnetic field present, all 
field equations arc nonlinear. 

10 On the Robinson- Trautman solutions 

Robinson- Trautman metrics are the general radiative vacuum solutions which 
admit a geodesic, shearfree and twistfree null congruence of diverging rays. 



In the standard coordinates the metric has the form |281] 



ds^ = 2r^p-^dCdC -2dudr - [ZilnP - 2r(lnP),„ - 2mr~^] du"^, (93) 

where C is a complex spatial (stereographic) coordinate (essentially 9 and 
If), r is the affine parameter along the rays, m is a retarded time, m is a 
function of u (which can be in some cases interpreted as the mass of the 
system), A = 2P'^{d'^ /d(dC), and P = P{uX,C) satisfies the fourth-order 
Robinson- Trautman equation 

Z\Z\(lnF) + 12 m (lnP),„ -4m,„ = 0. (94) 

The best candidates for describing radiation from isolated sources are the 
Robinson- Trautman metrics of type //with the 2-surfaces given by u, r = 
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constant and having spherical topology. The Gaussian curvature of can 
be expressed a.s K = AlnP. If K = constant, we obtain the Schwarzschild 
solution with mass equal to K^^ , 

These spacetimes have attracted increased attention in the last decad e - 
most recently in the work by Chrusciel, and Chrusciel and Singleton |282|. In 
these studies the Robinson- Trautman spacetimes have been shown to ex- 
ist globally for all positive "times" , and to converge asymptotically to a 
Schwarzschild metric. Interestingly, the extension of these spacetimes across 
the "Schwarzschild-like" event horizon can only be made with a finite degree 
of smoothness. All these rigorous studies are based on the derivation and 
analysis of an asymptotic expansion describing the long-time behaviour of 
the solutions of the nonlinear parabolic equation 




Fig. 11. The evolution of the cosmological Robinson- Trautman solutions with a 
positive cosmological constant. A black hole with the horizon is formed; at 
future infinity J"^ the spacetime approaches a de Sitter spacetime exponentially 
fast, in accordance with the cosmic no-hair conjecture. 



In our recent work 
spacetimes with a positive cosmological constant A. The results proving the 
global existence and convergence of the solutions of the Robinson- Trautman 
equation ( p4[ ) can be taken over from the previous studies since A does not 
explicitly enter this equation. We have shown that, starting with arbitrary, 
smooth initial data at u = uq (see Fig. 11), these cosmological Robinson- 
Trautman solutions converge exponentially fast to a Schwarzschild-de Sitter 
solution at large retarded times {u — > oo). The interior of a Schwarzschild- 
de Sitter black hole can be joined to an "external" cosmological Robinson- 
Trautman spacetime across the horizon Ti.'^ with a higher degree of smooth- 
ness than in the corresponding case with yl = 0. In particular, in the extreme 
case with 9Am^ = 1, in which the black hole and cosmological horizons coin- 
cide, the Robinson- Trautman spacetimes can be extended smoothly through 



163,P83| we studied Robinson- Trautman radiative 
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7i+ to the extreme Schwarzschild-de Sitter spacetime with the same values 
of A and m. However, such an extension is not analytic (and not unique). 

We have also demonstrated that the cosmological Robinson- Trautman so- 
lutions represent explicit models exhibiting the cosmic no-hair conjecture: a 
geodesic observer outside of the black hole horizon will see, that inside his 
past light cone, these spacetimes approach the de Sitter spacetime exponen- 
tially fast as he approaches the future (spacelike) infinity For a freely 
falling observer the observable universe thus becomes quite bald. This is what 
the cosmic no-hair conjecture claims. As far as we are aware, these models 
represent the only exact analytic demonstration of the cosmic no-hair con- 
jecture under the presence of gravitational waves. They also appear to be 
the only exact examples of black hole formation in nonspherical spacetimes 
which are not asymptotically flat. Hopefully, these models may serve as tests 
of various approximation methods, and as test beds in numerical studies of 
more realistic situations in cosmology. 



11 The boost-rotation symmetric radiative spacetimes 

In this section we would like to describe briefly the only explicit solutions 
available today which are radiative and represent the fields of finite sources. 
Needless to say, we cannot hope to find explicit analytic solutions of the 
Einstein equations without imposing a symmetry. A natural first assumption 
is axial symmetry, i.e. the existence of a spacelike rotational Killing vector 
d/dif. However, it appears hopeless to search for a radiative solution with 
only one symmetry. We are now not interested in colliding plane waves since 
these do not represent finite sources; we wish our spacetime to be as "asymp- 
totically flat as possible" . The unique role of the boost-rotation symmetric 
spacetimes is exhibited by a theorem, formulated precisely and proved for 



the vacuum case with hypersurface orthogonal Killing vectors in |284|, and 



generalized to electrovacuum spacetimes with Killing vectors which need not 



be hypersurface orthogonal in |279| (see also references therein). This the- 
orem roughly states that in axially symmetric, locally asymptotically flat 
spacetimes (in the sense that a null infinity satisfying Penrose's requirements 
exists, but it need not necessarily exist globally), the only additional symme- 
try that does not exclude radiation is the boost symmetry. 

In Minkowski spacetime the boost Killing vector has the form 

d d 

(boost = z— + t—, (95) 
ot az 

so that orbits of symmetry to which the Killing vector is tangent are hy- 
perbolas — = B = constant, x, y — constant. Orbits with B > are 
timelike; they can represent worldlines of uniformly accelerated particles in 
special relativity. Imagine, for example, a charged particle, axially symmet- 
ric about the z-axis, moving with a uniform acceleration along this axis. 
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The electromagnetic field produced by such a source will have boost-rotation 
symmetry. 

Figure 12 shows two particles uniformly accelerated in opposite directions 
along the z-axis. In the space diagram (left), the "string" connecting the par- 
ticles is also indicated. In the spacetime diagram, the particles' worldlines are 
shown in bold. Thinner hyperbolas represent the orbits of the boost Killing 
vector ( p5| ) in the regions > where it is spacelike. In Figure 13 the 
corresponding compactified diagram indicates that null infinity cannot be 
smooth everywhere since it contains four singular points in which particles' 
worldlines "start" and "end" . Notice that in electromagnetism the presence 
of two particles, one moving along z > 0, the other along z < 0, makes the 
field symmetric also with respect to inversion z — > — z. The electromagnetic 
field can be shown to be analytic everywhere, except for the places where the 
particles occur. These two particles move independently of each other, since 
their worldlines are divided by two null hypersurfaces z = <, z = —t. This is 
analogous to the boost-rotation symmetric spacetimes in general relativity 
that we are now going to discuss. 

Specific examples of solutions representing "uniformly accelerated parti- 
cles" have been analyzed for 35 years, startin g wi th the first solutions o f thi s 
type obtained by Bonnor and Swaminarayan |285|, and Israel and Khan |286|. 
In a curved spacetime the "uniform acceleration" is understood with respect 
to a fictitious Minkowski background, and the "particles" mean singularities 
or black holes. For a more extensive description of the history of these spe- 
cific solutions discovered before 1985, see [ ^87| . From a unified point of view, 
boost-rotation symmetric spacetimes (with hypersurface orthogonal Killing 
vectors) were defined and treated geometrically in |28S|. We refer to this de- 
tailed work for rigorous definitions and theorems. Here we shall only sketch 
some of the general properties and some applications of these spacetimes. 

The metric of a general boost- rotation symmetric spacetime in "Cartesian- 
type" coordinates {t,x,y,z} reads: 

ds^ = -Y^-^ [{e^x^ + e-i'y^)dx^ + 2xy{e^ - e-'')dxdy] + 
X -\- y 

+ ^^{eV + e-^x')dy^ + -^{e^z^ ~ en')dz' - 
x'^ +2; z'^ — 

5^-^ \2zt{e^ - e^dzdt + (e^z^ - eh^)dt^] , (96) 

where ^ and A are functions of ~ x^ + j/^ and z^ — . As a consequence 
of the vacuum Einstein equations, the function must satisfy an equation 
of the form which is identical to the flat-space wave equation; and function 
A is determined in terms of /i by quadrature. Now it can easily be seen that 
the metric (^6|) admits axial and boost Killing vectors which have exactly the 
same form as in Minkowski space, i.e. the axial Killing vector d/dLp and the 
boost Killing vector (|5|). In fact, the whole structure of group orbits in boost- 
rotation symmetric curved spacetimes outside the sources (or singularities) 
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is the same as the structure of the orbits generated by the axial and boost 
KiUing vectors in Minkowski space. In particular, the boost Killing vector 
( p5[ ) is timelike in the region > . The invariance of a metric (or of 
any other field) in a time-direction (determined in a coordinate-free manner 
by a timelike Killing vector) means stationarity, and of course, we could 
hardly expect to find radiative properties there. Intuitively, the existence of 
a timelike Killing vector in the region > P is understandable because 
there (generalized) uniformly accelerated reference frames can be introduced 
in which sources are at rest, and the fields are time independent. 




Fig. 12. Two particles uniformly accelerated in opposite directions. The orbits of 
the boost Killing vector (thinner hyperbolas) are spacelike in the region > z'^. 



However, in the other "half" of the spacetime, > z^, the boost Killing 
vector (95) is spacelike (see the lines representing orbits of the boost Killing 
vector in Fig. 12). Hence in this region the metric ( |9^ ) is nonstationary. 
Here we expect to discover radiative properties. Indeed, it can be shown that 
for t'^ > z'^ + the metric ( p6| ) can locally be transformed into the metric 
of Einstein-Rosen cylindrical waves. Although locally in the whole region 
> z^ the metric (^) can be transformed into a radiative metric, the global 
properties of the boost-rotation symmetric solutions are quite different from 
those of cylindrical waves. Again, we have to refer to the work |28S] for a 
detailed analysis. Let us only say that the boost-rotation symmetric solutions, 
if properly defined - with appropriate boundary conditions on functions A and 
\i - always admit asymptotically fiat null infinity J at least locally. Starting 
with arbitrary solutions A and /i, and adding suitable constants to both A and 
[i (Einstein's equations are then still satisfied), we can always guarantee that 
even global J exists in the sense that it admits smooth spherical sections. 
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For the special type of solutions for A and /i, complete J satisfies Penrose's 
requirements, except for four points in which the sources "start" and "end" 
(cf. Fig. 13). In aU cases one finds that the gravitational field in smooth 
regions of the null infinity is radiative [ 
term of the Riemann curvature tensor, proportional to r^^ (where — 
yO^ + z^), is nonvanishing and has the same algebraic structure as the Riemann 
tensor of plane waves. This is fully analogous to the asymptotic properties 
of radiative electromagnetic fields outside finite sources. Recently, general 
forms of the news functions have been obtained for electrovacuum spacetimes 
with boost-rotation symmetry and with Killing vectors which need not be 
hypersurface orthogonal [^79| . 
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Fig. 13. The Penrose compactified diagram of a boost-rotation symmetric space- 
time. Null infinity can admit smooth sections. 

It is well known that in general relativity the "causes" of motion are always 
incorporated in the theory - in contrast to electrodynamics where they need 
not even be describable by Maxwell's theory. In a general case of the boost- 
rotation symmetric solutions there exist nodal (conical) singularities of the 
metric distributed along the z-axis which can be considered as "strings" , and 
cause particles to accelerate. They reveal themselves also at J. However, the 
distribution of nodes can always be arranged in such a manner that admits 
smooth regular sections as mentioned above. 

In exceptional cases, when J7 is regular except for four points, either 
the particles are "self-accelerating" due to their "inner" multipole structure, 



271; 



In particular, the leading 
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which has to include a negative mass; or there are more particles distributed 
along z > (and symmetrically along z < 0) with the signs and the magni- 
tudes of their masses and accelerations chosen appropriately. (For the concept 
of a negative mass in general relativity, and the first discussion of a "chas- 
ing" pair of a positive and a negative mass particle, see classical papers by 
Bondi, and Bonnor and Swaminarayan [285|.) An infinite number of differ- 
ent analytic solutions representing self-accelerating particles was constructed 
explicitly |290|. Although a negative mass cannot be bought easily in the 
shop (as Bondi liked to say), these solutions are the only exact solutions of 
Einstein's equations available today for which one can find such quantities 
of physical interest as radiation patterns (ang ular distribution of gravita- 
tional radiation) , or total radiation powers 1 28^ ] . From a mathematical point 
of view, these solutions represent the only known spacetimes in which ar- 
bitrarily strong (boost-rotation symmetric) initial data can be chosen on a 
hyperboloidal hypcrsurface in the region > z^, which will lead to a com- 
plete smooth null infinity and a regular timelike future infinity. With these 
specific examples one thus docs not have to require weak-field initial data 
as one has to in the work of Fricdrich, and Christodoulou and Klaincrman, 
mentioned at the end of Section 1.3. 

The boost-rotation symmetric radiative spacetimes can be used as test 
beds for approximation methods or numerical relativity. Bicak, Reilly and 
Winicour [291| found the explicit boost-rotation symmetric "initial null cone 
solution" , which solves initial hypcrsurface and evolution equations in "ra- 
diative" coordinates employed in the null cone version of numerical relativity. 
This solution has been used for checking and improving numerical codes for 
computing gravitational radiation from more realistic sources; a new solution 
of this type has also been found [292 . Recently, the specific boost-rotation 



symmetric spacetimes constructed in 1 29C | were used as test beds in th e sta n- 



dard version of numerical relativity based on spacelike hypersurfaces |293|. 

There exist "generalized" boost-rotation symmetric spacetimes which are 
not asymptotically flat, but are of considerable physical interest. They de- 
scribe accelerated particles in asymptotically "uniform" external fields. One 
can construct such solutions from asymptotically flat boost-rotation symmet- 
ric solutions for the pairs of accelerated particles by a limiting procedure, in 
which one member of the pair is "removed" to infinity, and its mass param- 
eter is simultaneously increased [294|. Since the resulting spacetimes are not 
asymptotically flat, their radiative properties are not easy to analyze. Only 
if the external field is weak will there exist regions in which the spacetimes 
are approximately flat; and here their radiative properties might be investi- 
gated. So far no systematic analysis of these spacetime has been carried out. 
Nevertheless, they appear to offer the best rigorous examples of the motion 
of relativistic objects. No nodal singularities or negative masses are necessary 
to cause an acceleration. 

As an eloquent example of such a spacetime consider a charged (Reissner- 
Nordstrom) black hole with mass M and charge Q, immersed in an electric 
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field "uniform at infinity" , characterized by tlie field-strength parameter E. 
An exact solution of the Einstein-Maxwell equations exists which describes 
this situation |131|. ft goes over into an approximate solution obtained by 
perturbing the charged black hole spacetime by a weak external electric field 
which is uniform at infinity [295|. One of the results of the analysis of this 
solution is very simple: a charged black hole in an electric field starts to 
accelerate according to Newton's second law, Ala — QE, where all the quan- 
tities can be determined - and in principle measured - in an approximately 
flat region of the spacetime from the asymptotic form of the metric. Recall 
T. S. Eliot again: "There is only the fight to recover what has been lost / 
And found and lost again and again." 

These types of generalized boost-rotation symmetric spacetimes ("gener- 
alized C-metrics") have been used by Hawking, Horowitz, Ross, and others 
[ p96[ in the context of quantum gravity ~ to describe production of black hole 
pairs in strong background fields. 

Recently, we have studied the spinning C-metric discovered by Plebahski 
and Demianski |297|. Transformations can be found which bring this metric 
into the canonical form of spacetimes with boost-rotation symmetry |298|. 
The metric represents two uniformly accelerated, rotating black holes, either 
connected by conical singularity, or with conical singularities extending from 
each of them to infinity. The spacetime is radiative. No other spacetime of 
this type, with two Killing vectors which are not hypersurface orthogonal, is 
available in an explicit form. 



12 The cosmological models 

In light of Karl Popper's belief that "all science is cosmology" , it appears 
unnecessary to justify the choice of solutions for this last section. As in the 
whole article, these will be primarily vacuum solutions. On the other hand, 
in light of the light coming from about 10^^ galaxies, each with about 10^^ 
stars, it may seem weird to consider vacuum models of the Universe. Indeed, 
it has become part of the present-day culture that spatially homogeneous 
and isotropic, expanding Friedmann-Robertson- Walker (FRW) models, filled 
with uniformly distributed matter, correspond well to basic observational 
data. In order to achieve a more precise correspondence, it appears sufficient 
to consider just perturbations of these "standard cosmological models" . To 
explain some "improbable" features of these models such as their isotropy 
and homogeneity, one finds an escape in inflationary scenarios. These views 
of a "practical cosniologist" are, for example, embodied in one of the most 
comprehensive recent treatise on physical cosmology by Peebles . 

Theoretical (or mathematical) cosmologists, however, point out that more 
general cosmological models exist which differ significantly from a FRW 
model at early times, approach the FRW model very closely for a certain 
epoch, and may diverge from it again in the future. Clearly, the FRW uni- 
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verses represent only a very special class of viable cosmological models, 
though the simplest and most suitable for interpretations of "fuzzy" cos- 
mological observational data. 

Simple exact solutions play a significant role in the evolution of more 
general models, either as asymptotic or intermediate states. By an "interme- 
diate state" one means the situation when the universe enters and remains 
in a small neighbourhood of a saddle equilibrium point. A simple example 
is the Lemaitre matter-filled, homogeneous and isotropic model with a non- 
vanishing cosmological constant (see e.g. |2^), which expands from a dense 
state (the big bang, or "primeval atom" in Lemaitre's 1927 terminology), 
passes through a quasistatic epoch in which all parameters are close to those 
of the static Einstein universe (cf. Section 1.2), and then the universe ex- 
pands again. An "asymptotic state" means close either to an initial big bang 
(or possibly a final big crunch) singularity, or the situation at late times in 
forever expanding universes. It is easy to see that at late times in indefinitely 
expanding universes the matter density decreases, and vacuum solutions may 
become important. However, as we shall discuss below, vacuum models play 
an important role also close to a singularity, when the matter terms in Ein- 
stein's equations are negligible compared to the "velocity terms" (given by 
the rate of change of scale factors) or to the curvature terms (characteriz- 
ing the curvature of spacelike hypersurfaces). In particular, the pioneering 
(and still controversial) work started at the end of the 1950s by Lifshitz and 
Khalatnikov, and developed later on by Belinsky, Khalatnikov and Lifshitz, 
has shown that the fact that the presence of matter does not influence the 
qualitative behaviour of a cosmological model near a singularity has a very 



general significance (see [299 and [ pOO| for the main original references, and 
[ [139| for a brief review) . 

In gaining an intuition in the analysis of general cosmological singulari- 
ties, the class of spatially homogeneous anisotropic cosmological models have 
played a crucial role. These so called Bianchi models admit a simply transi- 
tive 3-dimensional homogeneity group. Among the Bianchi vacuum models 
there are special exact explicit solutions, in particular the Kasner and the 
Bianchi type II solutions, which exhibit some aspects of general cosmological 
singularities. The Bianchi models have also had an impact on other issues in 
general relativity and cosmology. 

Much work, notably in recent years, has been devoted to the class of 
both vacuum and matter-filled cosmological solutions which are homogeneous 
only on 2-dimensional spacelike orbits. Thus they depend on time and on 
one spatial variable, and can be used to study spatial inhomogeneities as 
density fluctuations or gravitational waves. The vacuum cosmological models 
with two spacelike Killing vectors, sometimes called the Gowdy models,f^ are 



In fact, by Gowdy models, one more often means only the cases with closed group 
orbits, with two commuting spacelike othogonally-transitive Killing vectors (the 
surface elements orthogonal to the group orbits are surface- forming). 
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interpreted as gravitational waves in an expanding (or contracting) universe 
with compact spatial sections. We shall discuss these two classes separately. 



12.1 Spatially homogeneous cosmologies 

The simplest solutions, the Minkowski, de Sitter, and anti de Sitter space- 
times, which have also been used in cosmological contexts (cf. Section 1.3), 
are A- dimensionally homogeneous. As noted in Section 8.1, the vacuum plane 
waves (see equations (51), (^), (|6l|)) are also homogeneous spacetimes; and 



since they can be suitably sliced by spacelike hypersurfaces with expand- 
ing normal congruence, they can become asymptotic states in homogeneous 
expanding cosmologies. There exist several important non- vacuum homoge- 
neous spacetimes, for example, the Einstein static universe (cf. Section 1.2), 
and Godel's stationary, rotating universe (see e.g. ^,301|), famous for the 



first demonstration that Einstein's equations with a physically permissible 
matter source are compatible with the existence of closed timelike lines, i.e. 
with the violation of causality. 

Here we shall consider models in which the symmetry group does not make 
spacetime a homogeneous space, but in which each event in spacetime is con- 
tained in a spatial hypersurface that is homogeneous. The standard FRW 
models represent a special case of such models (they admit, in addition, an 
isotropy group SO{i) at each point). The general spatially homogeneous so- 
lutions comprise of the Kantow ski- Sachs universes and a much wider class 
of Bianchi models. By definition, the Bianchi models admit a simply tran- 
sitive 3-dimensional homogeneity group Ga. There exist special "locally ro- 
tationally symmetric" (LRS) Bianchi models which admit a 4-dimensional 
isometry group G4 acting on homogeneous spacelike hypersurfaces, but these 
groups have a simply transitive subgroup G3. In contrast to this, Kantowski- 
Sachs spacetimes admit G4 (acting on homogeneous spacelike hypersurfaces) 
which does not have any simply transitive subgroup G3; it contains a mul- 
tiply transitive G3 acting on 2-dimensional surfaces of constant curvature, 
G4 — M X 50(3). A special case of the vacuum Kantowski-Sachs universe 
is represented by the Schwarzschild metric inside the horizon (with t and r 
interchanged). There has been a con tinui ng interest in the Kantowski-Sachs 
models since their discovery in 1966 [t302| , to which, as the authors acknowl- 
edge, J. Ehlers contributed by his advice. Some of these mode ls ha d already 



appeared in the PhD thesis of Kip Thorne in 1965 (see also 303 for mag- 



netic Kantowski-Sachs models). Here, however, we just refer the reader to 



[304 305] for their classical description, to [306| for a canonical and quan- 
tum treatment, and to [ |307| for the latest discussion of the Kantowski-Sachs 
quantum cosmologies. 

Although the 3-dimensional Lie groups which are simply transitive on 
homogenous 3-spaces were classified by Bianchi in 1897, the importance of 
Bianchi's work for constructing vacuum cosmological models was only dis- 



covered by Taub in 1951 |197|, when the Taub space (cf. Section 7) was first 
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given. It is less known that at approximately the same time, if not earlier, 
the first explicit spatially homogeneous expanding and rotating cosmological 
models with matter (of the Bianchi type IX) were constructed by Godel,^ 
who first presented his results at the International Congress of Mathematics 
held at Cambridge (Mass.) from August 30 till September 5, 1950. 

An exposition of Bianchi models has been given in a number of places: 
in the account on relativistic cosmology by Heckmann and Schiicking [|3G8[ 
(complementing the chapter on exact sol ution s by Ehlers and Kundt pW), 



i n th e monographs of Ryan and Shepley 1 304 , and Zel'dovich and Nov ikov 
[ ^09| , in several comprehensive surveys by MacCallum (see e.g. 305 and [t^lOj 
for his latest review containing a number of references), most recently, in the 
book on the dynamical system approach in cosmology (in the Bianchi models 
in particular) edited by Wainwright and Ellis |31l[ ; and, first but not least, in 
the classics of L andau and Lifshitz |l39| . The Hamiltonian approach initiated 
by Misner ||312| in 1968, and used in, amongst other things, the construction 
of various minisuperspace models in quantum gravity, has been reviewed by 



Ryan 266 ; for more recent accounts, see several contributions to Misner's 



Festschrift |313 . An interesting framework which unifies the Hamiltonian 
approach to the solutions which admit homogeneous hypersurfaces either 
spacelike (as Bianchi models) or timelike (as static spherical, or stationary 
cylindrical models) was recently developed by Uggla, Jantzen and Rosquist 



in |314] (with 115 references on many exact solutions). Herewith we shall only 
briefly introduce the Bianchi models, note their special role in understanding 
the character of an initial cosmological singularity, and mention some of the 
most recent developments not covered by the reviews cited above. 

The line element of the Bianchi models can be expressed in the form 

ds^ = -dt^ + gab{t) uj'' w^ (97) 

where the time-independent 1-forms {— E^dx"), a = 1,2,3, are dual to 
time-independentg spatial frame vectors (often an arbitrary time -variable 
t is introduced by dt — N{t) dt, N being the usual lapse function). Both 

Godel's profound ideas and results in cosmology, and their influence on later 
developments have been discussed in depth by G. Ellis in his lecture at the Godel 
'96 conference in Brno, Czech Republic, where Godel was born in 1906 (78 years 
before Godel, Ernst Mach was born in a plac e wh ich today belongs to Brno). 



In the extended written version of Ellis' talk [301| it is indicated that Godel's 
work also initiated the investigation of Taub. This may well be true with Godel's 
paper on the stationary rotating universe, but Taub's paper on Bianchi models 
was received by the Annals of Mathematics on May 15, 1959, i.e. before Godel's 
lecture on expanding and rotating models at the Congress of Mathematics took 
place. 

The gravitational degrees of freedom are associated with the component (scalar) 
functions gab{i) - the so called metric approach. Alternatively, in the orthonormal 
frame approach, one chooses gab{t) = Sat and describes the evolution by time- 
dependent forms oj". In still another approach one employs the automorphism 
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and Ea are group- invariant, commuting with the three KiUing fields which 
generate the homogeneity group. They satisfy the relations 

duj"" = - hjl^ J" A cj", (98) 
[E,,Eb] =C^fcE^ (99) 

where d is the exterior derivative and C^^ are the structure constants of the 
Lie algebra of the homogeneity group. The models are classified according to 
the possible distinct sets of the structure constants. They are first divided 
into two classes: in class A the trace C^^ — 0, and in class B, C^^ 7^ 0. In 
class A one can choose C^^ = n^^^Cobc (no summation over a), and classify 
various symmetry types by parameters n^"-* with values 0, ±1. In class B, in 
addition to n'^"\ one needs the value of a constant scalar h (related to C^^) 



to characterize types VI/j and VII^ (see e.g. |31l[ ). 

The simplest models are the Bianchi I cosmologies in class A with n*^"^ = 
0, i.e. C^^ = 0, so that all three Killing vectors (the group generators) com- 
mute. They contain the standard Einstein-de Sitter model with flat spatial 
hypersurfaces (curvature index fc = 0). In the vacuum case, all Bianchi I 
models are given by the well-known 1-parameter family of Kasner metrics 
(found in 1921 by E. Kasner and in 1933 by G. Lemaitre without considering 
the Bianchi groups) 

ds^ = -dt^ + t'^P'dx^ + t^P^dy^ + t^P'^dz"^, (100) 

where 

Pi+P2+P3 = l: pl+pl+pl = i- (101) 

These metrics were first used to investigate various effects in anisotropic cos- 
mological models. For example, in contrast to standard FRW models with 
"point-like" initial singularities, the Kasner metrics can permit the so called 
"cigar" and "pancake" singularities. To be more specific, consider the congru- 
ence of timelike lines with unit tangent vectors orthogonal to constant time 
hypersurfaces, and define the expansion tensor Oap by 6af3 — cTafS + ^Ohap, 
where hafi = gap + Jt-q"-/? is a projection tensor, (Jap — ni^a.;p) — \Qhap is the 
shear, and 9 = 9^ . Determining the three spatial eigenvectors of O^p with the 
corresponding eigenvalues 9i {i = 1,2,3), one can define the scale factors U 
by the relation 9i = {dli/dt)/li , and the Hubble scalar H = ^(^i -I- 6*2 -I- ^3). 
In the FRW models, all k at the big bang singularity. In the Kasner 
models at t ^ one finds that either two of the k go to zero, whereas the 
third unboundedly increases (a cigar); or one of the h tends to zero, while 
the other two approach a finite value (pancake). Also there is the "barrel" 
singularity in which the two of the h go to zero, and the third approaches 



of the symmetry group to simplify the spatial metric gab (see 31C,311] for more 
details). 
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a finite value. There is an open question as to whether some other possibili- 
ties exist Even in the perfect fluid Kasner model, the approach to the 
singularity is "velocity- dominated" - the "vacuum terms" given by the rates 
of change of the scale factors dominate the "matter terms" (curvature terms 
vanish since the Kasner models are spatially flat). 

The general vacuum Bianchi type II cosmologies (with one n^'^^ = +1, 



and the other two vanishing), discovered by Taub in |197|, contain two free 
parameters: 

= -A^dt^ + A-^ t^P^dx + Apibz dyf + A^{t^P^ dy^ + t^Pa ^^2), (102) 
where 

A^ = l + hh^P\ Pi+P2+P3 = l, pI+pI+pI = 1- (103) 



If we put the parameter b — 0, the metrics (102) become the Kasner solutions 



(100). Near the big bang the general Bianchi type II solution is asymptotic 



to a Kasner model. In the future it is asymptotic again to a Kasner model. 



but with different values of parameters pi (see e.g. [311|). This fact will be 
important in the following. 

The general Bianchi type V vacuum solutions are also known - these are 



given by the 1-parameter family of Joseph solutions |311 . The type V models 
are the simplest metrics in class B (with all n^"-' = but C^^ — 2aj^5°j, at, = 
constant), and are the simplest Bianchi models which contain the standard 
FRW open universes (fc — —1). The Joseph solutions are asymptotic to the 
specific Kasner solution in the past, and tend to the "isotropic Milne model" 
in the future. This is intuitively understandable since open FRW models, 
as they expand indefinitely into the future with matter density decreasing, 
also approach the Milne model. As is well known, the Milne model is just an 
empty flat (Minkowski) spacetime in coordinates adapted to homogeneous 
spacelike hypersurfaces (the "mass hyperboloids"), with expanding normals 
(see e.g. p8[): 

ds^ = -dT^ + T^[{l + p^)-^dp^ + p\de^ + sin^ 9 dip^)] , (104) 

with T = <(1 — u^)^/^, p = u{l — v?), u = rjt < 1, where t, r, 9, ip are standard 
Minkowski (spherical) coordinates. Because of its significance as an asymp- 
totic solution and its simplicity, the Milne model has been us ed fr equently 
in pedagogical expositions of relativistic cosmology (see e.g. p8| , ^14|) as well 
as in cosmological perturbation theory and quantization (see e.g. |?15| and 
references therein). The Milne universe is also an asymptotic state of other 
Bianchi models such as, for example, the intriguing Lukash vacuum type VII^ 
solution ]31fi| |, which can be interpreted as two monochromatic, circularly po- 
larized waves of time-dependent amplitude travelling in opposite directions 
on a FRW background, with flat or negative curvature spacelike sections. As 
was noticed earlier, some indefinitely expanding Bianchi models approach the 



homogeneous plane wave solutions. Barrow and Sonoda |317| studied the fu 



ture asymptotic behaviour of the known Bianchi solutions in detail by using 
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nonlinear stability techniques; and in [311| dynamical system methods were 
used. 

From the late 1960s onwards the greatest amount of work was probably 
devoted to the Bianchi type IX vacuum models, baptized the Mixmaster 



universe^ by Misner [pl2| . Type IX models are the most general class A 
models with all parameters n^"-* = +1. They are the only Bianchi universes 
which recollapse. If a perfect fluid is permitted as the matter source, the 
non-vacuum type IX solutions contain the closed FRW models (fc = +1) 
with space sections having spherical topology. As a Bianchi I space admits a 
group isomorphic with translations in a 3-dimensional Euclidean space, the 
group of type IX spaces is isomorphic to the group of rotations. None of 
the pairs of three Killing vectors commute. A general Bianchi IX vacuum 
solution is not known, but a particular solution is available: the Taub-NUT 
spacetime, or rather, its spatially homogeneous anisotropic region - the Taub 
universe (see Section 7.2). This fact was, for example, employed in an attempt 
to understand the limitations of the minisuperspace methods of quantum 
gravity: by reducing the degrees of freedom to a general Mixmaster universe 
and then further to the Taub universe one can see what such restrictions 
imply |318| . 

The dynamics of general Bianchi cosmologies - and of the Mixmaster 
models in particular - close to t he bi g bang singularity has been approached 
with essentially three methods |31l| : (i) piecewise approximation methods, 
(ii) Hamiltonian methods, and (iii) dynamical syste m metho ds. In the first 



method, used primarily by Russian cosmologists (cf. |29E, 300| ), the evolution 



is considered to be a sequence of periods in which certain terms in the Ein- 
stein equations dominate whereas other terms can be neglected. The Hamil- 



tonian methods appeared first in the "Mixmaster paper" by Misner |312| 



were reviewed by Ryan |266|, and more recently by Uggla in |311]. With the 
Hamiltonian (canonical) approaches, minisuperspace methods entered gen- 
eral relativity (cf. Section 9.2 on midisuperspace for cylindrical waves). In 
this approach, infinitely many degrees of freedom are reduced to a finite 
number: the state of the universe is described by a "particle" moving inside 
and reflecting instantaneously from the moving potential walls, which approx- 
imate the time-dependent potentials in the Hamiltonian. In the third method 
one employs the fact that Einstein's equations in the case of homogeneous 
cosmologies can be put into the form of an autonomous system of first-order 
(ordinary) differential equations on a finite dimensional space This is 
of the form dx/dt — f(x), with x G J?" representing a state of the model 
(for example, the suitably normalized components of the shear a, the Hubble 



The name comes from the fact that, in contrast to a standard FRW model, which 
has a horizon preventing the equalization of possible initial inhomogeneities over 
large scales, the horizon in a type IX universe is ab sent, so that mixing is in prin- 



ciple possible. However, as was shown e.g. in |309|, "repeated circumnavigations 



of the universe by light are impossible in the Mixmaster model" 
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scalar 0, and parameters related to 71^°"^ can serve as the "components" of x). 
A study of the orbits x{t) indicates the behaviour of the model. Dynamical 
system methods are the focus of the book |311|. They also are the main tools 
of the monograph |319|. 

In the case of the Bianchi IX models (either vacuum or with perfect fluid), 
all three methods imply (though do not supply a rigorous proof) that an ap- 
proach to the past big bang singularity is composed of an infinite sequence of 
intervals, in each of which the universe behaves approximately as a specific 
Kasner model (IOC). The transition "regimes" between two different subse- 
quent Kasner epochs, in which the contraction proceeds along subsequently 
differ ent axes, is approximately described by Bianchi type II vacuum solutions 
(102). This famous and enigmatic "oscillatory approach to the singularity" 
(or "Mixmaster behaviour") has rightly entered the classical literature (cf. 
e.g. 1 18, 19 139 ). It indicates that the big bang singularity (and, similarly, a 
singularity formed during a gravitational collapse) can be much more com- 
plicated than the "point-like" singularity in the standard FRW models. This 
oscillatory character has been suggested not only by the qualitative methods 
mentioned above, but also by extensive numerical work (see e.g. |311, 320| |). 
So far, however, it has resisted a rigorous proof. 

In the "standard" picture of the Mixmaster model it is supposed that the 
evolution of the Bianchi type IX universe near the singularity can be approxi- 
mated by a mapping of the so called Kasner circle onto itself. This is the unit 
circle in the S^S^ plane, where S± — a±/H describes the anisotropy in the 
Hubble flow (cf. e.g. Fig. 6.2 in |311|). Each point on the circle corresponds 
to a specific Kasner solution with given fixed values of parameters pi satis- 
fying the conditions (101). There are three exceptional points on the circle 
- those at which one of the pi — +1, and the other two vanish. From each 
non-exceptional point Pi on the Kasner circle there lead s a 1-dimensional 
unstable orbit given by the vacuum Bianchi II solution (102), which joins 
Pi to another point P2 on the circle, then P2 is mapped to P3 , etc. This 
"Kasner map" in the terminology of 311 , called frequently also the BKL 



(Belinsky-Khalatnikov-Lifshitz) map, describes subsequent changes of Kas- 
ner epochs during the oscillatory approach to a singularity. Recent rigorous 
results of Kendall |321| show that for any finite sequence generated by the 
BKL map, there exists a vacuum Bianchi type IX solution which reproduces 
the sequence with any required accuracy.^ 

The vacuum Bianchi IX models have been extensively analyzed in the 
context of deterministic chaos and their stochasticity, attracting the interest 
of leading experts in these fields [320 322 1. Above all, it is the numerical work 
which strongly suggests that it is impossible to make long-time predictions of 



A. Rendall (private communication) reports that the main points of the BKL 
picture for homogeneous universes have been rigorously confirmed in a recent 
work of H. Ringstrom (to be published). 
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the evolution of the system from the initial data, which is the most significant 
property of a chaotic system. 

Most recently, interest in the Bianchi cosmologies with (homogeneous) 
magnetic and scalar fields has been revived. Following his previous work with 
Wainwright and Kerr on magnetic Bianchi VIq cosmologies |32E], LeBlanc 
[ p24| has shown that even in Bianchi I cosmologies one finds an oscillatory 
approach towards the initial singularity if a magnetic field in a general di- 
rection is present. (The points on the Kasner circle are now joined by Rosen 
magneto- vacuum solutions.) Hence, Mixmaster-like oscillations occur due to 
the magnetic field degrees of freedom, even in the absence of an anisotropic 
spatial curvature (present in the vacuum type IX models) - the result an- 
ticipated by Jantzen [|325[ in his detailed work on Hamiltonian methods for 
Bianchi cosmologies wit h m agnetic and scalar fields. Similar conclusions have 
also been arrived at in |326] for magnetic Bianchi II cosmologies. (LeBlanc's 
papers contain some new exact magnetic Bianchi solutions and a number of 
references to previous work.) Interestingly, in contrast to the magnetic field, 
scalar fields in general suppress the Mixmaster oscillations when approaching 
the initial singularity [327, 328 1. 

The theory of spatially homogeneous, anisotropic models is an elegant, in- 
triguing branch of mathematical physics. It has played an important role in 
general relativity. The classical monograph of Zel'dovich and Novikov |309|, 
or the new volume of Wainwright and Ellis |311| analyze in detail the possi- 
ble observational relevance of these models: they point out spacetimes close 
to FRW cosmologies (at least during an epoch of finite duration) which are 
compatible with observational data. For the most recent work on Bianchi 
VII/i cosmologies which are potentially compatible with the highly isotropic 
microwave background radiation, see ]329( | (and references therein). Never- 
theless, the present status is such that, in contrast to for example the Kerr 
solution, which is becoming an increasingly strong attractor for practical as- 
trophysicists (cf. Section 4.3), the anisotropic models have not really entered 
(astro) physical cosmology so far. Peebles, for example, briefly comments in 
[ P8| : "The homogeneous anisotropic solutions allowed by general relativity 
are a very useful tool for the study of departures from the Robertson- Walker 
line element. As a realistic model for our Universe, however, these solutions 
seem to be of limited interest, for they require very special initial conditions: 
if the physics of the early universe allowed appreciable shear, why would it 
not also allow appreciable inhomogeneities?" 

An immediate reaction, of course, would be to point out that the FRW 
models require still more "special initial conditions" . However, there appears 
to be a deeper reason why the oscillatory approach towards a singularity may 
be of fundamental importance. Belinsky, Khalatnikov and Lifshitz 1 299, 300 1 
employed their piecewise approximation method, and concluded 30 years ago 
that a singularity in a general, inhomogeneous cosmological model is spacelike 
and locally oscillatory: i.e. in their scenario, the evolution at different spatial 
point decouples. At each spatial point the universe approaches the singular- 
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ity as a distinct Mixmaster universe. This view, often criticized by purists, 
appears now to be gaining an increasing number of converts, even among the 
most rigorous of relativists. As mentioned above, the homogeneous magnetic 
Bianchi type VIo models, investigated by LeBlanc et al., show Mixmaster 
behaviour. The Bianchi VIq models have, as do all Bianchi models, three 
Killing vectors, but two of them commute. The models can thus be general- 
ized by relaxing the symmetry connected with the third Killing vector; one 
can so obtain effectively the inhomogeneou s (in one dimension) Gowdy-type 



spacetimes. Weaver, Isenberg and Berger [ p30[ , following this idea of Ren 



dall, analyzed these models numerically, and discovered that the Mixmaster 
behaviour is reached at different spatial points. The numerical evidence for 
an oscillatory singularity in a generic vacuum U{1) symmetric cosmologies 
with the spatial topo logy of a 3-torus has been found still more recently by 



Berger and Moncrief 331 



Before turning to the Gowdy models, a last word on the "oscillatory ap- 
proach towards singularity". I heard E. M. Lifshitz giving a talk on this issue 
a couple of times, with Ya. B. Zel'dovich in the audience. In discussions after 
the talk, Zel'dovich, who appreciated much this work (its detailed description 
is included in ||309| ), could not resist pointing out that the number of oscil- 
lations and Kasner epochs will be very limited (to only about ten) because 
of quantum effects which arise when some scale of a model is smaller than 
the Planck length Ipi 10~'^^cm. This, however, seems to make the scenario 
still more intriguing. If this is confirmed rigorously within classical relativity, 
how will a future quantum gravity modify this picture? 



12.2 Inhomogeneous cosmologies 

Amo ng all of the known vacuum inhomogeneous models, the Gowdy solutions 
[ |332[ have undoubtedly played the most distinct role. They belong to the 
class of solutions with two commuting spacelike Killing vectors. Within a 
cosmological context, they form a subclass of a wider class of G2 cosmologies 
— as are now commonly denoted models which admit an Abelian group G2 
of isometrics with orbits being spacelike 2-surfaces. A 2-surface with a 2- 
parameter isometry group must be a space of constant curvature, and since 
neither a 2-sphere nor a 2-hyperboloid possess 2-parameter subgroups, it 
must be intrinsically flat. If the 2-surface is an Euclidean plane or a cylinder, 
then one speaks about planar or cylindrical universes. Gowdy universes are 
compact - the group orbits are 2-tori T^. 

The metrics with two spacelike Killing vectors are often called the general- 
ized Einstein-Rosen metrics as, for example, by Carmeli, Charach and Malin 
[ [333| in their comprehensive survey of inhomogeneous cosmological models 
of this type. In dimensionless coordinates (i, z, x^, x^), the line element can 
be written as {A, B — 1,2) 



(105) 
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where L is a constant length, F and ^ab depend on t and z only, and thus 
the spacelike Killing vectors are ^^^C" = (0, 0, 1, 0), (2)^° = (0, 0, 0, 1). 

The local behaviour of the solutions of this form is described by the gra- 

1 /2 

dient of the "volume element" of the group orbits W = (|det(7AB)|) • 
Classical cylindrical Einstein-Rosen waves (cf. Section 9) are obtained if W^a 
is globally spacelike. In Gowdy models, W,a varies from one region to an- 
other .0 

Considering for simplicity the polarized Gowdy models (when the Killing 
vectors are hypersurface orthogonal), the metric (|105| ) can be written in di- 
agonal form (cf. equations (|62|), (|63|), and (|6^), (|66[) in the analogous cases 



of plane and cylindrical waves) 

in which U{t,z) and 7(i, z) satisfy wavelike dynamical equations and con- 
straints following from the vacuum Einstein equations; the function W{t, z), 
which determines the volume element of the group orbit, can be cast into 
a standard form which depends on the topology of i = constant spacelike 
hypersurfaces S. 

As mentioned above, in Gowdy models one assumes these hypersurfaces 



to be compact. Gowdy |332] has shown that E can topologically be (i) a 
3-torus T'"^ = S'^ ® S'^ ® and W — t (except for the trivial case when 
spacetime is identified as a Minkowski space), (ii) a 3-handle (or hypertorus, 
or "closed wormhole") (8) with W = sin z sin or (iii) a 3-sphere S^, 
again with W — sin z sin t. (For some subtle cases not covered by Gowdy, see 
[ t334| .) As the form of W suggests, in the case of a topology, the universe 
starts with a big bang singularity at t = and then expands indefinitely, 
whereas in the other two cases it starts with a big bang at i = 0, expands 
to some maximal volume, and then recoUapses to a "big crunch" singularity 
at t = TT. One can determine exact solutions for metric functions in all three 



cases in terms of Bessel functions [335|. Hence, for the first time cosmological 



models closed by gravitational waves were constructed. Charach found Gowdy 



universes with some special electromagnetic fields [ 336 1 , and other generalized 
Gowdy models were obtained. We refer to the detailed survey |333| for more 
information, including the work on canonical and quantum treatments of 
these models, done at the beginning of the 1970s by Berger and Misner, and 
for extensive references. 

Let us only add a few remarks on some more recent developments in which 
the Gowdy models have played a role. Gowdy-type models have been used 
to study the propagation and collision of gravitational waves with toroidal 
wavefronts (as mentioned earlier, 2-tori are the group orbits in the Gowdy 

The same is true in the boost-rotation symmetric spacetimes considered in Sec- 
tion 11: the part > of the spacetimes, where the boost Killing vector is 
spacelike, can be divided into four different regi ons, in two of which vector W^a 



is spacelike, and in the other two timelike - see |288] for details. 
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cosmologies) in the FRW closed universes with a stiff fluid |337 . In the stan- 
dard Gowdy spacetimes it is assumed that the "twists" associated with the 
isometry group on vanish. In |33S| the generalized Gowdy models without 



this assumption are considered, and their global time existence is proved. 

As both interesting and non-trivial models, the Gowdy spacetimes have 
recently attracted the attention of mathematical and numerical relativists 
with an increasing intensity, as indicated already at the end of the previous 



section. Chrusciel, Isenberg and Moncrief 1 33£ ] proved that Gowdy spacetimes 
developed from a dense subset in the initial data set cannot be extended past 
their singularities, i.e. in "most" Gowdy models the strong cosmic censorship 
is satisfied. 

On cosmic censorship and spacetime singularities, espec ially in the context 



of compact cosmologies, we refer to a review by Moncrief [34C|, based on his 
lecture in the GR14 conference in Florence in 1995. The review shows clearly 
how intuition gained from such solutions as the Gowdy models or the Taub- 
NUT spaces, when combined with new mathematical ideas and techniques, 
can produce rigorous results with a generality out of reach until recently. To 
such results belongs also the very recent work of Kichenassamy and Rendall 
[ ^4l| on the sufficiently general class of solutions (containing the maximum 
number of arbitrary functions) representing unpolarized Gowdy spacetimes. 



The new mathematical technique, developed by Kichenassamy [342|, the so 
called Fuchsian algorithm, enables one to construct singular (and nonsingu- 
lar) solutions of partial differential equations with a large number of arbitrary 
functions, and thus provide a description of singularities. Applying the Fuch- 
sian algorithm to Einstein's equations for Gowdy spacetimes with topology 
T^, Kichenassamy and Rendall have proved that general solutions behave at 
the (past) singularity in a Kasner-like manner, i.e. they are asymptotically 
velocity dominated with a diverging Kretschmann (curvature) invariant. One 
needs an additional magnetic field not aligned with the two Killing vectors of 
the Gowdy unpolarized spacetimes in order to get a general oscillatory (Mix- 
master) approach to a singularity, as shown by the numerical calculations 
[ |330| mentioned at the end of the previous section. 

Much of the work on exact inhomogeneous vacuum cosmological mod- 



els has been related to "large perturbations" of Bianchi universes. In [343 
the authors confined attention to "plane wave" solutions propagating over 
Bianchi backgrounds of types I- VII. They found universes which are highly 
inhomogeneous and "chaotic" at early times, but are transformed into clearly 
"recognizable" gravitational waves at late times. 

Other types of metrics can be considered as exact "gravitational solitons" 
propagating on a cosmological background. These are usually obtained by ap- 
plying the inverse scattering or "soliton" technique of Belinsky and Zakharov 
[B43 to particular solutions of Einstein's equations as "seeds" . For example. 



Carr and Verdaguer |345| found gravisolitons by applying the technique to 



the homogeneous Kasner seed. Similarly to previous work |343 , their models 
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are very inhomogeneous at early times, but evolve towards homogeneity in a 
wavelike manner at late times. 

More recently, Belinsky |346|, by applying a two-soliton inverse scattering 
technique to a Bianchi type VIq solution as a seed, constructed an intriguing 
solution which he christened as a "gravitational breather", in analogy with the 
Gordon breather in the soliton theory of the sine-Gordon equation. Gravisoli- 
tons and antigravisolitons, characterized by an opposite topological charge, 
can be heuristically introduced and shown to have an attractive interaction. 
The breather is a bound state of the gravisoliton and antigravisoliton. Be- 
linsky suggests that a time oscillating breather exists; but a later discussion 
[ ^47| indicates that the oscillations quickly decay. Alekseev, by employing his 
generalization of the inverse scattering method to the Einstein-Maxwell the- 
ory, obtained e xact electrovacuum solutions generalizing Belinsky 's breather 
(see his review |34^ , containing a general introduction on exact solutions). 

Verdaguer |34£] prepared a very complete review of solitonic solutions 
admitting two spacelike Killing vector fields, with the main emphasis on cos- 
mological models. Among various aspects of such solutions, he has noted 
the role of the Bel- Robinson superenergy tensor in the interpretation of cos- 
mological metrics. This tensor and its higher-order generalizations has also 
been significa ntly used in estimates in the proofs of long-time existence theo- 
rems |39| , |340[ . Recently, differential conservation laws for large perturbations 
of gravitational field with respect to a given curved background have been 
fomulated |35C], which found an application in solving eq uatio ns for cosmo- 
logical perturbations corresponding to topological defects [351|. They should 
bring more light also on various solitonic models in cosmology. 



13 Concluding remarks 

It is hoped that the preceding pages have helped to elucidate at least one 
issue: that in such a complicated nonlinear theory as general relativity, it is 
not possible to ask relevant questions of a general character without finding 
and thoroughly analyzing specific exact solutions of its field equations. The 
role of some of the solutions in our understanding of gravity and the universe 
has been so many-sided that to exhibit this role properly on even more than 
a hundred pages is not really feasible ... 

Although we have concentrated on only (electro) vacuum solutions, there 
remains a number of such solutions that have also played some role in various 
contexts, but, owing to the absence of additional space and time, or the 
presence of the author's ignorance, have not been discussed. Tomimatsu-Sato 
solutions and their generalizations, static plane and cylindrical metrics, and 
some algebraically special solutions are examples. 

In his review of exact solutions, Ehlers wrote 35 years ago that "it 
seems desirable to construct material sources for vacuum solutions" , and 30 
years later Bonnor 164], in his review, expressed a similar view. In the above 
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we have noted only some of the thin disk sources of static and stationary 
spacetimes in Section 6. To find physically reasonable material sources for 
many of the known vacuum solutions remains a difficult open task. In or- 
der to make solutions of Einstein's equations with the right-hand side more 
tractable, one is often tempted to sacrifice realism and consider materials, 
again using Bondi's phraseology, which are not easy to buy in the shops. 
Nevertheless, there are solutions representing spacetimes filled with matter 
which would certainly belong in a more complete discussion of the role of 
exact solutions. 

For example, one of the simplest, the spherically symmetric Schwarzschild 
interior solution with an incompressible fluid as matter source, modelling 
"a star of uniform density" , gives surprisingly good estimates of an upper 
bound on the masses of neutron stars; on a more general level, it supplies 
an instructive example of relativistic hydrostatics [[l8[ . Many other spherical 
perfect fluid solutions are listed in j6^. The proof of a very plausible fact 
that any equilibrium, isolated stellar model whi ch is nonrotating must be 



spherically symmetric, was finally completed in 352 and 353 . Physically 



more adequate spherically symmetric static solutions with collisionless mat- 
ter described by the Boltzmann (Vlasov) equation have been studied [101| 
(yielding, for example, arbitrarily large central redshifts); and some of their 
aspects have been recently reviewed from a rigorous, mathematical point of 
view [354|. Going over to the description of matter in terms of physical fields, 
we should mention the first spherically symmetric regular solutions of the 
Einstein- Yang-Mills equations ( "non- Abelian solitons" discovered by Bart- 
nik and McKinnon |355| in 1988), and non- Abelian black holes with "hair", 
which were found soon afterwards. They stimulated a remarkable activity in 
the search for models in which gravity is coupled with Yang-Mills, Higgs, and 
Skyrmion fields. Very recently these solutions have been surveyed in detail 
in the review by Volkov and Gal'tsov | ]356{ . 

The role of the standard FRW cosmological models on the development 
of relativity and cosmology can hardly be overemphasized. As for two more 
recent examples of this influence let us just recall that the existence of cosmo- 
logical horizons in these models was one of the crucial points which inspired 
the birth of inflationary cosmology (see e.g. ||2^) ; an d the very smooth char- 
acter of the initial singularity has led Penrose |102] to formulate his Weyl 
curvature hypothesis, related to a still unclear concept of gravitational en- 
tropy. Homogeneous but anisotropic Bianchi models filled with perfect fluid 
are extensively analyzed in |311]. Very recent studies of Bianchi models with 
collisionless matter [|357| reveal how the matter content can qualitatively alter 
the character of the model. 

A number of Bianchi models approach self-similar solutions. Perfect fluid 
solutions admitting a homothetic vector, which in this case implies both 
geometrical and physical self-similarity, have been reviewed most recently by 



Carr and Coley |358|. In their review various astrophysical and cosmological 



applications of such solutions are also discussed. Self-similar solutions have 
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played a crucial role in the critical phenomena in gravitational collapse. Since 
their discovery by Choptuik in 1993, they have attracted much effort, which 
has revealed quite unexpected facts. In | |358[ these phenomena are analyzed 
briefly. For a more comprehensive review, see [359|. 

Self-similar, spherically symmetric solutions have been very relevant in 
constructing exam ples of the formation of naked singularities in gravitational 
collapse (see 358 for a brief summary and references). In particular, the in- 
homogeneous, spherically symmetric Lemaitre-Bondi-Tolman universes con- 
taining dust have been employed in this context. Solutions with null dust 
should be mentioned as well, especially the spherically symmetric Vaidya so- 
lutions: imploding spherical null-dust models have been constructed in which 
naked singularities arise at their centre (see |]3^ for summary and references). 

The Lemaitre-Bondi-Tolman models are the most frequently analyzed in- 
homogeneous cosmological models which contain the standard FRW dust 
models as special cases (see e.g. [ p8[|3^ ). In his recent book Krasihski [360| 
has compiled and discussed most if not all of these exact inhomogeneous cos- 
mological solutions found so far which can be viewed as "exact perturbations" 
of the FRW models. 

Many solutions known already still wait for their role to be uncovered. 
The role of many others may forever remain just in their "being" . However, 
even if new solutions of a "Kerr-like significance" will not be obtained in 
the near future, we believe that one should not cease in embarking upon 
journeys for finding them, and perhaps even more importantly, for revealing 
new roles of solutions already known. The roads may not be easy, but with 
todays equipment like Maple or Mathematica, the speed is increasing. Is 
there another so explicit way of how to learn more about the rich possibilities 
embodied in Einstein's field equations? 

The most remarkable figure of Czech symbolism, Otokar Bfezina (1868- 
1929) has consoling words for those who do not meet the "Kerr-type" metric 
on the road: "Nothing is lost in the world of the spirit; even a stone thrown 
away may find its place in the hands of a builder, and a house in flames may 
save the life of someone who has lost his way...". 
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